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This article aims at discussing the cosmological constant problem at a pedagogical but fully 
technical level. We review how the vacuum energy can be regularized in flat and curved space- 
time and how it can be understood in terms of Feynman bubble diagrams. In particular, we show 
that the properly renormalized value of the zero-point energy density today (for a free theory) is 
in fact far from being 122 orders of magnitude larger than the critical energy density, as often 
quoted in the literature. We mainly consider the case of scalar fields but also treat the cases of 
fermions and gauge bosons which allows us to discuss the question of vacuum energy in super- 
symmetry. Then, we discuss how the cosmological constant can be measured in cosmology and 
constrained with experiments such as measurements of planet orbits in our solar system or atomic 
spectra. We also review why the Lamb shift and the Casimir effect seem to indicate that the 
quantum zero-point fluctuations are not an artifact of the quantum field theory formalism. We 
investigate how experiments on the universality of free fall can constrain the gravitational properties 
of vacuum energy and we discuss the status of the weak equivalence principle in quantum mechanics, 
in particular the Collela, Overhausser and Werner experiment and the quantum Galileo experiment 
performed with a Salecker-Wigner-Peres clock. Finally, we briefly conclude with a discussion on the 
solutions to the cosmological constant problem that have been proposed so far. 
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I. INTRODUCTION 

This review article aims at discussing the cosmolog- 
ical constant problem [1—15]. This question is central 
in modern physics because its resolution would certainly 
mean a very important step forward in our attempts to 
understand physics beyond the current standard model. 

The history of the cosmological constant problem is a 
long and rich one [5], Its premises were in fact already 
present immediately after the birth of quantum field the- 
ory. Indeed, from the Heisenberg uncertainty principle, 
we know that the ground state energy of the quantum 
mechanical oscillator cannot be zero because the poten- 
tial and the kinetic energies cannot vanish at the same 
time. Since a (free) quantum field can be viewed as an 
infinite collection of harmonic oscillators, it immediately 
comes that its ground state energy must be infinite. Of 
course, this is not yet the cosmological constant problem 
because gravity does not enter the stage. But, clearly, 
potentially, a very severe problem is already present. 

The infinity mentioned above is the first infinity that 
one encounters in quantum field theory [16-22]. How- 
ever, the presence of this infinity did not prevent the 
founding fathers of quantum field theory to develop the 
theory since, in absence of gravity (a working assump- 
tion of quantum field theory justified by the weakness 
of gravity), only differences in energies arc observable. 
Therefore, in this context, this infinity can be merely ig- 
nored and the rest of the formalism can be worked out 
without any problem. As is well-known, other types of 
infinities appear and, this time, there is no way to get 
around them. Treating and taming them is the goal of 
rcnormalization. The very impressive agreement between 



high precision measurements in accelerators and the pre- 
dictions of quantum field theory in presence of radiative 
corrections is the proof that renormalization is able to 
correctly regulate these infinities. 

Therefore, we see that these two types of infinities are 
treated very differently. The problem of the zero-point 
energy is just avoided while the problem of the radiative 
corrections is directly and explicitly addressed. Clearly, 
one cannot help thinking that the first problem is in fact 
swept under the carpet. And, indeed, as soon as gravity 
is turned on, it badly strikes back. 

At this point, it is worth noticing the following. One 
should not get the impression that the zero-point en- 
ergy cannot be rcnormalized. As a matter of fact, as 
will be discussed in this review, it can be made perfectly 
finite. However, this finite, renormalized, value of the 
zero-point energy seems to be too large to be compati- 
ble with the observations. Again, this is very different 
from the usual case of quantum field theory where the 
finite part extracted from a divergent expression always 
leads to a good agreement with the experiments. At first 
sight, the cosmological constant problem is therefore nei- 
ther the presence of a new infinity nor our inability to 
regularize it but rather the apparent failure of the renor- 
malization scheme to produce, at the quantitative level, 
a finite vacuum energy compatible with the observational 
data. 

In fact, a failure of renormalization is not the only 
logical possibility. The vacuum fluctuations could also 
be a mathematical artifact of the quantum field theory 
framework and/or their gravitational properties could be 
abnormal. However, when one tries to test (theoretically 
or experimentally) these ideas, it seems that one never 
encounters a problem of the kind mentioned above. This 
fact must also be considered as a part of the cosmological 
constant problem which, therefore, appears to have many 
different ramifications. We see that this problem is in fact 
a very deep problem which lies at the cross roads between 
different branches of physics (gravitational physics, quan- 
tum field theory, cosmology etc . . . ) . In brief, it has to 
do with the gravitational properties of the quantum vac- 
uum. 

The importance of this question has recently been re- 
inforced by the discovery that the expansion of our Uni- 
verse is accelerated [23, 24]. According to the standard 
lore, this could be the first observational evidence that 
the quantum vacuum is able to curve space-time. On the 
other hand, the effect of a cosmological constant is visible 
only on large scales and, therefore, it seems to be prob- 
lematic to check this result elsewhere than in cosmology. 
But, if true, the cosmological constant problem repre- 
sents a unique situation where some aspects of quantum 
gravity are at play and where, at the same time, cor- 
responding observational signatures are not hopelessly 
beyond our technical capabilities. This makes it a valu- 
able opportunity to go beyond our current understanding 
of theoretical physics. The only other situation where 
the three above mentioned aspects (quantum mechan- 
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ics, gravity and the possibility to detect sizable observa- 
tional effects) are mixed is the theory of cosmological in- 
flation [25-34] (for a review, see Refs. [35-37]) where the 
Cosmic Microwave Background (CMB) anisotropics, first 
detected by the COBE (COsmic Background Explorer) 
satellite [38], arc supposed to originate from the quantum 
fluctuations of the inflation and gravitational fields (for 
the observational status of inflation, see Ref. [39]). 

In this review article, we go through all the issues men- 
tioned before. Our prejudice is to present the question at 
the technical level, in full details, in order for the paper 
to be reasonably self-contained. The price to pay for this 
approach is that, sometimes, we are led to consider prob- 
lems that are well treated in standard textbooks and/or 
are not, strictly speaking, directly concerned with the 
cosmological constant problem. On the other hand, this 
is also an opportunity to review many different subjects 
belonging to different types of physics. 

The outline of the article is the following one. In 
Sec. II, we discuss how the cosmological constant A is 
introduced in the Einstein field equations and we explain 
why the vacuum energy also participates in the value of 
A. In Sec. Ill, we treat the classical cosmological con- 
stant problem, i.e. we show that, in the presence of grav- 
ity, it is no longer true that only differences of energy 
are observable. We illustrate this discussion by means of 
the electro-weak phase transition. In Sec. IV, we demon- 
strate that, even if the classical ground state is tuned 
to zero, the quantum vacuum fluctuations still give a 
contribution to the cosmological constant. We discuss 
how the corresponding quantum vacuum energy can be 
calculated and argue that the method often used in the 
literature and which consists in introducing a cut-off, is 
not appropriate. We also investigate how the quantum 
vacuum energy can be expressed in terms of Feynman di- 
agrams (bubble diagrams) and investigate the case where 
interactions are present. In Sec. V, we focus on the so- 
called bubble diagrams. We discuss them in the context 
of quantum field theory but also in ordinary quantum 
mechanics. We then show that they have very different 
properties than the other loop diagrams. In Sec. VI, we 
use the Gaussian effective potential approach to calcu- 
late the vacuum energy in a situation where non per- 
turbative effects are present. Since all the calculations 
are carried out in the case of a scalar field, we treat in 
Sec. VII the case of other types of fields (spinor and vec- 
tor fields). This also gives us the opportunity to discuss 
the idea of super-symmetry. In Sec. VIII, we present the 
cosmological constant problem in a more rigorous way. 
In particular, we show that a flat space-time calculation 
is a good approximation although a fully consistent ap- 
proach should be formulated in curved space-time. One 
concludes the first part of this review article with Sec. IX 
where we estimate the vacuum energy and find a value 
very far from the often quoted "122 orders of magni- 
tude" . In Sec. X, we explain, from a theoretical point 
of view, how the value of the cosmological constant was 
recently measured in cosmology. In particular, we dis- 



cuss the hypothesis that arc implicitly assumed in order 
to obtain this result. In Sec. XI, wc discuss whether the 
cosmological constant can be measured elsewhere than 
in cosmology and argue that no other experimental con- 
text can compete with cosmology. However, we also show 
that other experiments can put constraints on the vac- 
uum energy that are interesting from the point of view 
of the cosmological constant problem. In Sec. XII, we 
review the experiments (Lamb shift, Casimir effect) that 
seem to prove the existence of the vacuum fluctuations. 
In Sec. XIII, we investigate whether there are observa- 
tions that can probe their gravitational properties. In 
Sec. XIV, we discuss whether the weak equivalence prin- 
ciple still holds in quantum mechanics since this has ob- 
vious implications for the question of the vacuum weight. 
Finally, in Sec. XV, wc present our conclusions and very 
briefly review the solutions to the cosmological constant 
problem that have been proposed so far. Let us mention 
that this article is part of a more general review on the 
dark energy question, see Refs. [40-43]. 

Before concluding this introduction, we would like to 
make the following remark. This article aims at dis- 
cussing what the cosmological constant problem is. This 
article does not aim at discussing what solutions to this 
problem have been proposed even if, as already men- 
tioned, in Sec. XV, we say a few words on this topic 
(and choose, in a totally arbitrary way which only reflects 
the author prejudices, to focus more on some of them). 
The justification for the above mentioned choice is that 
the question of the gravitational properties of vacuum 
zero-point fluctuations is already a highly non-trivial one 
as the variety of the subjects listed above demonstrate. 
Therefore, it seems to us better to deeply understand 
what is at play before trying to propose a way out. Of 
course, it is because this problem is very difficult and has 
many different ramifications in many different branches 
of physics that it is so interesting. 

II. THE COSMOLOGICAL CONSTANT 

In this first section, we introduce the cosmological con- 
stant as a free parameter in the classical action of the 
gravitational field. This action, together with the action 
describing matter, can be expressed as 

S = ^- I d i X^](R~2A B ) + Smarter *] , (1) 

2k J 

where K = 8ttG/c 4 = 87r/m P1 = 1/Mp\, m P , and M P1 be- 
ing the Planck mass and the reduced Planck mass respec- 
tively The first term is the standard Einstein-Hilbert 
action and represents the gravitational part of the total 
action given by Eq. (1). The cosmological constant ap- 
pears in the second term of the above expression. As 
announced, at this level, it is merely a new parameter 
of the total action and for this reason we write it as A B 
for "bare cosmological constant". It has the dimension 
of the inverse of a square length. It is compatible with 
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general covariance and is of course compatible with a 
conserved energy momentum tensor since V^g^ = 0. 
Therefore, this term appears to be totally natural from 
the relativistic point of view and there is a priori no rea- 
son to discard it. Since, according to the standard lore 
of field theory, everything which is not forbidden should 
be considered, the cosmological constant should clearly 
be included in our description of the gravitational field. 
Finally, the third term in the above equation denotes the 
matter action where \1/ represents a generic matter field 
that we do not need to specify at this stage. Variation of 
the total action with respect to the metric tensor leads 
to the Einstein equations of motion which read 



where the stress-energy tensor is defined by 



2 5S a 



5g^ v 



(2) 



(3) 



As already mentioned, A B appears to be just a parameter 
of the model. At this level, the only thing that can be 
done is to try to constrain it using various observations. 
But, clearly, in the present framework, there is no way to 
calculate its value from more fundamental considerations. 

However, as originally shown by Sakharov [44] , the na- 
ture of the discussion is crucially changed when one takes 
into account quantum field theory. The point is that the 
stress energy tensor of a field placed in the vacuum state 
must be given by 



(4) 



where p vac is the constant energy density of the vacuum. 
This equation is valid for all the fields present in the 
Universe, as will be shown below. There are several ways 
to prove the above result. Firstly, one can use the fact 
that, in flat space-time (in Minkowski space-time), the 
only invariant tensor is rj^ . Since the vacuum state must 
be the same for all observers, one necessarily has {T^ v ) oc 
t]^. In curved space-time, this means that 



(Tpv) = -p vac (£,x) 5/J 



(5) 



and from the fact that the stress-energy tensor must be 
conserved, we reach the conclusion that p vac must be 
a constant. Therefore, one obtains the formula shown 
above. 

Another way to obtain the same result is to consider 
a specific example, for instance a scalar field $. The 
corresponding action reads 



(6) 



where U($) is the potential. Using the definition (3), the 
corresponding stress-energy tensor can be written explic- 
itly as 



(7) 



From this expression, one sees that the scalar field is 
in fact a perfect fluid. Now, the vacuum state is the 
minimum energy state. Clearly, in order to minimize the 
energy we have to consider a situation where the kinetic 
energy vanishes and where the field sits at the minimum 
of its potential. In this case, the stress energy tensor 
reduces to 



(8) 



which has exactly the expected form with p vac = V($ min ). 

In fact what was shown before is that there are at least 
two sources for the vacuum energy. There is a "classical" 
contribution given by Eq. (8) which originates from the 
value of the potential at its minimum. The corresponding 
vacuum energy will be calculated in Sec. III. There is also 
a "quantum-mechanical" source given by Eq. (5) which 
originates from the zero point fluctuations of the ground 
state. This problem will be treated in Sec. IV. 

Having established the form of the stress energy mo- 
mentum tensor in the vacuum, one can now proceed with 
the Sakharov argument. The next step consists in assum- 
ing that the equivalence principle applies to the zero- 
point fluctuations (here, and from now on, we mean to 
the "classical" and "quantum-mechanical" contributions 
as discussed in the last paragraph), that is to say that 
the zero-point fluctuations gravitate. Clearly, this is not 
a trivial step but, after all, the vacuum fluctuations are 
just a specific type of energy and, in general relativity, all 
forms of energy gravitate. Therefore, a consistent way of 
writing the Einstein equations when quantum field the- 
ory is taken into account seems to be 

R» v - \Rg» v + K9»» = «r M T ttcr + k(t^) , (9) 

where on the right hand side the first contribution comes 
from ordinary matter while the second one represents the 
contribution originating from the vacuum (of course, it 
is to be understood that we sum up the contributions 
coming form all the fields present in the universe) . Using 
the form of (T^ v ) established before, one arrives at 



f 



where 



A 



A c 



(10) 



(11) 



Therefore, we conclude that the effective cosmological 
constant is the sum of the bare cosmological constant 
and of a contribution originating from the vacuum fluc- 
tuations. The effective cosmological constant A cff is the 
quantity that one can observe and constrain when tests 
of the Einstein equations are carried out. 

As we now discuss, the problem is that ft/0 vac is made 
of several terms which are all huge in comparison with 
the observed value of A « . 
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III. THE CLASSICAL COSMOLOGICAL 
CONSTANT PROBLEM 

The cosmological constant problem is in fact a multi 
facets question as the quantity A off receives contributions 
from different origins. As was demonstrated in the last 
section on the example of a scalar field, a classical con- 
tribution to A off comes from the value of the potential 
at its minimum. If this is not zero, then it affects the 
value of the effective cosmological constant. This point 
is discussed in this section and is exemplified with the 
electroweak transition. Another contribution comes from 
the quantum zero point fluctuations but this will be dis- 
cussed in the next sections. 



A. Phase Transition 

It is possible to study the dynamics of a phase tran- 
sition with the help of the following simple model. Let 
us consider a scalar field $ in interaction with another 
scalar field $ such that 

V(<$>, *) = V(<5>) + |$ 2 * 2 , (12) 

where g is a dimensionless coupling constant and where 
the self-interacting potential is given by 

V(*) = V + ± ($ 2 - v 2 ) 2 , (13) 

A being a coupling constant describing the self- 
interaction and v the value of the scalar field at the mini- 
mum of its potential in absence of an interaction with ty. 
The quantity Vq denotes the classical off-set. If the field 
^ is in thermal equilibrium, then one is entitled to replace 
^> 2 with ( v I /2 )r where the average is taken in a thermal 
state with temperature T. Since (^ 2 )t oc T 2 [45], the 
effective potential becomes 

V cS ($) = V + ± (<S> 2 - v 2 ) 2 + |T 2 $ 2 , (14) 

where we have slightly redefined the coupling constant g 
in order to take into account the proportionality constant 
between the thermal average and the temperature. This 
potential can also be expressed as 

V cS m = V a + + ^ (g - l) $ 2 + V, (15) 

where we have defined T cr i = v\J X/g. One sees on the 
above equation that the interaction with the field ^ gives 
a temperature dependence to the effective mass of the $ 
field, namely 

m 2 ff (T) = Xv 2 (£r-l). (16) 

\ cri / 

As a consequence, when T > T cri , that is to say before 
the transition, the square of the effective mass is positive 



while, after the transition, when T < T CI1 , it becomes 
negative. In the first situation, the minimum is located 
at $ = and the corresponding value of the vacuum en- 
ergy is Vo + Au 4 /4, see Fig. 1. In the second situation, the 
minimum is located at $ = v and the off-set is simply 
given by Vq, see Fig. 1. We now clearly see the prob- 
lem. If we require the vacuum energy to vanish before 
the transition, we must choose Vq = — Aw 4 /4. But, then, 
after the transition, the vacuum energy is no longer zero 
and is given by the negative value p vac = —Xv /A. On 
the other hand, one can choose p vac to vanish after the 
transition. In this case, one must choose Vb = 0. But 
this means that it was not zero before the transition, 
Pvac = Xv A /A. This last option can maybe viewed as the 
preferred solution since we do not have direct observa- 
tional constraints on the matter content of the universe 
prior to Big Bang Nucleosynthesis (BBN). Anyhow, we 
see that the choice of a parameter in the potential allows 
us to change the vacuum energy at the classical level, 
hence the name "classical cosmological constant prob- 
lem" . But, clearly, the off-set cannot be zero before and 
after the transition. Of course, in order to really establish 
that there is a problem, the question remains to estimate 
this vacuum energy for realistic phase transition and to 
compare it with the observational constraints. In the 
next subsection, we turn to the first issue and consider 
the case of the electroweak phase transition. 



B. The Electro- weak Phase Transition 

In this section, we calculate the vacuum energy "in- 
duced" by the electroweak phase transition. In order 
to achieve this goal, we first must recall some basics 
facts about the standard model of particle physics, see 
Refs. [16-22]. 

The Higgs field is a doublet of complex scalar fields. 
It is charged under the group U(1) Y x SU(2) L character- 
ized by two coupling constant g' and g. The Higgs field 
Lagrangian reads 

£Higg S = - 0D M £) f D^S-F(S,E t ), (17) 
where the covariant derivative can be expressed as 

Z^E = d M E + i|-Y H B„£ + igT a W^, (18) 

where a runs from one to three. The quantities and 
W° are the gauge bosons, Y H is the Higgs weak hyper- 
charge and T a arc the generators of the SU(2) group (see 
below for more details about these quantities). The po- 
tential is chosen to be 

^Et^^StE + ^stE) 2 , (19) 

and clearly resembles Eqs. (14) and (15) for T < T cri . 
Indeed, in this expression, the quantity m 2 is negative 







J I I— I I I— I I l_l I I I I I I I I I I I I I I I I I I u 



-I I— I I I— I I I— I I I— I I I— I I I— I I I— I I I— I I I— I I u 



-3 -2 



-1 





<0> 



FIG. 1: The effective potential given by Eq. (14). Before the transition, for T > T cr i, the minimum of the potential is located 
at the origin and the vacuum energy is given by Vb + Ai> 4 /4. After the transition, for T < T cr i, the minimum is located at $ = v 
and the corresponding vacuum energy has changed and now equals Vb- It is clear that Vb can always be chosen such that the 
vacuum energy vanishes either before or after the transition. It is equally clear that one cannot choose the parameters of the 
potential such that p vac is zero before and after the phase transition. 



(as explained before, m 2 is in fact the square of the ef- 
fective mass) . This also means that the minimum of the 
potential is not located at (E) = but at 

(2)=* = ^ (20) 

obtained from the condition dV/d (E^E) = 0. The quan- 
tity v is in fact the new vacuum expectation value of the 
Higgs after spontaneous symmetry breaking. In particle 
physics, one studies the theory after the transition, when 
the Higgs field has stabilized at its new minimum. But, 
in reality, as explained before, we must consider that the 
transition is a dynamical process. Before the transition 
we had m 2 > and the minimum was located at (E) = 
and after the transition, m 2 has become negative and, 
consequently, the new vacuum is given by v. Returning 
to the minimum after the transition, it is easy to show 
that the value of the potential at (E) = v is given by 

v(m=v) = ~, (2i) 

and is negative. Clearly, this is because we have chosen 
Vq = in Eq. (19). As already explained in Sec. Ill A, 



this means that the vacuum energy vanishes before the 
transition and becomes negative after, see also Fig. 1. 
This means that we have considered the situation corre- 
sponding to the left panel in Fig. 2. Of course, this is 
arbitrary and we could also have considered the situa- 
tion corresponding to the right panel in Fig. 2. Again, 
as explained in the previous subsection, in this case, the 
vacuum energy would vanish after the transition but not 
before. 

The numerical value of the vacuum energy given by 
Eq. (21) is fixed by the electrowcak physics in particu- 
lar by the parameters m and A. Therefore, in order to 
compute p vac explicitly, we must explain how m and A 
are evaluated. For this purpose, we now briefly recall the 
main features of the standard model of particle physics. 
Let us come back to the Higgs Lagrangian given at the 
beginning of this section, see Eq. (17). In the expression 
of ^Higgsj the quantities T a are the generator of SU(2)l 
and are given by T a = a a /2 where a a are the Pauli ma- 
trices 

- -(;;)■•»- (■ D ■ *>= (j " ) ■ m 

Using the gauge invariance, one can always rewrite the 



-3-2-10 1 2 3 

(H) 



-3-2-10 1 2 3 



FIG. 2: Effective potential of the Higgs boson before and after the electroweak phase transition. The left panel corresponds to 
a situation where the vacuum energy vanishes at high temperature. As a consequence p vac is negative at temperature smaller 
than the critical temperature. This is the situation treated in the text where the quantity — m 4 /(4A) is explicitly calculated. On 
the right panel, the off-set parameter Vb is chosen such that the vacuum energy is zero after the transition. As a consequence, 
it does not vanish at high temperatures. 



complex doublet as 



H 



where H is a real scalar field, the "Higgs boson" . There- 
fore, the covariant derivative (18) can be expressed as 



(23) 



9 -w 3 



and the Higgs Lagrangian (17) becomes 

£ H i ggs = -\d^Hd»H - V (£, £*) - ( v 2 + V2vH + ^H 2 ) f^-W^W 1 ^ + ^-W 2 W 2 » + ^-W*W 3 » 



S-Y 2 B»B»- 9 -i-Y n B»W 3 »), 



where the scalar potential can be written as 

/ +s Aw 4 loo Aw o A A 
V (E, Et) = -— + -Xv H + - T H* + -H\ 

We now introduce the physical gauge bosons. As is well-known, they arc defined by 

W+ - W, 



B„ 



sin OZf, + cos OA^ , Wl = cos BZ^ + sin 6A^ , W 2 = M M , W 1 

i\2 



W+ + w- 



V2 



(24) 



(25) 



(26) 



(27) 



where 9 is the Weinberg angle, represents the photon and Z^, the massive gauge bosons. Expressing the 
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Higgs Lagrangian (25) in terms of these new fields, one obtains the following expression 



1 



£ H ig g s = -faHd^H - V (£, E*) - 

„2 



<2vH 



1 



-H' 



9 2, 
— COS I 

4 



9S 



Y„ sin 9 cos ( 



y —w+w~^ 

4 r 



4 



9 . 2 
— sin i 
4 



/2 / 

■9 T.r2 2 /1 99 ■ n 

— Y cos Y„ sin 6 cos i 

4 2 



( y cos6»sin( 



o' 2 

— Y 2 cosflsinf 
2 



35' 



Y„ cos 



33 ^ • 2 /) 
— F H sin (9 



(28) 



Of course £mggs * s n °t the only term in the standard model Lagrangian. There are also the terms describing the 
kinetic terms of the gauge fields 



Ck- 



Gauge 



= --B, 1V B^ 



(29) 



where is the B field strength, namely B^ v = d^B u — 
dyBp and where a similar definition applies to the other 
gauge fields. 

Finally, one must also specify how leptons are de- 
scribed. In the standard model of particle physics, a 
fermion is represented by a four-component Dirac spinor 
^f, see also Sec. VII A where the zero-point fluctuations 
of spinors are calculated. Let us now recall what chiral- 
ity is. For this purpose, we introduce the matrix 7 s such 



that ( 7 5 ) = I and ( 7 5 ) 
two chiral projectors by 



t 



7 



5 . This allows us to define 



P» = 



h 



'75 



P, = 



I4 - 75 



(30) 



where I4 denotes the identity four by four matrix. We 
now define the Dirac matrices, see Sec. VII A for a more 
general approach. In the chiral representation, they are 
given by 



-I a 
-I 2 

I 2 

-I; 



a 1 
-a 1 



(31) 



and the two projectors can be written as 



Given a Dirac spinor \& (we denote the four-component 
spinors by a capital Greek letter and the two-component 
spinors by an ordinary Greek letter) 



(33) 



one can thus define the two states of chirality (right and 
left handed) by means of the two following expressions 



* = ( V R 



* T = R f 



(34) 



The quantities ip n L are two-component spinors known 
as Weyl spinors. They play a very important role in 
super-symmetry, see Sec. VII D. Notice that we also have 
\]> = \j/t-j,o _ _ (^t^tj. This implies that the (mass- 
less) Dirac Lagrangian can be expressed as 



C 



Dii 



(35) 

where we have defined the covariant Pauli matrices = 
(I 2 ,ct 4 ) and a* = (la.-tr*). 

After these brief reminders, we now consider the first 
family of the standard model leptons (the other two fam- 
ilies can be treated in the same fashion and the quarks 
must be considered separately), i.e. the electron and the 
neutrino. Following the usual convention, we denote by 
e R the right-handed electron Weyl spinor and by e L the 
left-handed one. With regards to the neutrino, only the 
left-handed particle is present in the standard model and 
will be denoted as . Then, we have the following prop- 
erties. The particles e L and are charged under U(1) Y 
but also under SU(2) L , i.e. they are considered as a com- 
plex doublet of two-component spinors 



(36) 



On the other hand, e R is charged under U(1) Y but is 
postulated to be a SU(2) L singlet. Therefore, concretely, 
the Lagrangian can be written as 



-'Leptons 



-ih\o» + ijY L B^ + igT a W^j L e 



id^ + i^B^ j 



(37) 



where L\ means (^,e[). Notice the factor 1/2 in front 
of the gauge boson B^ which is the usual convention [this 
factor was already present in Eq. (18)]. The quantities 
Y L and Y R are the leptonic weak hyper-charges. Let us 
focus on the first term in the above expression. The term 
between the parenthesis is a two by two matrix that can 
be written in terms of the physical gauge bosons, using 



9 



the definitions (27). Similarly, the second term can be Lagrangian for the leptons can be rewritten as 
expressed in terms of the fields and Z^. Then, the 



C 



Leptons 



[ |-y L sin6»- | cos 6> ) 



I —Y^ cos( 



?sin<9 1 A, 



i \ —Y, sin 9 H — cos( 
2 L 2 



( |-Y L cos0- | sine 



(38) 



Therefore, the total Lagrangian is now given by £ = 
£mgg S + ^K-Gaugo + ^Leptons- This Lagrangian depends 
on various free parameters that we now discuss and 
choose. Firstly, the neutrino is not electromagnetically 
charged and from Eq. (38) one has that 



g'Y L cos 6 + gsm.0 = 0. 



(39) 



Secondly, the electron carries the charge e. This means 
that one must recover the following Lagrangian, 



-i*7" (5 M + iQAp) *, 



(40) 



where Q is the electromagnetic charge, Q = e in our 
case. Working out the interaction part of the above La- 
grangian, one obtains 

and, looking at Eq. (38), one deduces that 



(41) 



—Y cost 
2 



9 ■ , 
— sin ( 

2 



= — Y„ cos( 
2 



(42) 



As a consequence, comparing Eq. (39) and (42), one can 
establish that 2F L =5^. It is conventional to choose 
Y h = — 1 and therefore Y R = -2. Then, Eq. (39) 
leads to g' = gtan0, from which one can write cos 9 = 
g/\Jg 2 + g' 2 and sin6> = g'/\/g 2 + g' 2 . Finally, Eq. (42) 
implies that gsinO = e. Everything can be summarized 
by the Gcll-Mann formula 



Q 



Is + l 



(43) 



if the leptonic isospin doublet is assigned a "weak isospin" 
1 = 1/2 such that the neutrino has a third component 
I3 = 1/2 while the electron is such that J3 = — 1/2. Of 
course the right handed particles are such that 1 = 0. 
Returning to Eq. (28) and using the above equations, we 
see that the photons remains mass-less if Y a = 1. This 
also allows us to determine the mass of the gauge bosons 
and one arrives at 



2 

w± 


2 2 

9 v 
2 ' 


m| 





— sin I 
2 



" L w± 
cos 2 6 



(44) 
(45) 



One also easily proves that the coefficient of the cross- 
term Z^A^ vanishes. Finally, looking at Eq. (26), one 
establishes that the Higgs mass is given by 



Xv 2 



(46) 



and, from the expression of the potential (26), one finds 
the expression of its value at the minimum 



V((E) 



11 r 



X 4 ~ 4A ' 



Tf. (47) 



in agreement with what we had before. Using the expres- 
sion of the Higgs mass, one can re-write this expression 
as 



Pvac = =v) = -\m 2 H v 2 



(48) 



We see that in order to calculate the vacuum energy, 
we need to know v and the Higgs mass. Let us first 
discuss how the vacuum expectation value of the Higgs 
can be obtained. For this purpose, we focus on the term 
in the electroweak Lagrangian describing the interaction 
between leptons and charged gauge bosons (the so-called 
charged currents). From the above considerations, see 
Eq. (38), they are given by 



¥ = -^ I /t^W-e L + -^et^W> L . (49) 



V2 



V2 



It is easy to verify by an explicit calculation that it can 
be re-expressed in terms of Dirac spinors^s 



(51) 



1 Let us recall that they are defined by the following expression 

*«=(£)' (50) 
and a similar expression for 



10 



which is similar to the Fermi theory for the weak inter- 
action described by the following Lagrangian 



r - Gf t 

^ Fermi ' 



(52) 



with J\ = J 



_ M 
A 



(i.e. an electronic and a muonic 



part) where the current is given by j( e ' = ^ e j\(l — 

75)^^ an( i a similar definition for j{ . If, for instance 
we consider the decay of muon into an electron, a muonic 
neutrino and an anti electronic neutrino, described by the 
graph 




(53) 



the two theories will lead to the same leading contri- 
bution at small energies if the following identification is 
made 



%m 2 w 



(54) 



where p CTl ~ 10~ 47 GcV 4 is the critical energy density 
today. As will be discussed in the following, this result 
is in fact in contradiction with various observations. 

To conclude this section, let us recall that one can al- 
ways adjust the vacuum energy today to zero by tuning 
the parameter Vq. This is clearly not a very satisfactory 
method and this implies that the vacuum energy density 
was huge prior to the electroweak phase transition. It is 
also worth noticing that this problem is not specific to 
the electroweak transition. Exactly the same discussion 
could be presented for the Quantum Chromo Dynamics 
(QCD) transition [3, 8] (for instance) and would lead to 
the same problem since 



10~ 2 GeV 4 



io 4 W 



(59) 



Moreover, as we are now going to discuss in detail, even 
if we "solve" the classical cosmological constant by tun- 
ing the vacuum energy to zero, see Eq. (8), the problem 
reappears at the quantum level. Indeed, we have seen 
that the zero-point quantum fluctuations also give a con- 
tribution to the cosmological constant, see Eq. (5). The 
goal of the next section is to estimate the corresponding 
energy density. 



Indeed, the appearance of the gauge boson mass comes 
from the propagators of W^. At small momenta com- 
pared to the mass, only mw remains and the above 
graph is in fact equivalent to the graph corresponding 
to the Fermi theory. Then, using the expression of the 



VF-bosons mass established before, m 2 ^ 
arrives at 

2 V2 



g 2 v 2 /2, one 



(55) 



Notice that there is a factor y/2 difference with respect 
to the textbook [20] because we have defined the Higgs 
as v + Hj \/2 instead of (v + H) j v2. The Fermi constant 
is given by G F ^ 1.16 x 10~ 5 (GeV) -2 . Therefore, 



175 GeV. 



(56) 



or \/2v ~ 246 GeV. We see that we have been able to 
calculate the Higgs vacuum expectation value v. As a 
consequence, Eq. (48) reads 



\/2™| 
16 G 2 , ' 



(57) 



The mass of the Higgs boson is not known although, 
at the time of writing, there are reasons to believe that 
rn H < 129 GeV at 95% CL [46] (there are even reasons 
to believe that m H ~ 125 GeV but this is not yet es- 
tablished at a sufficient statistical level at the time of 
writing). Therefore, we find that 



-1.2 x 10 8 GeV 4 



-10 55 p c 



(58) 



IV. THE QUANTUM-MECHANICAL 
COSMOLOGICAL CONSTANT PROBLEM 

We have seen that, even if we solve the classical cos- 
mological constant problem and find a convincing reason 
to put the minimum of the potential to zero, p vac is a 
quantity which also receives contributions from the zero- 
point fluctuations of all the quantum fields present in the 
Universe. We now discuss this "quantum-mechanical" 
cosmological constant problem. In principle, the corre- 
sponding vacuum energy density can be estimated from 
first principles. The goal of this section (and of Sec. VI) 
is to carry out this task. We will first discuss how 
the Lorentz invariancc affects the calculation of the vac- 
uum energy, see the next subsection IV A. Then, in sub- 
Sec. IV B, we will show how this question can be formu- 
lated in terms of Feynman diagrams. In sub- Sec. IV C, 
we will also evaluate the vacuum energy not only for a 
free theory as usually done but also in the more realis- 
tic case where interactions are present. Finally, in order 
to understand better the origin of the vacuum energy 
density, in sub- Sees. VA and VB, we will address the 
same question but in the context of ordinary quantum 
mechanics. This will allow us in sub- Sec. VC and VD 
to discuss in more detail the properties of the so-called 
bubble diagrams from which the cosmological constant 
originates. As mentioned above, all the fields present in 
the universe participate in the vacuum energy. There- 
fore, we need to do the calculation for scalar, fermion 
and vector fields. Here, in a first step, we carry out the 
calculation for a real scalar field. The other types will be 
treated later on in Sec. VII. For simplicity, we will also 
consider that the metric is flat and that the fields live in 
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Minkowski space-time. A sloppy justification is that the 
main contribution to the vacuum energy density comes 
from modes with very high momenta, corresponding to 
scales at which the curvature of space-time is negligible. 
In fact, we will come back to this very important ques- 
tion in more detail in what follows and will try to discuss 
this point in a more rigorous way, see Sec. VIII. 



A. The Zero-Point Energy Density 

Let us consider a simple real free scalar field with the 
potential = m 2 <I> 2 /2 where to is the mass of the 

scalar particle. In flat space-time, the equation of motion 
is nothing but the Klein-Gordon equation, namely 



$ + S* 3 didj$ - m 2 $ = 0. 



(60) 



Since the scalar field is free, the equation of motion is 
linear and, as a consequence, one can Fourier expand 
x) as 



$(t,aO = 



1 



d 3 fc 



c k e 



—iu)t-\-ik-x 



+c t e fa,t-ifc.» 



with 



i(k) — sj k 2 + to 2 . 



(61) 



(62) 



The four-dimensional momentum has been written as 
k^ = (k°,k) and the integration is performed over its 
spatial part only. We have used the notation \k\ = k. Its 
time component is given by k° = u. As usual, the field 
has been quantized by considering c& and c k (the so- 
called annihilation and creation operators) as quantum 
operators satisfying the following commutation relation 



c fc ,4, 



6^ (k - k') . 



(63) 



This equation is in fact equivalent to the standard com- 
mutation rule between the field operator and its conju- 
gate momentum. 

We are now in a position where one can compute 
(0|T M „|0). From Eq. (7), we see that we must first eval- 
uate the mean values of various quantities depending on 
the field operator and its derivatives. Straightforward 
calculations lead to 



(0|$'|0) 
(0\S ii d i $d j $\0) 
(0|$ 2 |0) 



1 



(2tt) 3 
1 

(2^)3 

1 



d 3 fc 

2uj(k) 
d 3 fc 
2uj(k) 
d 3 fc 

(2tt) 3 J 2u(k) ' 



u, 2 (fc), 



(64) 
(65) 
(66) 



These equations can then be used to determine the en- 
ergy density and pressure of the vacuum. From the above 
expression (7) of the stress-energy tensor, one has 



H 



_$2 + + I m ^$ 



(67) 



where % denotes the Hamiltonian density. This leads 
to the following expression for the energy density p = 
u^u^T^v measured by a fundamental observer character- 
ized by its four-velocity vector u 1 * = (1, 0) 



(p) = <0|u o u o r oo |0) 



1 1 



(27r) d 2 



d 3 fcw(A;). 



(68) 



In the same manner, one can easily estimate the pressure. 
One obtains 



(P) = 



I$2 _ ItfjQ QQ Q _ I m 2$2 

2 6 2 



0), (69) 



where J_v= g^ v + u^u v is a projector. Upon using the 
expressions (64), (65) and (66), one arrives at the follow- 
ing expression 



(P) 



1 1 /"jS, fc2 



(ky 



(70) 



(2nf 6 
Notice that we also have 

(0|T W |0) - --^l J d 3 kki = 0, (71) 

and, hence, we do not need to consider to off-diagonal 
elements. This completes our calculation of the vacuum 
stress-energy tensor. 

From the above considerations, using Eq. (11), we de- 
duce that 



A eff =A B + -^-3 y"d 3 fci w (fc). 



(72) 



The obvious problem with the previous calculation is 
that the integrals expressing the energy density (p) and 
the pressure (p) blow up in the ultra-violet regime (and, 
therefore, strictly speaking, the effective cosmological 
constant is in fact infinite). This is of course a well- 
known problem in quantum field theory. Usually, this 
type of divergences is ignored (or sometimes removed by 
mean of the normal ordering product) because, in ab- 
sence of gravity, only differences of energy can be de- 
tected [16, 17]. However, in the present context, this is 
clearly no longer the case since the absolute value of the 
vacuum energy weighs and, therefore, can be measured. 
This is why regularizing these infinities becomes a crucial 
issue that we now discuss in some details. 

The common method found in the literature is to in- 
troduce a cut-off at k = M, where the physical interpre- 
tation of M is that this is the scale at which the effective 
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theory used before breaks down [2, 4, 5, 7, 9]. The quan- 
tity M is not known: it could be the Planck scale, the 
string scale or even the super-symmetric breaking scale. 
As we will see later on, this is in fact not important be- 
cause the cosmological constant problem is present re- 
gardless of the precise value of the cut-off. In this ap- 
proach, the energy density becomes 

rM 



(p) 



1 

M 4 
16tt 2 



dkk 2 \/ k 2 + rn 2 



(73) 



77l z 
IT 1 



1 m 2 
1+ 2M 2 



1 m 4 / M 

2 M 4 n \ m 




M 4 
16-k 2 



1 



m 
~M 2 



(74) 



(75) 



where, in the last expression, we have expanded the exact 
expression in terms of the small parameter m/M . We see 
that the divergence is quartic in the cut-off scale. Let us 
now perform the same calculation for the pressure. Using 
Eq. (70), one obtains 

rM k A 
dk- 



Vfc 2 



(76) 



1 M 4 



3167T 2 

3 m 4 
+ 2~M 4 

1 M 4 



M 2 



3 to 2 
~ 2M 2 
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m 2 \ 
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" M 2 J 



3 16tt 2 



1 



m 
M 2 



(77) 



(78) 



It is clear from the previous expressions that (p)/(p) 7^ 
— 1 which indicates that the stress energy tensor is not 
of the form cx —pg^ v . In the limit m — > 0, as can be 
easily shown from Eqs. (75) and (78), the equation of 
state is in fact (p)/(p) = 1/3. This would mean that the 
zero point fluctuations do not behave like a cosmological 
constant but rather like radiation. On the other hand, 
this result seems to be strange since we know from fun- 
damental considerations that the stress-energy tensor of 
a system in its ground state must be proportional to the 
metric tensor. Clearly, something has gone wrong is our 
calculation. 

Interestingly enough, let us also note in passing that if 
one only focuses on the logarithmic terms, namely 



32tt 2 



M M , 
H- + -\/l + TT?5 h ( 79 ) 



(80) 




then the equation of state of the vacuum is correctly re- 
produced (even if the energy density can be negative) . In 
fact, as we explain below, this is not a coincidence [47]. 



It is easy to check that the previous considerations are 
not an artifact of the fact that we have introduced a 
sharp cut-off. Indeed, if rather than a sharp cut-off in 
the upper limit of the integrals, one introduces a smooth 
exponential cut-off, one obtains (for simplicity, one only 
considers the case m = 0), 



(P) = TZo 



1 

4^ 2 
1 

i^ 2 



dfc k A e- 



1 r(4) 

4tt 2 a 4 



6M 4 



(81) 
(82) 



where we have written a = 1/M and where T(z) is the 
Euler's integral of the second kind [48, 49]. One re- 
covers our familiar quartic divergence. Doing the same 
manipulation for the pressure, one obtains again that 

(P)/(P) = 1/3. 

In fact, what has gone wrong is that we have used 
schemes of regularization that do not respect the sym- 
metry of our underlying theory [47, 50, 51]. Indeed, the 
Lorentz invariance is broken when we impose a cut-off on 
the spatial momentum only. As exemplified by Eq. (82), 
it is easy to show that this is a feature of any regulariza- 
tion scheme breaking Lorentz invariance and not a spe- 
cific property of the particular cut-off scheme used above 
(and usually used in the literature). In fact, this prop- 
erty is well-known in the context of gauge theories. If 
one regulates divergent graphs with a method that is not 
gauge invariant, then one obtains incorrect results [52]. 
Here, the situation is very similar. 

Therefore, in order to meaningfully evaluate the zero- 
point energy density, one must use a regularization 
scheme that does not break Lorentz invariance [47, 50, 
51]. An obvious choice is dimensional regularization 
and we now explore this route. In order to use this 
method, we first reformulate our problem in a d di- 
mensional Minkowski space-time. Obviously, the Klein- 
Gordon equation (60) remains unchanged except that the 
indices of the Kroncckcr symbol now run from to d — 1 . 
The Fourier expansion of the field becomes [compare with 
Eq. (61)] 



1 



(27T) 



(<f-l)/2 



1 „t iut—ik-. 



(83) 



and, as a consequence, the energy density can be ex- 
pressed as 



(P) - ^4 / d-M*) (84) 
dkk d ~ 2 d d - 2 nuj(k), (85) 



A-d 



(2n) (d - 1) 2 Jo 

where we have introduced a scale p, in order for the equa- 
tion to be dimensionally correct. This equation should 
be compared with Eq. (68). Then, using the fact that 
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the angular integrals can be easily performed, namely 
Jd d - 2 n = 27r( d - 1 )/ 2 /T(d/2 - 1/2), where, again, T(z) 
is the Euler's integral of the second kind [48, 49], one 
obtains 



(p) 



T(-d/2) fm 



2(47r) (<i-l)/2r(-l/2) 



(86) 



Performing the same calculation for the pressure, one 
arrives at a similar expression 



0') 



4-d 



1 



(2 7 r) (d - 1) 2(d-l) 



T(-d/2) fm 



4(47r) «-D/a T(l/2) Vm, 
But r(— 1/2) = — 2r(l/2) and, as a consequence, 
(p) = "(P>- 



(87) 
(88) 

(89) 



This time, as expected, we have obtained the equation 
of state of the vacuum. The fact that the regularization 
method used satisfies the Lorentz symmetry has led us 
to the correct result. 

The consistency of the above result can be checked in 
a different manner [47]. Using Eqs. (84) and (87) for 
the energy density and the pressure, one sees that they 
satisfy the following differential equation 



(p) — 2m 



2 d(p) 
dm 2 



(d-l)(p). 



Using the fact that (p) = — (p) one obtains 



dm 2 



(90) 



(91) 



which can be integrated to give (p) oc m d in full agree- 
ment with Eq. (86). Therefore, we conclude that the 
above considerations are all consistent. 

Let us now analyze in more details the structure of the 
divergences. For this purpose one writes d = 4 — e, where 
e is supposed to be a small parameter, and we expand 
(p) in terms of e. The Euler function is such that 



1 



1 



1 



-2 + e/2-l + e/2e/2 



-r i 



(92) 



with r(l + e/2) ~ r(l) + eT(l)/2, r'(i) = #(i)r(i) 
where ^> is the Polygamma function [48, 49] and ^(1) = 
-7 ~ 0.57772 the Euler-Mascheroni constant [48, 49]. If, 
in addition, we use the two following expansions 



(47T) 



-3/2+e/2 



1 



4-£ 



(4tt) 3 / 2 
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(93) 
(94) 



then one arrives at the following expression for the energy 
density of the vacuum 



(p>~- 



64tt 2 



2 3 

— I 7 — In , 

e 2 ' V 47rp 2 



(95) 



Using a MS rcnormalization scheme convention [MS 
means "modified minimal subtraction" and simply con- 
sists in subtracting the pole ~ 1/e together with the 
accompanying terms 7 and ln(47r), see chapter 11 of 
Ref. [16]; in fact, here, we slightly modify the convention 
such that the term 3/2 is subtracted as well], the regu- 
larized, finite, energy density of the vacuum now reads 



(P) 



64tH 



■In 



(96) 



This result is very different from the result obtained by 
imposing a sharp cut-off. In fact, one recovers the struc- 
tures of the logarithm terms found in Eqs. (79) and (80) 
since it is well-known that the dimensional regulariza- 
tion scheme removes the power law terms. As discussed 
around Eqs. (79) and (80), this is the reason why those 
terms lead to the correct vacuum equation of state. We 
see that the result is not proportional to the cut-off scale 
to the power four, as usually claimed on dimensional 
grounds, but to the mass of the particle to the power 
four, which is also dimcnsionally correct. Physically, for 
instance, this means that the photon does not contribute 
to the vacuum energy density contrary to what a sharp 
cut-off calculation would predict. Therefore, the order of 
magnitude of the final result can a priori be very differ- 
ent. We also notice that the sign of (p) can change ac- 
cording to m > p or m < p. Of course, it remains to be 
seen whether this could solve the cosmological constant 
problem and we will argue below that, unfortunately, the 
answer to this question is no. It should also be kept in 
mind that the above calculation is valid for free fields. It 
is clear that a more realistic version should take into ac- 
count the interactions and we will also consider this issue 
in the following. However, before addressing these vari- 
ous points and in order to gain further physical insights, 
we would like to investigate again the above issues, this 
time from the Feynman diagrams point of view. We turn 
to this question in the next sub-section. 



B. The Vacuum and the Feynman Diagrams 

We have seen in the last sub-section that the vac- 
uum contribution to the cosmological constant is given 
in terms of divergent integrals. A priori, there is no rea- 
son to be surprised since divergent integrals are routinely 
found in quantum field theory. We have the method of 
rcnormalization at our disposal to get rid of these infini- 
ties and the so accurate predictions of particle physics 
are here to remind us that this is a very convincing way 
of taming these divergences. However, we will see that 
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the divergent integrals that we encounter when we calcu- 
late the vacuum energy are in fact of a different nature 
than the divergent integrals met when we calculate, say, 
a cross section. In order to discuss this point in more 
detail and to understand better the structure of these di- 
vergences, we now briefly return to some textbook con- 
siderations [1G-22]. 

In quantum field theory, the basic object is the generat- 
ing functional defined by the following path integral (we 
still restrict our considerations to a simple scalar field) 

Z [J] = N [ X>$ cxp { i I d 4 x [C + J(x)®(x)] 1 , (97) 



where Af is a number chosen such that Z\ J = 0] = 1. The 
quantity C is the Lagrangian density while J represents 
a source coupled to the scalar field. From the generat- 
ing functional, one can compute the various correlation 
functions 

G<") (an •••£„) = (n|T[$(ari)..-$(a; n )]|fi} (98) 

1\" 8 n Z[J] 

i 



8J(xi) ■ ■ ■ 8J(x n ) 



(99) 



where |fi) denotes the vacuum state of the theory. We 
distinguish it from |0) which denotes the vacuum state 
of the theory without interaction. In this last case, the 
generating functional can be computed exactly. It reads 
(notice that we have Zq [0] = 1 and that the "0" subscript 
indicates that we refer to the free theory) 



Z [J] = exp 



-\ J d i xd i yJ{x)D F (x-y)J{y) 



(100) 



where D F is the Feynman propagator defined by the fol- 
lowing expression 



D F (x 1 -x 2 ) = (n\T[$(x 1 )*(x 2 )]\n), 



(101) 



T denoting the time ordered product. As is well-known, 
the propagator can be conveniently expressed in Fourier 
space as 



£> F (Xl - X 2 ) 



d 4 /fc 



(2tt) 4 J k 2 + 



(102) 



In presence of interactions, one can write C = £o + Ant 
where £; n t represents the self-interacting part of the La- 
grangian density. In this case the generating functional 
can be written as 



Z[J] = AAcxpji J d A xC int 



1 8 



i SJ(x) 



Z [J], (103) 



where, in order to have the correct normalization, the 
factor Af must be given by 



Af 1 = exp 



d 4 x £ int 



1 



i 8J(x) 



Z [J] 



(104) 



As explained in the standard quantum field theory text- 
books, the above expressions (103) and (104) are partic- 
ularly well suited for a perturbativc expansion. For in- 
stance, if we want to calculate the generating functional 
for the case of a quartic theory = A<I> 4 /4!, where 
A is the coupling constant, we are led to the following 
expansion in A 



Z[J] = AT 



1 



A 



d 4 x 



1 8 A Z [J] 
i 4 SJ(x) 4 



(105) 



As usual, the various terms that appear in this expansion 
can be graphically represented by Feynman diagrams. 
The Feynman rules of the theory in real space are given 
by 



D(xi - x 2 ) 



iX / <Tx 



x 2 , 




i d xJ(x) 



(106) 



(107) 



(108) 



and a loop for a propagator evaluated at the same point 
in space-time. This completes our quick reminder of the 
basics of field theory that we need in the following. 

Now, we compute the energy density of a free scalar 
field (the interacting case is treated in the next subsec- 
tion) using the tools that we have just presented. In fact, 
it is interesting to calculate the trace of the stress-energy 
tensor. Upon using Eq. (7), one obtains 



(T) = WT, l/ ) = -(p) + (d-l)(p). 



(109) 



Then, since the energy density and the pressure are given 
by Eqs. (84) and (87), the above relation can be re- 
cxprcssed as 



(T) = -n 



4-d 



(2ir) d - 1 2uj(k) 



(110) 



But, from Eqs. (101) and (102), and now specifying to 
the case d = 4, one has 



DM 



d 4 fc 



dfc°d 3 fc 



(2tt) 4 J k 2 + m 2 (2tt) 4 J -(fc°) 2 +w 2 ' 

(111) 

where we have distinguished the integration over the time 
and space components of the momentum. Using the iden- 
tity / dk°/[-(k ) 2 +uj 2 } = -tt/O), see also Eq. (155), 
this leads to 



D p (0) 



d 3 fe 



(2tt) 3 2w' 



or, in other words, 



j=o 



(T) = m 2 D F (0). 



(112) 



(113) 
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But, if we now use the fact that (p) = —(p), Eq. (109) 
can also be written as (T) = —d(p). As a consequence, 
Eq. (113) implies that 



(p) = -—D F (0). 



(114) 



This equation can also be obtained by a direct calcula- 
tion. Indeed one has 



iji 



i-d 



dk°d d - 1 k 



(2ir) d J -(k°) 2 +uj 2 



r(-d/2) (m 



(115) 



(116) 



™ 2 2(4^) (d - 1)/2 E(-l/2) 

= -4^ ( n7 ) 

where, in the second step of the calculation, we have used 
dimensional regularization. Now, using Eq. (11) and the 
Feynman rule stipulating that a propagator evaluated at 



the same point is represented by a loop, one can re- write 
the cosmological constant as 



A ff = A B 



40 



(118) 



In other words, the cosmological constant is given by a 
"bubble diagram", i.e. a diagram that has no external 
leg. 

We have just established that the vacuum energy den- 
sity, or equivalently the cosmological constant, can be 
expressed in term of a bubble diagram. This type of di- 
agram plays a special role in quantum field theory and 
we will discuss this point in more detail in the following 
next sub-sections. In order to develop our intuition, let 
us first calculate explicitly the generating function in the 
case of a self-interacting scalar field. The calculation is 
straightforward and upon using Eq. (100), this leads to 



i- I d 4 x 



1 S 4 Z [J] 
i 4 6J{x) 4 



-ijA 3 J d 4 xD 2 (0) -6 J d i xD(0) 



d 4 yD(x-y)J(y) 



d A x 



d 4 yD(x - y)J(y) 



>Z Q [J] + 



( 



1 

4! 



3 0O+6 



V 




Z [J}. 



(119) 



(120) 



One notices the appearance of a "bubble diagram" (the first one), similar (i.e. with no external leg) but not identical 
to the one encountered in the calculation of the cosmological constant. We also have a diagram with a loop (the 
second one) but with external legs and a diagram with no loop (the last one). All the diagram with loops (that is to 
say the two first ones) are divergent in quantum field theory. Nevertheless, the bubble diagram and the loop diagram 
play in fact a very different role as we are going to see. Let us now calculate the normalization Af, see Eq. (104). It 
uses the same expression as before but with J = which means that all the diagram containing a source disappear. 
As a consequence, we have 



J\f~ 



1 



3 
4! 



oo 



(121) 



We can now gather everything and compute the generating functional. Upon using the expression (105) and the 
results (120) and (121), one arrives at 




(122) 



I 

We see that the bubble diagram has canceled in the fi- nal expression. Since all the predictions (cross-sections 
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etc ...) follow from the calculation of the generating 
functional, we see that the bubble diagrams never con- 
tribute, in absence of gravity, to any observables. On the 
contrary, when the gravitational field is turned on, the 
bubble diagrams become important and can affect the 
predictions of the theory. In summary, the cosmologi- 
cal constant problem can be viewed as the appearance of 
a new type of divergent graphs. As demonstrated with 
Eq. (96), these graphs can be regularized but, as will be 
discussed below, the regularized value seems to be incom- 
patible with some astrophysical observations. 

C. The Vacuum in Presence of Interactions 

So far we have calculated the vacuum energy density 
for a free field. It is natural to ask what happens when 
one considers interaction [50] since we know that a re- 
alistic model necessarily contains such terms. For this 
purpose, let us therefore consider again a simple self- 



interacting scalar field with & n t = A<E> 4 /4!. As a conse- 
quence, this adds the following contribution to the energy 
density of the vacuum 

A " = a^'Hir^^O) (123) 

- OO m 



where, on the first line in the above expression, we have 
used the fact that the vacuum wave-functional is a Gaus- 
sian. As expected, this new term is also given by a bubble 
diagram. However, one should not forget that, when in- 
teraction is present, one must also renormalize the mass. 
As is well-known, at one loop, the renormalization of the 
mass comes from the tadpole diagram. More precisely, 
one has to calculate the correction to the propagator. It 
is given by the following expression 



(n|T[*(a:i)#(x2)]|n) = xim- 
or, in terms of explicit expressions [see Eq. (102)] 



1 



•<2 



<n|T[*(a;i)$(a; 2 )]|n> = D F ( Xl - x 2 ) 



~2 



d 4 xD F ( Xl - x)D F (0)D F (x-x 2 ) 





r d 4 P 


(2tt) 4 j 


P 2 + TO 2 




r d 4 P 


(2tt) 4 j 






f 


(2tt) 4 j 


p 2 + m 2 



2 pW (2tt) 4 



1 



d 4 p 



&iPii'{ x i ^2) 



(p 2 + m 2 ) 2 



a tPp'(if-xf) 



d 4 p 



1 



A 



2 - r 



where the last equation is equivalent to the previous one 
at leading order in the coupling constant. We see that 
this implies the following redefinition (or renormaliza- 
tion) of the mass 



I 

the field, can be re-expressed as 
™2 A 



(P) 



777 



F (o) 



= m 2 -^ F (0) 



(130) 



(125) 

(126) 
(127) 
(128) 
(129) 



(131) 



As a consequence, the vacuum energy which is now the 
sum of the free contribution (114) and of the contribu- 
tion (123), which originates from the self- interaction of 



where, in the second line, we have used Eq. (130). The 
extra contribution coming from the interacting term is 
exactly canceled by the renormalization of the mass [50, 
53]. In other words, one sees that the presence of the 
interaction has not modified the expression of the vacuum 
energy density (114), at least at one loop, provided the 
final result is expressed in terms of the renormalizcd mass 
rather than the bare one. 
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THE BUBBLE DIAGRAMS 



The previous sections have shown the central role 
played by the bubble diagrams in the vacuum energy 
problem. These diagrams do not affect the physical pre- 
dictions when the gravitational field is turned off and 
this is why, in standard quantum field theory, there are 
usually ignored. In this section, we wish to discuss in 
more detail this type of diagram. For this purpose, we 
will analyze them in ordinary quantum mechanics where 
the technical aspects are less complicated and where the 
interpretation is easier [54-56]. 



A. Quantum Mechanics as a Field Theory 

Let us start with the simplest quantum-mechanical 
system, namely the one-dimensional harmonic oscillator. 
Let us denote by x(t) the position of the particle evolv- 
ing into a parabolic potential. Then, the corresponding 
Lagrangian is given by 



L(x, x) 



1 



1 



-mx 



9 9 

muj x 



(132) 



where m is the mass of the particle and to is the effective 
frequency. Then, let us define the "scalar field" 4>(i) by 



*(t) = y/mx(t) , 



(133) 



which implies that $ has dimension —1/2. As a conse- 
quence, the Lagrangian now reads 



L = $- 

2 2 



(134) 



This Lagrangian is similar to the one one encounters in 
field theory for a free scalar field. This means that all the 
techniques of quantum field theory can in fact be used in 
ordinary quantum mechanics. Of course, this does not 
lead to genuine new results since both approaches are 
equivalent but this can shed a new light of some physi- 
cal results. This is the strategy that we use here applied 
to the vacuum energy problem. We now quickly remind 
how the formalism of quantum field theory can be imple- 
mented in quantum mechanics. We will also compare this 
approach to the standard approach of quantum mechan- 
ics. The first step is to define the conjugate momentum. 
It is given by 



r)T 

II(i) = -=: = *(*). 
and this leads to the following Hamiltonian 



H = 



U 2 



-<f> 



(135) 



(136) 



The next step is to introduce the creation and annihila- 
tion a and c operators. They are defined by the exprcs- 



n 



'2uj 



(137) 
(138) 



As usual one can check that [$, II] = i implies [c, = 1. 
As it is standard in field theory, we now work in the 
Heisenberg picture where the operators are time de- 
pendent. A straightforward calculation shows that the 
Hamiltonian can be expressed as 



H = — (cc 1 + c^c) 



(139) 



This allows us to calculate the time evolution of the op- 
erator c. Indeed, the Heisenberg equation reads 



*di = ^ C,H ^ = wc ' 



from which wc deduce that 



c(t) = c(0> 



c e 



(140) 



(141) 



As a consequence, the field operator can be written in 
terms of the creation and annihilation operators in a way 
which is very similar to Eq. (61), namely 



*(*) = 



1 



c e 



(142) 



Of course we do not have an integral over the Fourier 
modes because, in quantum mechanics, we deal with a 
single harmonic oscillator as opposed to an infinite col- 
lection of harmonic oscillators in field theory. 

At this level, the theory is known exactly. It is clear 
that the most interesting part is in fact the treatment 
of the interactions which, in field theory, is based on the 
calculation of the Feynman diagrams. Since this calcula- 
tion is itself based on the calculation of the propagator, 
one must now evaluate it for our simple system. It is 
defined by an equation similar to Eq. (101), namely 



D{hM) = (0|T [#(ii)*(t 3 )] |0) , 



(143) 



where we recall that the time ordered product is defined 
by 



r[$(tx)$(i 2 )] 



$(ti)$(t 2 ), h>t 2 
$(t 2 )*(ti), t 2 >t 1 . 



Then, the calculation proceeds straightforwardly and one 
obtains for t± > t 2 



D(h,t 2 ) = (0 



1 



c e 



c e 



-iu(ti —t 2 ) 



—iujt2 i „t^^i; 



2co 



) (144) 
(145) 
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FIG. 3: Contours in the complex time plane used to calculate the Feynman propagators of the one-dimensional harmonic 
oscillator. As usual, there are two poles on the real axis that can be avoided by going along the two small red circles. The final 
contour is closed by a large blue circle of radius R. 



and, in a similar way if i 2 > ii 



c e iw > 2 + c, 



1 



1 



_ iw(ti-t 2 ) 

' 2w 



o e J 

0^ (146) 
(147) 



The two previous equations can be summarized into a 
single one 



D(t 1 ,t 2 ) = 7 le-^I* 2 -* 1 I 
2w 



(148) 



We see that the system is so simple that we obtain an ex- 
plicit expression of the propagator in real space (in fact, 
in field theory, there is also an expression of the prop- 
agator in real space but it is not so often used because 
of its complexity - it is given in terms of ordinary and 
modified Bessel functions -). Nevertheless, in order to 
compare with the standard approach and Eq. (102), it is 
interesting to obtain the expression of the propagator in 
Fourier space. The Fourier transform of a function f(t) 
is defined by 



/(*) 



1 



27T J- 



f(E) e M dE 



(149) 



The calculation is easy and well-known. Here, for con- 
venience, we quickly remind how it proceeds. For that 



purpose, let us consider the function 



(150) 



It has two poles at z = ±cj on the real axis, see Fig. 3. 
We assume that t > and consider the following contour: 
along the real axis with two small circles around the 
two poles (the circle around z = — uj goes into the lower 
part of the complex plane while the circle around z = +lu 
goes into the upper part); and we complete by a big circle 
going into the lower part of the complex plane. In 
this case, only the pole z = +lj is inside the contour 
r = lUT±U T^, see Fig. 3. We have 



lim \zf(z)\ 



zc 



lim 

|z|->oo Z" 



(151) 



since Q{z) < in the lower part of the complex plan. 
According to the first Jordan lemma, this means that 
the integral on the big circle vanishes. In the same way, 
we have lim z _>.o \zf(z) \ = and the second Jordan lemma 
tells that the integral on Tf is also zero. Then, using the 
fact that f(z) = P(z)/Q(z) and that, in this case, the 
residue can be written as Res(/, z) = P(z)/Q'(z), one 
obtains 



Res(/, z 



-iujt 



2oj 



(152) 
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As a consequence 

f(z)dz 



lim 

R— ¥+00 



I 

2^ 



1 

i 

2^ 



2w ' 



f(z)dz 

O 

-2z7rRcs(/, z = 

(153) 



where the minus sign in the second line comes from the 
fact that the contour is clockwise (i.e. opposite to the 
direct direction). 

Now if t < 0, then the big circle is taken in the upper 
plane. As a consequence, the pole inside the contour is 
now the one located at z = —u>. Therefore, 



1 

2^ 



f(z)dz 



1 

27 



f(z)dz 



— x +2i7rRes(/, z = -co) 



2uj 



(154) 



The previous considerations allow us to conclude that 

e — iw|t| 



I 

2^ 



+ 00 Q-iEt 



E 2 



■dE = 



2uj 



D(t), 



(155) 



that is to say we have established the Fourier transform 
of the propagator. It should be obvious that the above 
equation is the counterpart of Eq. (102). As in field the- 
ory, we will show that it can be used to develop a per- 
turbative method which allows us to take into account 
interactions. 



B. Perturbations Theory in Quantum Mechanics 

Before discussing our main subject, we need to explain 
how the Feynman diagrams can be used in quantum me- 
chanics [55, 56]. Usually, in quantum mechanics, the 
perturbation theory is not based on the Feynman dia- 
grams [57]. This gives the impression that the techniques 
of field theory are very different. In fact, as we now 
demonstrate, these techniques can also be utilized in or- 
dinary quantum mechanics [54] and, obviously (as could 
have been guessed from the very beginning), they lead to 
the same physical predictions. In the following, we will 
discuss this question and compare the two approaches 
in some details. We believe that this can improve our 
physical understanding of the vacuum energy problem. 

In order to achieve the above mentioned task, we need 
to recall how interactions are treated in quantum me- 
chanics. For this purpose, let us consider a system obey- 
ing the following Schrodinger equation 



1 dt 



H(p 3 ,x s )\V(t)) s , 



(156) 



where H(p s , x s ) is the Hamiltonian operator. The above 
equation is written in the Schrodinger representation, 
hence the subscript "S" . Then, one can define the evo- 
lution operator U(t,t{) which relates the state vector at 
the time t to the state vector at some initial time t\ < t 
by the following expression 



|¥(t)> s = ET(t,ti)l*(*i)> E 
This operator obeys 

.dU(t,ti) 



dt 



HU(t,U), 



(157) 



(158) 



and, provided the Hamiltonian is not explicitly time de- 
pendent, it can also be written as 



U(t,ti) = e 



_ lH ( t - ti ) 



(159) 



Of course, in this representation, the operator p s and 
x s are time independent since only the state can change 
with time. 

One the contrary, in the Hcisenberg representation, 
the state vector does not evolve and can be expressed in 
terms of the Schrodinger state vector evaluated at some 
reference time t = to, 



I*)h = l*(*0)> E 



(160) 



Notice that the reference time is arbitrary and is not 
necessarily t[. In this picture, the operators are time- 
dependent and, for instance, x H (t) is given by 



x B (t) = Ui(t,t )x a U(t,t ), 



(161) 



such that x H (*o) = x s- I n this picture, the time evolu- 
tion of the operators is computed from the Heisenberg 
equation which reads 



dt 



ouHt,t ) f dU(t,tp) 
x s U(t, to) + iu'(t, t )x s 



dt 



dt 



[x s ,H], 



(162) 
(163) 
(164) 



where in the second line we have used that —idW/dt = 

Let us now consider the case where interactions are 
explicitly present. We assume that the corresponding 
Hamiltonian can be split into two parts, one correspond- 
ing to the free theory and one corresponding to the in- 
teractions 



H = Hn 



int ■ 



(165) 



We also assume that H lnt is a small perturbation. The 
free evolution, due to the free Hamiltonian Hq, is de- 
scribed by the evolution operator Uq. This operator is of 
course different from U which describes the full evolution 
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with the interactions taken into account. Clearly, Uq can 
be expressed as 



U (t,ti) 



-iH (t-ti) 



(166) 



Then the idea is to factor out the free evolution (which is 
supposed to be known) and to focus on the time evolution 
due to the perturbation on top of the free evolution. For 
this purpose, we define the interaction representation by 
(again, the representations are supposed to coincide when 
t = to) 



Xl {t) = U^(t,to)x s Uo(t,t Q ) 



e iH (t-t ) x e -iH (t-t )_ 



(167) 
(168) 



Clearly x 1 (t) is different from x H (t) because its time evo- 
lution is "generated" only by the free Hamiltonian and 
not by the total one. Let us now calculate the equation 
of motion of x 1 (t). One has 



— = iH qXi 
= i [-Hq,£i] 



e iH (t-t ) x 



-iH e- tHo{t - to ll69) 
(170) 



Therefore, the time evolution of x I (t) is controlled by 
the free Hamiltonian. This evolution is supposed to be 
known. In the same manner, we define the state vector 
in the interaction picture by 



|*(t)> t = ^(4, *o)|*(t)) E 



(171) 



As a consequence, in the interaction picture, both the 
state vectors and the operators are time dependent. The 
equation obeyed by the state vector is given by 



' dt 



--Hoi*}, +e iHo(t - to 'i^ (172) 

= -H c lffo(t - to) |*) s 

+e iH o(t-to)^ Ho + H . nt ^ s (173) 

= e iHo(t - to) H int \^) s (174) 
= c lH ^ t - t ^H int e^ Ho{t ' ta '>\^) I (175) 
= Hi|*>i, (176) 

where H I = UQ(t,to)Hi nt Uo(t,to) in agreement with the 
definition (167). The time evolution of is therefore 
generated only by the interacting Hamiltonian. Then, 
it is natural to define the evolution operator in the in- 
teraction picture. This operator is associated with the 
Schrodinger equation (176). Its definition reads 



(177) 



Of course U 1 is different from Uq (and from U). A 
straightforward calculation shows that it can also be ex- 
pressed as 



U r {t,ti) = Texp 



■J AtH^t) 



(178) 



Another expression for U l is based on the following con- 
siderations. From Eqs. (177) and (171), we have 

!*(*)>, = fi)|*(ti)>i (179) 
= (7 I (t,i i )(7 t (t i ,io)|*(ti)) s (180) 
= Ut(t,t )\y(t)) s . (181) 

Therefore, one obtains 

|*(t)) B = C/ (t,to)C/ I (i,ti)^o t ( i i^o)|*(ii)) s (182) 
= tf(f,ti)|*(*i)>s> (183) 
from which we deduce that 

U(t,ti) = Uo&tap^tjuUtuta), (184) 

or 



JH (t-t ) p -iH(t-ti) -iH {ti-to) 



(185) 
(186) 



This expression is very natural since it expresses nothing 
but the fact that the full evolution of the system is just 
the combination of the free evolution and of the evolution 
due to the interacting term. Let us also remark that one 
can write 

Xl {t) = U^(t,t )x s Uo(t,t ) (187) 
= U^(t,to)U(t,to)x u U\t,to)U (t,t ) (188) 



JH (t-t ) -iH(t-t ) , 



" u, x a (t)e"" 



[/.(Moj^W^fMo). 



iH{t-t ) -iH {t-t ) 



(189) 
(190) 



Again, this result is very intuitive and expresses the fact 
that, in the interacting picture, the free evolution of the 
operators is factorized out. 

After these preliminaries, one can now derive the Cell- 
Man Law equation which is at the heart of the pertur- 
bative treatment presented in what follows. Let |n : ) be 
the energy eigenstates of the full Hamiltonian, i.e. the free 
Hamiltonian plus the Hamiltonian describing the interac- 
tions. We denote the true vacuum state by |0,) = in 
agreement with the conventions adopted in the previous 
subsections. Then, one can always expand the vacuum 
state of the free theory according to 



|o) = Ek>^i°>< 



n=0 



|Q)(Q|0) + E|n r )(n I |0). 



(191) 
(192) 



As a consequence, applying the operator c lH ( t + t o) to 
the above relation, one obtains 

e - t H(t+t )^ = e -iE (t+to) |fi)(fi|0) 
oo 

+ ^e- 4i5 " ( ' + ' o) |ni>(«i|0>. (193) 

n=£0 
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Then, we take the limit t —> oo(l — ie). This kills all the 
terms on the right hand side but the first one. Indeed, we 
have E n > Eq (n ^ 0) and this means that e~ lEot goes 
to zero less rapidly than all the other exponential terms. 
Therefore, in this limit, Eq. (193) reads e - lH( - t+t ^\0) ~ 
e~ iE °( t+t °'\Q)(£l\Q) and, by inverting this formula, this 
allows us to express the interacting vacuum as 



iE (t+t ) 

|ft) = lim , , ; 

' ' t->co(l-ie) (ft|0) 

e iE {t+t ) 

= lim — i ; — 

t^oo(l-ie) (ft|0) 



-iH(t+t ) 



(194) 



e -iH{t+t ) e iH (t+t )l 



e iE (t+t„) 

^-u ) -m~ Ui{t0, - m 



(195) 



In the second line, we have used the fact that H o \0) = 
since we can always define the ground state as the 
state associated to the zero of energy. This means 
that e lif °( t+t °)|0) = |0) and explains the appearance of 
the new factor e lH a(t+ta) _ j n ^he same fashion, using 



\n) 



JE (t-t ) 



|ft), one has 



(0|e 



-iH(t-ta) 



lim (0 1 ft) (ft I e" 

i— >oo(l— ie) 



-iE (t-t ) 



,(196) 



from which one deduces that 

„iEo(i-to) 

(ft| = lim (0|c 

t-s-oo(l-ie) (0|ft) 
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Therefore, the norm of the vacuum state is given by 

„2iE t 



(198) 
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t->-oo(l-«0 |(0|ft)| 

(199) 

In the following, we take (ft|ft) = 1 although this choice 
is not mandatory (but, then, when we compute a corre- 
lation function, it would be necessary to divide the cor- 
responding expression by (ft|ft)). 

We are now in a position where one can calculate the 
two-point correlation function of the field operator [recall 
that x(t) is directly proportional to the field operator, see 
Eq. (133)]. Using the previous results, one obtains 
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(202) 

(203) 

(204) 



where, in order to go from the first to the second line, 
we have used Eq. (199) and the normalization (ft|ft) = 1 
to express the term |(0|ft)| 2 . Eq. (204) is the Gell-Mann 
Law equation [16]. It will be the basis for our perturba- 
tivc treatment of quantum mechanics. It plays the same 
role as Eq. (98) since it allows us to perturbatively evalu- 
ate the various correlation functions in terms of Feynman 
diagrams (the Gell-Mann Law equation can be general- 
ized to higher correlation function if needed). 



C. Quantum Mechanics and Bubble Diagrams 

In order to mimic the self-interacting field theory stud- 
ied before, we now consider a quartic oscillator in quan- 
tum mechanics [55, 56], i.e. a system where the potential 
is given by V(x) ~ x 2 + x 4 . More precisely, this means 
that the interacting Hamiltonian can be expressed as 

#int = ^$ 4 , (205) 

where we recall that $ is given by Eq. (133). Let us 
then compute the two-point correlation function. For 
this purpose, as explained before, we use the Gell-Mann 



22 



Law equation. The numerator reads 



T 



$(t a )*(ti) exp (-i J drH, 



0), (206) 



where, for simplicity, we do not write the limits (but there 
will be reestablished when needed). It is also important 



to remember that the field arc written in the interact- 
ing picture (and, again for simplicity we do not write the 
subscript "I", they can be considered as free fields). Ob- 
viously, one cannot calculate exactly this quantity so we 
do it perturbatively by expanding the exponential. At 
first order in the coupling constant A, we obtain 



T $(t 2 )$(ti) 



dr$ 4 (r) 



0) = (0|T[*(t a )$(ti)]|0) 
A 

-*4! 



dr (0 |T [$(i 2 )$(ii)$ 4 (r)] | 0) + • • • , (207) 
D(h,t 2 ) - i A J dr (0 |T [^(ta)*^)* 4 ^)] | 0) + • • • . 



(208) 

It should now be clear that the same expression can also be written in terms on Feynman diagrams. One arrives at 
the following expression 

3 q 

hm mt 2 +— ti m V • t 2 +0 (A 2 ) , 

(209) 
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$(t 2 )$(*i) exp -i / &tH 1 



> = h 



where we have defined the following Feynman rules [55, 
56] 



D(t 1 ,t 2 ) = h • mh 



iX / dr 




(210) 



(211) 



Clearly, these Feynman rules are nothing but the one- 
dimensional version of the standard Feynman rules (106) 
and (107) of quantum field theory. The only difference is 
that four dimensional space-time integrations are simply 
replaced with one-dimensional time integration. In fact, 
in terms of diagrams, the expressions in field theory or in 
quantum mechanics are exactly similar. In particular, we 
see the appearance of diagrams with loops and external 
legs (the tadpole) and also of diagrams with loops but no 
external legs, our bubble diagrams. The advantage of our 
approach now becomes clear. We are in a position where 
not only one can compute these diagrams in quantum 
mechanics but we can also investigate their deep meaning 
by comparing the result with the standard perturbativc 
theory of quantum mechanics where everything is well 
under control. 

Let us start with the tadpole diagram. Using 
Eqs. (210) and (211), the corresponding expression can 



be written as 



-iX J dr D(ti,T)D(r,T )D(t, t 2 ) 



(212) 



-^3 J drc- M|tl - r| ^ |r -* 21 (213) 

*A _, , /, i 
-r^D (t 2 ,h )( t 3 -ti -- 



(214) 



The most striking difference with the case of quantum 
field theory is that the tadpole diagram is now finite. In 
some sense, this is the signal that ordinary quantum me- 
chanics does not need renormalization: the perturbativc 
approach leads to finite result. Technically, this is due to 
the fact that we deal with one-dimensional integrals that 
have better convergence properties. Physically, this is 
due to the fact that a single quantum-mechanical oscilla- 
tor leads to finite predictions while an infinite collections 
of such systems apparently lead to infinite physical quan- 
tities. 

Now let us calculate the "bubble diagram" present in 
the second term of the above expansion (209). Having es- 
tablished that the tadpole diagram is now finite, it is par- 
ticularly interesting to see what happens for this type of 
diagram. Using again the Feynman rules (210) and (211), 
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it is straightforward to establish that 

OO = -* A / D 2 {t,t)At (215) 

= J dT -°°- ( 216 ) 

Therefore, even in quantum mechanics, the bubble dia- 
gram remains divergent. This shows that the divergent 
nature of the loop and bubble diagrams is, in some sense, 
different. They are both divergent in quantum field the- 
ory but one (loop diagrams with external legs) becomes 



finite in the limit of quantum mechanics while the other 
(bubble diagram) remains infinite. 

But this now leads to the following question. We know 
from the standard perturbation theory of quantum me- 
chanics that all the physical predictions are finite, even 
without renormalization. However, using the other ap- 
proach based on the Feynman diagrams, we have just 
seen that it leads to terms that can be divergent. How 
these two facts can be consistent with each other? The 
answer is of course that, as in field theory, the bubble 
diagrams exactly cancel out in the calculation of the cor- 
relation functions. This can easily be seen if one notices 
that 



OO 



3 12 / 12 

t 2 + -7 h • • t 2 +-r h • — ¥ — • t 2 = U • • t 2 + —r ti • — ¥ — • t 2 
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(217) 



and that 



exp 



drH, 



(219) 



Therefore (as expected) the infinities of the bubble diagram never play a role in ordinary quantum mechanics. As a 
consequence, the two-point correlation function is given by 

(n|*(i 2 )*(ti)|ft) = ii« »t 2 + ^ h • mt 2 +0 (X 2 ), (220) 



which is perfectly finite. Again, in terms of Feynman 
diagrams the above equation is strictly equivalent to the 
corresponding one in quantum field theory, see Eq. (125). 



D. Quantum Mechanics and Vacuum Energy 

Let us now discuss how the vacuum energy is evalu- 
ated in ordinary quantum mechanics. The calculation is 
straightforward when there is no interaction and, from 
Eq. (139), one obtains 

(H) = -j • (221) 

This is of course nothing but the ground state energy 
of an harmonic oscillator. It cannot be zero because, 



due the Heisenberg uncertainty principle, the kinetic and 
potential energy cannot vanish at the same time. It can 
also be rewritten as 

(H)=lu 2 D(0), (222) 

where we have used Eq. (148). This exactly corresponds 
to Eq. (114) with m = uj and 4 — > 1 (since we now work 
in one dimension rather than in four), except that there 
is a sign difference which comes from the fact that, in 
field theory, we have used the signature (- + + +). Of 
course, the big difference is that -D(O) is a finite quantity 
in quantum mechanics. In some sense, the vacuum en- 
ergy problem appears to be very simple, almost trivial. 
The ground state of a quantum-mechanical oscillator is 
finite but a quantum field is a system which can be de- 
scribed by an infinite collection of harmonic oscillators 
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and, as a consequence, its ground state energy is infinite. 

Let us now consider a case where an interaction is 
present as we did in Sec. IV C. For simplicity, we consider- 
again the case of a quartic oscillator with Vint = A$ 4 /4!. 
Upon using the standard perturbation theory in quantum 
mechanics, one can calculate the displacement of the en- 
ergy levels. They are given by the standard formula [57] 



E, n — E„ 



(n\V ir 



(223) 



where \n) are the unperturbed eigenvectors. It is easy to 
show that they can be expressed as 



AE n = E n — E n 



A 



32w 2 



(2n 2 + 2n + l) , (224) 



which means that, for the first two levels, the shift in 
energy is given by 



AE 



A 



32w 2 



AEi = - 



5A 
32w 2 



(225) 



We can now use the other method, based on the Feynman 
diagrams. It is easy to obtain that 



A^o = ^ 4 > = !<* 2 > 2 = ^ 2 (0) 



8/di 



oo 



(226) 



Clearly, we have just obtained the strict equivalent of 
Eq. (124), the term J d A x being replaced by J dt as could 
have been guessed in the context of quantum mechanics. 
Then, using the fact that D(0) = l/(2w), one obtains 
that 



32w 2 



(227) 



that is to say we recover the result established before, see 
Eq. (225). Therefore, we have shown explicitly that the 
two approaches lead to the same result. It is also inter- 
esting to notice that AEq is given by a bubble diagram 
but that this divergent graph is canceled by the infinite 
term J dt so that a finite quantity is left. 

The previous considerations also lead to the following 
question: before we had to deal with the renormalization 
of the mass. How does this aspect appear in the present 
context? In this case, we deal with the same diagrams, 
namely 



1 



(n\T[<f>(ti)<f>(t 2 )]\Cl) = x x m «x 2 +-xi« — ¥ — *X2+..., (228) 

or, in terms of explicit expressions [compare with Eqs. (125) and (126)] 

(fi|T[0(ii)0(i 2 )]|fi) = Bp(*i-*2)-y j drD(t\ — t)D f (0)D f (t t 2 ) (229) 



iX If e~ iB(tl "* 2) 
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r 



and this leads to 

<4» = ^ 2 + ^(0), (231) 

which is exactly Eq. (130) (again, there is a sign difference 
which originates from the sign difference in the propaga- 
tor, see above). In quantum mechanics, the propagator 
is finite and therefore 

^ren = ^ + ^, (232) 

or 



Then, one can express the spectrum in terms of the renor- 
malized "mass" w ren as we have done in quantum field 
theory although this is not necessary here since the result 
is anyway finite. Straightforward calculations lead to 

E n = uj ien (n + ( 2 « 2 - 2n - 1) . (234) 

It does not come as a surprised that, as a function of the 
renormalizcd mass, the spectrum is a different function 
of n, see Eq. (224). However, from the above expression 
one also obtains 

E.-Eo^ uw (235) 



A 



(233) 



In presence of an interaction, the difference between the 
first excited state and the fundamental level has been 
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"modified" (or, rather, is expressed differently). But, 
observationally, we define or measure the spectrum with 
respect to E-y — Eq. In other words what we define as u> 
is in fact E\ — Eq. Therefore, renormalizing u> consists 
in preserving this definition even in presence of an inter- 
action. Obviously, E\ — Eq is not the bare uj (the one 
which appears in the Lagrangian) when the an-harmonic 
term is present. What we do here is to define w as the 
observed one, namely as E\ — Eo, and then calculate the 
spectrum in terms of that observational quantity. What 
is done in quantum field theory is clearly exactly similar 
to the above described procedure except of course that 
the bare quantities are in fact infinite because D F (0) is in- 
finite. Therefore, this exercise illustrates nicely the deep 
meaning of renormalization. 

This concludes our discussion of the vacuum energy 
density in terms of Feynman diagrams. We have shown 
that the cosmological constant is in fact given in terms 
of very peculiar graphs, the so-called bubble diagrams. 
These diagrams are, in a sense, more divergent than the 
usual loop diagrams because they remain infinite even 
in the limit of quantum mechanics contrary to the last 
ones. In non gravitational physics, this is not a prob- 
lem because the bubble diagrams always cancel out in 
the equations describing observable quantities. When 
the gravitational field is turned on, we face the tasks 
of renormalizing these diagrams which is more difficult 
than in the standard situation because of their bad be- 
havior mentioned before. Given that astrophysical ob- 
servations seem to indicate that the vacuum energy is 
non vanishing (see below), the details of the renormaliz- 
ing procedure become a crucial issue. Not only we have 
to extract a finite quantity from a divergent graph but 
this finite quantity must be in agreement with the obser- 
vations. As argued before, the technique which consists 
in imposing a spatial cut-off is certainly not consistent 
(and, as well-known, gives a very large contribution to the 
vacuum energy). A more reasonable approach seems to 
be dimensional regularization since this satisfies Lorentz 
invariance and is therefore consistent with the vacuum 
equation of state. As shown before, this typically leads 
to Eq. (96). In the next section, we present a new argu- 
ment, based on the Gaussian effective potential, which 
seems to support this conclusion. 



VI. THE GAUSSIAN EFFECTIVE POTENTIAL 

The Gaussian effective potential is a non perturbative 
approach to quantum field theory [58-66]. Originally 
developed in the context of quantum mechanics, it is 
has been generalized to field theory. Much less known 
is the fact that the Gaussian effective potential method 
has something to say about the vacuum energy problem 
and the aim of this section is to present these consid- 
erations [59]. In the next subsection, we quickly recall 
the main idea and how this approach can be applied to 
ordinary quantum mechanics. Then, we explain how it 



can be implemented in quantum field theory and, finally, 
we discuss its implication for the cosmological constant 
problem. 



A. The Gaussian Effective Potential in Quantum 
Mechanics 

The main idea underlying the Gaussian effective poten- 
tial approach is the following one [58]. At the classical 
level, a system is described by a potential. Once this one 
is specified, one can compute the evolution of the sys- 
tem. When quantum effects are taken into account, the 
behavior of the system will be modified and, by means of 
the formalism of quantum mechanics, one can work out 
the corresponding physical predictions. The idea of the 
Gaussian effective potential is to find an effective poten- 
tial which, at the same time, can be used as if it were 
a "classical" potential (that is to say can be used in the 
framework of Newtonian dynamics) and takes into ac- 
count the quantum effects. At first sight, this approach 
seems equivalent to what is known in the literature as 
the effective potential [67, 68]. In this method, one re- 
quires the wave-function mean value to be centered at 
a fixed point $o- Then, the effective potential gives the 
exact ground state energy of the system in presence of 
quantum corrections. However, in some situation, the 
result can be quite artificial. Consider for instance a sit- 
uation where the wave-function possesses two large peaks 
on each side of $o- In such a case, the effective potential 
does not necessarily lead to a correct description of what 
happens in the vicinity of $o because its shape could 
just be determined by some average procedure between 
the two distant peaks. Typically, one would obtain a con- 
vex potential centered at $o while our physical intuition 
rather indicates that the effective potential should be a 
double- well potential, the two wells corresponding to the 
two peaks of the wave-function. Another problem is that 
the standard effective potential is usually computed per- 
turbativcly in h. Clearly, when the quantum corrections 
become important, this approach breaks down. Since the 
Gaussian effective potential does not suffer from the two 
above mentioned drawbacks, we will use it in this review. 
When possible, we will discuss the difference between the 
two approaches in order to stress their differences and 
similarities [69]. 

We now turn to the precise definition of the Gaussian 
effective potential. This definition is based on a vari- 
ational principle. Here we consider a one-dimensional 
system as we have already done in Sec. VA. In particu- 
lar, this means that the scalar field is given by Eq. (133). 
Then, the Gaussian effective potential is defined by [58] 

V G ($ ) = nun (* \H\ , (236) 

where H is the Hamiltonian of the system and where the 
wave-function is taken to be a simple Gaussian (hence 
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the name Gaussian effective potential) 

(SlY 4 e -n(*-* ) 2 /(2ft). (237) 

The parameter O is related to the width of the wave- 
function. Here the idea is to find the value of f2 which 
minimizes the energy of the system. In practice, the 
mean value of the Hamiltonian that appears in Eq. (236) 
can be expressed as 



(H) = 



d$#* ($) 



h 2 d 2 

2 d$2 



(238) 

To illustrate how the calculation of the Gaussian effective 
potential works on a concrete example, we choose the 
following potential 



v ; 2 



A$ 4 
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that is to say the potential of an an-harmonic oscillator. 
Another motivation for this choice is that the compari- 
son with field theory is easy since the prototypical model 
possesses the same potential. Upon using Eq. (237) in 
Eq. (238), one arrives at 
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We must now determine the optimal O. As explained 
above, this is achieved by requiring d (H) /<9Q = 0, which 
leads to the equation 



O 3 - 



12A* 2 ,) 



n - 6h\ = o. 



(241) 



This equation allows us to determine Q in terms of $o ■ In 
the following, we write this result as f2o = f2($o)- Upon 
using the above relation in the expression of (H), one 
obtains the following expression for the Gaussian effective 
potential 



V G (So) = \m 2 <S>l 



2 



3fi 2 A 

4O ' 



(242) 



Of course, as already emphasized before, one should not 
forget that Oo is a function of $o- 

A complete study of the above potential has been pro- 
vided in Ref. [58] and we will not repeat this analysis. 
Here, we just want to illustrate that the Gaussian ef- 
fective potential is an accurate method to estimate the 
ground state of a system (even in the strong coupling 
limit) since this is directly relevant for the cosmologi- 
cal constant problem. Let us consider the case where 
m 2 > (the case of the double-well potential, reminis- 
cent of symmetry breaking in field theory has been stud- 
ied in Ref. [58]). The ground state of the an-harmonic 
oscillator is given by Eq. (242) 
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where Oo is solution of 

Q 3 - m 2 fl - 6h\ = 0. 



(244) 



The discriminant of this cubic equation can be expressed 
as 



A ($ = 0) = 4m b 1 - 



972 1 

"16" i 5 



(245) 



where the quantity £ in the above relation is given by 

™2 



(2A) 2 / 3 



(246) 



The awkward coefficient 972/16 in Eq. (245) originates 
from our will to work with the definition of £ used in 
Ref. [58]. The cubic equation (244) can be explicitly 
solved. We deal with different branches according to 
the sign of the discriminant. This sign is determined 
by 6im = (972/16) 1 /e ~ 1.9827. If £ < £iim (£ = repre- 
sent the quartic oscillator) we are in the strong coupling 
regime which is of particular interest in order to demon- 
strate the usefulness of the Gaussian effective potential. 
Therefore, let us focus on this situation. In this case, the 
exact solution of Eq. (244) can be written as 
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Together with Eq. (243), the above equation gives an ex- 
cellent approximation for the ground state. For instance, 
if to = 1 and A = 10, then £ 2 ~ (1/20) 2 / 3 ~ 0.1357 and 
flo = 4. As a consequence, this gives V G ($o = 0) — 1-53, 
a result that is accurate at ~ 1.75% according to Ref. [58]. 

This simple example has shown that the Gaussian ef- 
fective potential can be an efficient tool to calculate the 
ground state of a quantum system even in a non pertur- 
bative regime. Therefore, it appears as an interesting ap- 
proach for the vacuum energy problem. However, before 
turning to the calculation of the Gaussian effective po- 
tential in field theory, it is interesting to compare it with 
the standard effective potential in more detail. This is 
the purpose of the next subsection. 



B. Comparison with the One Loop Effective 
Potential 

In order to explain how the effective potential is ob- 
tained, we need to quickly return to the basics of quan- 
tum field theory. In Eq. (97), we have defined the gener- 
ating functional Z[J\. Let us now define another gener- 
ating functional W\J] by [16, 17] 



Z[J] 



(248) 
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Following the procedure in Eq. (98), one can define new 
correlation functions according to 



5 n W [J] 



5J(xi) ■ ■ ■SJ(x n ) 



(249) 
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As is well-known, these functions represent in fact the 
connected n-points functions. Then, one introduces the 
classical field $ r 



SW[J] 



8J(x) 



(250) 
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where we remind that J(x) is a source. This definition 
can be easily justified. Indeed, from Eqs. (97) and (248), 
one sees that 
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Therefore, we see that the classical field is in fact the vac- 
uum expectation value of the field operator in the vacuum 
state of the theory. When the source vanishes, the clas- 
sical field does not necessarily goes to zero, a typical and 
well-known example being of course symmetry breaking. 

The next step consists in introducing a new generating 
function, the so-called effective action, defined by 



r[$ c ] = w[J] - J & A xJ{x)<s> c {x) 



(253) 



Clearly, r[$ c ] does not depend on the source since 
5T/5J = 0. Let us now calculate the functional deriva- 
tive of r with respect with the classical field. This gives 
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In order to interpret these equations and to better un- 
derstand their meaning, it is interesting to apply them to 
the case of a free theory where all the calculations can be 
carried out explicitly. From Eq. (100), we immediately 
see that (the subscript zero indicates that the calculation 
is performed for the free theory) 

*Wq[J] =-\j d 4 xd 4 yJ(x)D F (x - y)J(y). (256) 

Upon using Eq. (250), one obtains 

M z ) = ~\ J d 4 xD F {z-x)J{x). (257) 



But iD F is the Green function of the Klein-Gordon oper- 
ator. As a consequence, the above equation implies that 



(ifd^-m 2 ) $ c (x) = J(x), 



(258) 



i.e. the classical field obeys the classical equation of mo- 
tion (hence its name). Then, using Eqs. (256) and (257) 
into the definition (253), the effective action can be ex- 
pressed as 
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d 4 xd 4 yJ(x)D F (x-y)J(y). (259) 



The next step is to use the equation of motion of the 
classical field (258) to rewrite the effective action as 
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where d^' y means a partial derivative with respect to the 
coordinates at point x and y. Then, integrating by part 
the x dependent terms in the above integral leads to 
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Finally, using the fact that D F /i is the Green function 
satisfying —if* 'd ii d v D F + m 2 D F = S, one obtains 

r [$c] = \S d 4 xM?) (rTW v - m 2 ) Mx) 

= -\f d 4 x(7 7 ^a Al $ c a^ c + m 2 $ c ),(262) 

and one recovers the action of a free scalar field theory. 
In presence of interactions, r[$ c ] will no longer agree 
with the classical action. The quantum corrections will 
transform it into a complicated non local functional. In 
this situation, one typically expects an expression that 
can be written as 



r[#c] = -/< 



^A(^ c )l,^d^ c d v ^ c 
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where A($> c ) and B(Q C ) are functions (not functional) 
of $ c and the dots represent higher derivative terms. 
Clearly, the previous considerations justify the name "ef- 
fective action" for r[$ c ]. 

Another interesting aspect is that the effective action 
is also the generating functional of the irreducible corre- 
lation functions rW (also known as proper vertices). In 
other words, one has 
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Equivalcntly, one can express the effective action as a 
Volterra expansion, namely 



r[$ c 
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d 4 x 1 ---d 4 x n T^{x 1 ,--- ,x n ) 
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which is nothing but another way to write the func- 
tional (263). 

Now, suppose that we are interested in determining 
the effective action concretely. It is clear that this is 
a complicated calculation. As a first step, it is worth 
calculating the effective potential in Eq. (263). It is true 
that it does not represent all the information contained in 
the effective action but this would allow us to determine 
the location of the minimum of the system in presence 
of quantum corrections. Therefore, this is an interesting 
quantity. For this purpose, it is sufficient to consider a 
case where $ c is constant. Indeed, in this situation, the 
effective action reduces to T = — VV c s(^ c ), where V is 
the space-time volume. On the other hand, the Volterra 
expansion can be re-expressed as 



r[$ c 



71=0 



(266) 



space of the exponentials lead to Dirac functions in mo- 
mentum. As a consequence, one arrives at the following 
expression 

oo 1 

r[$c] = ^^-r (n) (fc l = o)$™, (269) 

TV. 

where, as already mentioned, V — (2ir) 4 S(0) is the space- 
time volume. One finally reaches the result that 

Kfr (0c) = - V ^r (n) (h = o)$?. (270) 

n=0 

Therefore, in order to calculate the effective potential, 
one just has to evaluate the proper vertices with vanish- 
ing momenta. At one loop, this gives rise to the following 



sum 



-=cr + a 



(271) 



It turns out that this series can be summed up as we are 
now going to show. 

In order to see how it works, let us calculate the second 
diagram in the above expansion. For this purpose, we 
start with the following diagram 



This expression can be further simplified if one takes the 
Fourier transform of the proper vertex, 



^ (^1 , • • • , x n ) 



d 4 fci d 4 fc» 
(2tt) 4 " ' (2^)* 



^ {ki , • • • , k n ) 



xe~ ikl ' Xl ■ ■ ■ Q~ ikn ' x " 



(267) 



In fact, for convenience, we introduce the coefficients 

— (n) 

r (fci, ■ ■ • ,k n ) that are easier to manipulate. They are 
defined by the following expression 

rW (fcx, • • • , k n ) = (2tt) 4 s (fcj + • • • kn) r (n) (ki, ■■■k n ), 

(268) 

Then, we insert the Fourier expansion (267) into the ex- 
pression (266) of the effective action. The integrals over 



V 



(272) 



According to the Feynman rules (106) and (107), this 
diagram is equivalent to the expression 

V = {-i\f J dudvD(t x -u)D(t 2 -u)D' 2 {u-v) 

xD(v -t 3 )D(v-t 4 ), (273) 

where u and v denotes the two internal vertices. Upon us- 
ing Eq. (155), this expression takes the form [here, we use 
the quantum- mechanical theory introduced in Sec. VA 
where the "mass" is denoted ui) 



^ 2 / , , / d£i dE 2 dE 3 dE 4 dE dE ie^E^-u) ie -tE 2 (t 2 -u) ie -iE 3 (v-t 3 ) 
V = \—iX) I dudv 



IC 



El 



E 2 



2tt 2tt 2tt 2tt 2tt 2tt 

IE(u-v) ^~ — iE(u—v) 

1 E 2 -lu 2 



El 



El 



El 



(274) 



The integrals over u and v are easily performed and leads to two Dirac functions. Then, a further integration over E 
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can be performed and we are left with the expression 

T) ( ^fO ^XfJP MJP P TP \ f AE l d ^ 2 dE3 dE * dE ie ~ lEltl ie ~ lE2t2 te ~ lE3t3 

V = (-iX) (2tt) S (E 1 + E 2 - E 3 - E 4 ) I — — »~Ev 2~?2 2^2 2 

J 2ir 2tt 2tt 2tt 2tt E( — uj 2 E£ — uj 2 E 2 — uj 2 E£ - uj 2 

i i 

-. (275) 



E 2 - uj 2 (E 3 + E±- E) 2 - lu- 

In order to compute the expression (271) (or rather the contribution of the second diagram to this expression), we 
have seen that we must amputate the diagram by appending the external leg propagators. This amounts to ignore 
the factors i/(E 2 — uj 2 ) in the above expression. Moreover, since the graph is given by the Fourier transform of T>, 
this means that only the expression inside the above quadruple integrals over Ei should be considered. Finally, the 
"momentum" (in the present context the energy) of the external legs should be put to zero (i.e. we require all the 
Ei's left in the expression of T> to vanish) and we must factorize out a Dirac function expressing energy conservation 
[times (2tt) 2 since we are considering a + 1 theory]. This means that Eq. (275) takes the following form 

^~>^ , .,,2 fdE 



= <- iA) ' liW^W^' (276) 



The fact that the propagator appears squared is of course 
not fortuitous. It comes from the fact that we have two 
"double legs" attached to the loop. For n "double legs" , 
it is clear that we would have the propagator to the power 
n and an overall factor (— i\) n . As a consequence, at one 
loop, Eq. (270) can be expressed as 



00 



E 

P =i 



$2p 

2p 



A4! 



dE 



1 



2tt [E 2 - lu 2 )p 



(277) 



where we have taken into account a symmetry factor and 
where we have noted n = 2p, p being the number of 
external points. The sum can now be performed and we 
obtain 



Veg ($ c ) 



2tt 



12A$2 

E 2 - uj 2 



~ {-m + 



12A$; 



(278) 
(279) 
(280) 



where we have re-established h and where we have taken 
into account the sign difference for uj (or m) due to the 
fact that we work with a ( — h ++) signature, see the 
remark after Eq. (222). This expression should be com- 
pared to Eq. (242). In fact the first term just comes 
from the fact that the method of the effective potential 
assumes that V e g (0) = 0. Therefore, in order to compare 
with the Gaussian effective potential, one must ignore 
this rescaling. As a consequence, the total potential (in- 
cluding the leading order in h°) reads 



V eS ($ c ) 



1 



ym 2 + 12A$2. 



(281) 



One sees that the vacuum energy is predicted to be 
K-ff(O) = hm/2 a result which is not realistic in the non 
perturbative regime. 



In fact the comparison can be made in an even more 
explicit manner if one remembers that the one loop Gaus- 
sian effective potential consists in retaining the dominant 
term is a systematic expansion in h. As a consequence, 
one should ignore the h 2 term in Eq. (240) and the h term 
in Eq. (241). In particular, this last approximation leads 
to the following explicit solution for the gap equation 



n 2 



m 2 + 12X^1 



d 2 V 



(282) 



If we then use Eq. (242) and neglects the h 2 term in that 
equation, one obtains 




o(h 2 ) 



+ 12A* 2 , 



(283) 
, (284) 



that is to say exactly Eq. (281). Wc conclude that the 
effective potential is in fact contained in the Gaussian 
effective potential approach. It just consists in neglecting 
higher power of H since this is a one loop approach. As a 
consequence, the effective potential cannot be reliable in 
the non perturbative regime where the quantum effects 
are strong. 

Having demonstrated that the Gaussian effective po- 
tential approach is an efficient tool, we now need to study 
how it works in quantum field theory. Then, we will able 
to use it in order to calculate the energy density of the 
vacuum state. 



C. The Gaussian Effective Potential in Field 
Theory and the Vacuum Energy 

We now turn to the main question of this section, 
namely the calculation of the Gaussian effective potential 
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in quantum field theory [59, 60, 63, 65]. For this purpose, 
we now consider a model similar to the one considered 
in the previous subsections, namely a free model plus a 
quartic self interaction 

U = \® 2 + 5^'**^* + \ m l<S> 2 + A B $ 4 , (285) 

where m B and A B are the bare mass and coupling con- 
stant, respectively. Then, we write the field operator as 
(we work in d space-time dimensions) 



$(t,as) = $ n 



1 



d d - 1 fe 



c k e 



—icjt-\-ik-x 



with 



(286) 



(287) 



The quantity $ represents a classical and constant field. 
The parameter SI denotes the mass of the excitations 
around this classical value. The idea at the basis of 
the calculation is the same as before, namely minimiz- 
ing the energy with respect to the parameter tt. Using 
the properties of the creation and annihilation operators, 
it is straightforward to calculate the mean value of the 
Hamiltonian in the vacuum state. We find 



n 2 



A B *$ + 6A B *gJ (fi) 



Jo(n) + im£$g 
f 3A B / 2 (S1), (288) 



where we have defined the integrals I p (Q) by 
1 f d^k 



i P (n) 



(27T) 



d-l 



2wo(fc) n 1 ' 



(289) 



Of course, the integrals I P (SY) can be divergent and this 
signals, as usual in field theory, the need for renormal- 
ization. In fact, these integrals can be performed exactly 
and one is led to 



1 r[(2-d-2p)/2] 



2(47r)( d " 



l)/2 



r[(l - 2p)/2] 



d+2p-2 



(290) 



It is clear that, for some values of d and p, the above 
expression does not exist and is only formal. 

As already mentioned, in order to calculate the Gaus- 
sian effective potential, one must vary (H) with respect 
to f2. Upon using Eq. (288) and dI p /dQ = (2j5-l)ftl p _i, 
one arrives at the following "gap" equation, the field the- 
ory counterpart of Eq. (241) 



12A B [$g+i (n)] 



(291) 



This is an algebraic equation for f2. Its solution can be 
written as f2 = ($o;m B , A B ), that is to say it depends 
on $o and on the bare parameters (in the following we 



no longer write the dependence on the bare parameters) . 
Then, the Gaussian effective potential is defined by in- 
serting this solution into Eq. (288). This leads to 



V a (^ ) = (H) [fi(* ),$o]. 



(292) 



Upon using the gap equation, one can also find a more 
compact expression, namely 

U G ($o) = h PI(* )] - 3A B / 2 [H($ )] + Im B $ 2 

(293) 



At this level, this expression is in fact still formal be- 
cause of the divergences of the integrals I p . Therefore, 
we must first tame those divergences. This is done by 
means of a renormalization of the parameters appearing 
in the Lagrangian of the model. 

Renormalization consists in parametrizing the physi- 
cal quantities in terms of new, observable parameters, 
different from the bare parameters that appear in the 
Lagrangian. We now define the renormalizcd mass by 



d 2 V r . 



d$ 2 



(294) 



*n=0 



Upon using Eq. (293) and the following relation [which 
comes from differentiating Eq. (291)] 



12A„ 



d$ n i + 6A B j_i(n)' 



(295) 



one finds that 
d 2 K 



d$ 2 



i\ + i2A B [$2 + i (n)] 

(12A B $ ) 2 I_i(n) 
l + 6A B /_i(fi) 



(296) 



As a consequence, putting $o — hi the above equation, 
the renormahzed mass can be expressed as 



to 2 =m% + 12A B / (^o), 



(297) 



where f2o denotes the solution of the gap equation (291) 
when $ = 0, i.e. O 2 = to 2 + 12A B Jo (fi)- It is obvious 
that Slo is not a function of $o but depends only on ra B 
and A B . In fact, this can also been seen from Eq. (291) 

since it tells us that m 2 = O . As a consequence, the 
above equation can be re-written as 

ml =ml -12A B 7 (mJ. (298) 

We are now in a position where we can calculate the 
Gaussian effective potential in terms of the renormahzed 
mass. For this purpose, we restart from Eq. (292) and 
substitute the expression (298). One obtains 



u G ($ ) = h (n) + 2 (< - n J Mn) + 2< $ o 

+A b $4 + 6Ab [j (n)-/ (m R )] ®l 
+3A B II (O) - 6A B h (O) h (m R ) . (299) 
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However, at this stage, the expression of the Gaussian ef- 
fective potential is not yet a function of the renormalized 
mass only as one still needs to express Io(p) and li(ft) 
in terms of Jo(m R ) and ii(m R ). There are several ways 
to do it, see Ref. [59]. Here we choose to use a method 
where a cut-off is temporarily introduced. We work in 
d = 4 dimensions to simplify the calculations. Let us 
exemplify the method with Iq. We have 



2b(fi) 



1 



2(2tt) e 
1 

4^2 



d^fcej 



12 



dkk 2 k 2 



n 2 



-1/2 



(300) 



(301) 



where we have performed the angular integrals and where 
we have introduced a cut-off M to regulate the otherwise 
divergent integral. This gives 



1 

4tt 2 



M 2 



ft 2 
2 



In (m + V^+AP) + 



As a consequence, one can write 



<r 2 



In ft 



(302) 



/„(«)- /o(m) 



(m 2 -f7 2 )/_i(m) + 



ft 2 m 2 
In ft In m 



4tt 2 



M 



71/ 



Vfi 2 + A/ 2 -— \/m 2 + M 2 + 



- (m 2 - ft 2 ) 



M 



\Jm 2 



M 2 



In 5 



ft 2 



In 



M + \/m 2 + M 2 



m + Vft 2 + m 2 



(303) 



r 



where we have used 
1 



/-i(ft) 



4tt 2 



-M 



Vft 2 + M 2 



In ft 



In (M + Vft 2 + M 2 ) 



(304) 



an expression which can be obtained with the same 
method as described above. Then we take the limit 
M — > +oo and Eq. (303) becomes 



/o(ft)-/ (m) 



1 



ft 2 



m 2 ) I-i(m) 



+ [x]nx-(x-l)\, (305) 
lbn z 



where x = ft 2 /m 2 . Notice that, at this stage of the 
calculation, there is no need to write x in terms of ft 
and/or to r . The above expression is valid for any ft and 
to. Following Ref. [59], we call the function between the 
square brackets in the above equation L%(x). Exactly in 
the same manner, one can also show that 



/x(ft) - h{m) = - (ft 2 - to 2 ) J (m) 



ft 2 



to 2 ) 2 /-i (to) + 



32tt 



:L 3 (x), (306) 



with L 3 (x) = [2x 2 \nx - 2(x - 1) - 3 (a; - l) 2 ]/4. Then, 
upon using Eqs. (305) and (306) in Eq. (299), one obtains 



y G ($o) = h (m R ) - 3A B / 2 (to r ) + -to 2 * 2 + X^ 4 + ^L 3 (x) + -X B m R (x - l) 2 1\ (to r ) 



+ -/_ 1 (m R )m R (z-l) (x-1) 



3A3 
2tt 2 



' 16tt 2 2 

L 2 (x) + 



£ 2 (o0 + 327r 2 -^f 
mi 



(307) 



with x = ft /to 2 . Finally, the gap equation (291) can 
also be re- written in terms of the renormalized mass. The 



corresponding expression reads 
(S-l)[l + 6A B /_ 1 (m a )] 



3A3 

47T 2 ^ 
+ 167T 2 



L 2 (x) 
$ 2 1 



(308) 
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We see that the dependence on <f>o of the Gaussian effec- 
tive potential is quite complicated. Indeed, solving the 
gap equation provides us with x = x($o) that we must 
then insert in Eq. (307). 

We shall now discuss the renormalization of the cou- 
pling constant A B . It is defined by 



A„ = 



1 d 4 K 



4! d$ 4 



(309) 



Proceeding in the same manner as before, one obtains 
the following equation for the quantity A R 



A, 



A, 



1 - 12A B J_ 1 (to r ) 
1 + 6A B /_i (m R ) ' 



(310) 



The renormalization of the coupling constant has been 
discussed at length in Refs. [59, 63] and we will not re- 
peat all the details since, here, our main concern is the 
vacuum energy. In fact we will focus on the case named 
"precarious" in Ref. [59]. Let us return to Eq. (310); one 
can invert this expression and express the bare coupling 
constant in terms of the rcnormalizcd one. One obtains 



A, : 



1 - 6A R 7_i (to r ) 




48A R 7_ 1 (m R ) 
[l-eA^rKJf 



Then, from Eq. (304), we see that J_i (to r ) is in fact a 
logarithmically divergent quantity when the cut-off M is 
sent to infinity. Therefore, we must treat J_i (to r ) as an 
arbitrarily large quantity. In this limit, if one chooses the 
plus sign in Eq. (311), one obtains 



A, 



i 



J_i (m R ) 



(312) 



As shown in Ref. [59] , this case is not viable. Intuitively, 
it is clear that a large negative coupling constant leads to 
a potential which is not bounded from below. Therefore, 
we are left with the minus sign for which one finds that 



A, : 




(313) 



In passing, it is worth signaling that this equation can 
be re-expressed in a different manner as follows. Let us 
first introduce the new characteristic scale M c (not to be 
confused with the cut-off M discussed before) defined by 



1 



1 



A R ~ 4tt 2 
Then, if we use the fact that 



In 



M c 



/_i(M c )-J_x(m R ) 



1 

8^ 



(314) 



(315) 



where L\(x) — lnx [an expression which can be de- 
rived with the same method that has led to Eqs. (303) 
and (306)], the bare coupling constant can be re- written 

as 



K = - 



1 



6/_i(M ) 



(316) 



Therefore, it has a negative but infinitesimal value which 
implies that the stability properties of the corresponding 
scenario are much better compared to the case (312). As 
argued in Ref. [59] , besides simplifying the expression of 
the renormalized coupling constant, the scale M c plays 
an important role in a complete study of the precarious 
case (even if, in the present context, this is not apparent 
since we focus on the vacuum energy only). Inserting the 
expression (313) [or Eq. (316)] into Eq. (307), one can 
re-write the effective potential as 



h (to r ) - 3A B / 2 (m R ) + -m 2 ^ 2 , 

4 4 
TO , s TO , N , 

+—*-L 3 (x) s- (x- 1 2 



16A„ 



-O 



I-i(m R ) 



(311) while the gap equation (308) becomes 



4tt 2 



L 2 (x) + 16tt 



<J>2 

2 v 



(317) 



(318) 



Let us notice that, in Eq. (317), we have left the quantity 
A B into the field independent term I\ (to r ) — 3A B /g (m R ) 
since it represents the zero point energy that will be 
treated in detail in the next sub-section VI D. Finally, 
if one uses the explicit form of the functions L 2 and L3, 
then the Gaussian effective potential takes the form 

^ G (d> ) = h (mj - 3A B / 2 (to r ) + l -mlx$>l 



1287T 2 

4tt 2 



2x 2 hix - 2(x - 1) - 3(x - l) 2 



and the gap equation can be re-expressed as 



(x-1) 1 + 



4^ 
A„ 



- 16tt" 



r ln" 



(319) 



(320) 



The two equations (319) and (320) give a complete ex- 
pression of the Gaussian effective potential for the renor- 
malized precarious "A<& 4 " theory [59, 63]. 

As we have done for the quantum mechanics case, it 
is interesting to compare the above result to the stan- 
dard effective potential approach. We know that this 
potential can be obtained from the Gaussian effective 
potential method by neglecting the terms proportional 
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to h 2 . Therefore, we restart from Eq. (288) and neglect 
the term 3A b /q since this corresponds to the standard 
one loop method (i.e. it is of order h 2 ). One obtains 

*o) = im + \{m 2 B ^n 2 )i {n) + l -m 2 ^l 

+ \ B $ 4 + 6KHlo(n) + 0(h 2 ). (321) 

As a consequence, this leads to the new following gap 
equation 

n 2 = ml + 12A B $^. (322) 
Clearly, this relation is the exact counter-part of 



Eq. (282). Then, one has to calculate the renormalized 
mass and coupling constant, one arrives at 

ml = m B + 12A B / (m B ), (323) 
A R = A B [l-18A B /_i(m B )]. (324) 

These equations should be compared to Eqs. (298) 
and (310). We notice that it is now m B that appears 
in the argument of the integrals rather than m R as in 
Eqs. (298) and (310). This is because Eq. (322) implies 
that r2($o = 0) = m B and not fi($o = 0) = m R as be- 
fore. The final step consists in using the two previous 
equation in Eq. (321). One obtains 



V G ($ ) = Mm B ) + 



-ml^l 



-216A B 4 

B m 4 



64tt 2 



1 + 12A, 



$0 



In 1 + 12A, 



12A — 

mt 



(325) 



the dots corresponding to the higher order terms in h. If 
one neglects those terms, then V G reduces to V c g. The 
corresponding expression is nothing but the Coleman- 
Weinberg potential. The previous analysis confirms that 
the effective potential is indeed a particular case of the 
Gaussian effective potential approach. An even more de- 
tailed comparisons of the two methods can be found in 
Refs. [59, 63]. Here, we do not pursue this issue and now 
come to the question of the vacuum energy problem. 



D. Vacuum Energy and Gaussian Effective 
Potential 



We are now in a position where the main problem of 
this section can be addressed. From Eqs. (307) and (317), 
we see that the vacuum energy density predicted by the 
Gaussian effective potential method (i.e. the term left 
after having taken 4>o = in those expressions) can be 
expressed as 



h(m R ) - 3A B I (m R ). 



(326) 



As we have already stressed many times, it is impor- 
tant to remember that this result is a priori non pcrtur- 
bative. Using the expression of the precarious coupling 
constant (313), one obtains 



h(m R ) 

1 4K) 

2/_ 1 (m R ) 
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2A/_i(m B ) 



Jo (mj. (327) 



The expressions of the integrals I p have been given in 
Eqs. (289) and (290). The three integrals that appear in 



the above expression are all divergent. Nevertheless, the 
ratio Io/I-i is finite and reads 



(328) 



Nevertheless, even after having inserted the above ra- 
tio into Eq. (327), this expression still contains singular 
terms. Explicitly, one has 



4(2 - d) 
1 



d-4 



R 20F(4^)( d - 1 )/ 2 



X R (2-d) 



(329) 



However, writing this expression for d = 4 + e and taking 
the limit e — > 0, we have the remarkable result that a 
second cancellation occurs and that the final expression 
is in fact finite. It can be expressed as 



1287T 2 
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/ m 2 
21n — | 
\M 2 



(330) 



where we have used Eq. (316). Even more remarkable is 
that the above equation is very similar to Eq. (96). In 
particular, the vacuum energy is proportional to fourth 
power, not of the cut-off, but of the mass of the cor- 
responding particle and the scale M c is present as a 
logarithmic correction. Therefore, we conclude that our 
analysis based on the Gaussian effective potential has to- 
tally confirmed the results presented in the previous sec- 
tions. We end this section by noticing that the vacuum 
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energy predicted by the Coleman Weinberg potential is, 
see Eq. (325), p vac = Zi(m B ) and is therefore different 
from the one obtained before, see Eq. (330). Moreover, 
it seems that this expression remains divergent. 



VII. CONTRIBUTION FROM OTHER FIELDS 

So far, for simplicity, we have only treated the case of 
a scalar field. However, in order to be more realistic, it 
is clear that one must also evaluate the vacuum energy 
for other type of fields. This question is the purpose of 
the present section. As we will see, the corresponding 
calculations are in fact very similar. In fact, the impor- 
tant question is the dependence of the overall coefficient 
which multiplies p mc with the spin of the particle. As 
is well-known, this is at the origin of a remarkable can- 
cellation which leads to the concept of super-symmetry. 
Therefore, in the following, we will pay special attention 
to this issue. 



A. Fermion Fields 

We start with the case of a Dirac spinor field describing 
a spin 1/2 particle [16-22]. For convenience, we quickly 
remind basics fact about the quantum field theory of a 
spin 1/2 field. Then, we turn to the calculation of p vac . 

The Lagrangian of the free model is [we treated the 
mass-less case in Eq. (35)] 



£ D irac = -* {i^d^ - m) 
where 7 M are matrices satisfying 



{^,Y} = -2.9 



/ ' ' ■-' 



(331) 



(332) 



The above equation is in fact valid for any metric tensor. 
Here, we just consider the fiat space-time case and, there- 
fore, wg will always consider that — ^t/. The quan- 
tity ^{t, x) represents a four components Dirac spinor 
and ^ = Vfrty) [this definition was already mentioned 
just after Eq. (33)]. Here, we use the Dirac representa- 
tion which means that the 7 M matrices can be taken to 
be equal to [recall that, in Eq. (31), we used the chiral 
representation] 



I 2 

-h 

I 2 ' 

h 



a 1 
-a 1 



(333) 



where <j l are the standard two-dimensional Pauli matri- 
ces. The equation of motion for the spinor is easily 
obtained from Eq. (331). Varying the Dirac Lagrangian, 
it is straightforward to show that 



which is of course nothing but the Dirac equation. Then, 
as we have done for the scalar field case, we now expand 
the ^ operator in Fourier modes according to 
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(2tt) 3 / 2 



dfe 
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E 



+ (<4) f v(k,r)e lult - lk - x 



b r k u(k,r)e 



(335) 



As appropriate for a spinor field, the operators b k and 
d k satisfy anti-commutation rules. The index r repre- 
sents the two polarization states. Inserting the expan- 
sion (335) into the Dirac equation (334) leads to two 
equations, namely 

( 7 °u; - fki - m) u{k, r) = 0, (336) 
( 7 °o;-7 l fc i +m)v(k,r) = 0. (337) 

Upon using the explicit form (333), the first equation can 
be re-written as 



(m — cj) I 



a 1 hi 



a% 2 (m + w)I 2 J " (fc ' r) 



0. 



(338) 



In order to have a non-trivial solution, the determinant 
of this matrix should vanish and this leads to 



;(fc) = \Jk 2 +■ 



(339) 



where we have used the following properties of the Pauli 
matrices, u % kia^kj = k 2 . On the other hand, Eq. (338) 
also implies that 



u(fc, r) 



v m + oj 



(340) 



where ip is an arbitrary two components spinor. Choosing 
the solution to be an eigenstate of the spin operator along 
the z-direction (in the rest frame of the particle) , the two 
components of u(k,r) can be expressed as 



t(fc, 1) = \fm + 1 



( 






(I 


a 1 




\m-\ 


- to 



i) 



(341) 



and 



u(k, 2) = \fm + i 



o % kj 
\ m + oj \ 1/ / 



(342) 



It is easy to check that this spinor is normalized according 
to 



(i^d^ - m) * = 0, 



(334) 



u{k,r)u(k, s) = 2mS s 



(343) 
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In the same manner, the solution t>(fc, r), corresponding 
to the other branch of the Fourier expansion, can be writ- 
ten as 



v(k, 1) = ym + 1 



and 



v(k,2) = \/m + > 



V 

/ a l k 



(344) 



m + uj 



J 



(345) 



The normalization of this spinor satisfies a relation al- 
most identical to that of u(k,r) except for the presence 
of a minus sign 



TT(fc, r)v(k, s) = —2m5 s 
Finally, we also have 

u(k, r)v(k, s) = 0, 



(346) 



(347) 



which can easily be obtained from the above expressions. 

All the previous considerations are standard textbooks 
calculations. We now turn to our main goal, namely the 
calculation of the vacuum energy. Therefore, the next 
step consists in calculating the energy momentum tensor. 
Either using the Noether procedure or a variation of the 
action with respect to vierbeins lead to the equation 

Tjui = --. (*7/A* + *7*fy* - 0|.*7** - fl/*7i.*) ■ 

(348) 

Of course, this expression is symmetric in fi and v as it 
should. Endowed with this formula, one can now evaluate 
the various components of the stress energy tensor in the 
vacuum state. Let us first calculate the energy density. 
One has 



(P) 



-(0 | (-2*tfib* + 2fib*ttf) | 0), 



(349) 



since 7 = —70 and (7 ) 2 = It- Inserting the Fourier 
expansion (335) into the above expression, one arrives at 

{p) = I dk f^ wt(fc ' sMfc ' s) ' (350) 

This equation can be further simplified. Upon using the 
explicit form of the spinors given above, one can show 
that 



S^v^ (k, s)v(k, s) = 4u)(k), 



(351) 



8 = 1 



and, therefore, we obtain the following expression for the 
vacuum energy 



(P) 



(2 



~ J dkcj(k). (352) 



Several comments are in order at this point. From the 
above expression, we see that we find minus four times 
the corresponding result for a scalar field, see Eq. (68). 
We have a factor four because we have two particles (the 
particle and its anti-particle) with two polarization state 
each. Of course, the most striking aspect of the above 
equation is that (p) is negative. It is therefore interesting 
to understand the origin of the minus sign in more de- 
tails. For this purpose, let us calculate the Hamiltonian 
operator associated to the Dirac field 



H 



(353) 



Inserting the Fourier expansion (335) of the spinor VP into 
the above equation, it is straightforward to arrive at 2 



H= dfe£>(fc) \(b%)H%-d%(d%y 



(357) 



However, as we have already recalled, a spinor field 
is a quantity which anti-commutes and, therefore, one 

has (d%, (dp)^\ = 5 (fc — p) 5 rs . As a consequence, the 

Hamiltonian (357) can be re-written as 



H = J dkf^ 



uj(k) [N b + N d - 6(0)} 



(358) 



where Nf, = (6|)^ b% and Nd = {d%)^ d% are the parti- 
cle and anti-particle number operators respectively. Of 
course, their mean value vanishes in the vacuum state 
which contains no particle. The zero point energy is thus 
given by the formally infinite term (5(0). Since this term 
appears with a minus sign in the Hamiltonian, the corre- 
sponding energy density is indeed negative. The origin of 
this minus sign is the ant i- commuting properties of the 
creation and annihilation operators. We conclude that 
the vacuum energy density is negative because we deal 
with fcrmions which are anti-commuting objects. 

Let us now turn to the calculation of the pressure. 



2 For convenience, let us notice that the following relationships are 
useful for this calculation 



u^(k, s)u(k, r) = 2u)5 r 
v^(k, s)v(k, r) = 2uS r 
i)t (fc, s)v(—k,r) = 0. 



(354) 
(355) 
(356) 



These three equations are responsible for the fact that only 
the quantity w(k) appears into the expression of the Hamilto- 
nian (357). 
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Upon using Eq. (348), one has 
1 



(P> 



(359) 



29^*) - ittV(^Tbflb* - <9 *7o*) (360) 



(361) 



From the above expression, we see that we have two terms 
to calculate. Using the definition of 'J', the first one reads 



dft 



(2^) 3 J v^fc) 

2 

x £V(fc,s)7V«(M)> (362) 



.s=l 

while the second one takes the form 

(o |#t 7 <ya^| o) = -(o (a^V*! o). (363) 

As a consequence, one finds that the vacuum pressure is 
given by the following expression 



M - 4 



dfe 



3 (2tt)3 J 

2 

x ^V(fc,,s)7VMM)- (364) 

s=l 

One must now evaluate the sum present in the above 
integral. We have seen that the spinors v obeys the Dirac 
equation (j°uj — j l ki + m) v(k, s) = 0. Using this last 
equation, one deduces that 

u'(fe, s)7°7 I fciw(fc, s) = uw^fc, s)u(fc, s) 

+ v t (fc,s)7 mw(A;,s). (365) 

Then, one can sum over the polarization and use the 
normalization of the spinors established before, see for 
instance Eqs. (343), (346) and (347), to obtain that 

2 

^i) t (fc,.s)7°7 l A: i w(fc,s) = Auj 2 

s=l 

2 

+ ^u t (fc,s)7°m«(fc,s). (366) 

s=l 

On the other hand, with the explicit expres- 
sions (344) and (345), it is easy to directly check that 
XL=i u ^(k, s)7°TOi>(k, s) = —Am 2 from which one can 
establish that 



t>t(fe, s)j°j l hv(k, s) = Auj 2 - Am 2 



Ak 2 



(367) 



which is exactly the sum needed to evaluate the pressure. 
Therefore, our final expression for (p) can be written as 



(P) = - 



1 



dfc- 



k' 2 



(2tt) 3 6 J uj{k) 



(368) 



We see that, again, we recover exactly the same result as 
for the scalar field but with the additional multiplicative 
factor —4, see Eq. (70). This means that all the previous 
considerations discussed in the scalar field case also apply 
to the fermion case. In particular, the regularization of 
the above divergent integrals would proceed in the same 
way. Introducing a cut-off would break Lorentz invari- 
ance and would lead to an incorrect equation of state. 
On the other hand, using, say, dimensional regulariza- 
tion produces the correct vacuum equation of state and 
leads to a regularized vacuum energy which is minus four 
times the regularized scalar field vacuum energy. We do 
not need to repeat all these calculations here and we now 
move to the case of a vector field. 



B. Proca Fields 

Having treated the case of a spin 1/2 field in the pre- 
vious sub-section, we now turn to the case of a Proca 
field [16-22]. A Proca field A^(t,x) is a massive vector, 
spin 1, field. The corresponding action can be written as 



S[A. 



1 p 



-m 2 A u A» 
2 M 



(369) 

where F^ = d ll A v — d^A^ and m is the mass of the 
particle. Again, we will restrict our considerations to 
flat space-time and, therefore, we will take for the metric 
determinant g — rj = 1. 

Before turning to the calculation of the vacuum stress 
energy momentum, we quickly remind some basics fact 
about the quantization of this type of fields. Varying 
with respect to A^, one obtains the equation of motion 



d^" - m 2 A v 



0. 



(370) 



Taking the derivative of this expressions and noting that 
dyd^F^ = (since this is the contraction of a symmetric 
expression with an anti-symmetric one), one arrives at 



M" = o, 



(371) 



provided that m =^ 0, which we assume in this section. 
As a consequence, the equation of motion of the vector 
field is given by 



d v d v A„ - m 2 A, 



0. 



(372) 



which is nothing but a Klein-Gordon equation for each 
component of A^. As usual, we Fourier expand the field 
in terms of creation and annihilation operators. This 
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leads to the following expression 



A„(t,x) = 



dk 



(2nf 2 J VMk)^ 



a iu)t-\-ik-x 



(373) 



where, as usual, one has w(fc) = \/k 2 + m 2 . Of course, 
the operators satisfy commutation rules since a vector 
field is a bosonic field. The sum over a is a sum over four 
vectors that span space-time and form a basis. At this 
point, satisfies the equation of motion but not yet 
the constraint d v A v = 0. Working out this constraint, 
it is easy to show that this amounts to e^fc M a£ = 0. If 
we now choose a frame where the wave-vector is aligned 
along the z-direction, then the basis can be taken to be 



(374) 




(375) 



where e 1 ' 2 are vectors perpendicular to the wave-vector 
(the time components of e° and e 3 carry a minus sign be- 
cause we show here the covariant vectors). But we have 
7^ despite the fact that this is the only choice 
such that e°e M3 = 0. Therefore, the only possibility left 
in order to satisfy the constraint is to take — 0. As a 
consequence, the Fourier expansion of the vector poten- 
tial is now given by 



Afi (t, x) 



dk 



a —iujt+ik-x 



(376) 



where the three vectors appearing in this decomposition 
are still given in Eqs. (374) and (375). 

Let us now come to our main subject and calculate the 
stress-energy tensor. There are several ways to carry out 
this calculation. Here, we choose to return to Eq. (369) 
and to vary this action with respect to the space-time 
metric tensor. Then, in order to work with the stress 
energy tensor in flat space-time, we will take g^ y = r}^ v 
but, of course, only after having performed the variation. 
Straightforward manipulations lead to the following ex- 
pression 



2 5S\A L 



1 W 

| \g Xa g K0 F XK F a p + X -m 2 
-g^F^F^+n^A^Av. 



(377) 



(378) 



We are now in a position where we can compute the 
various components of T^ v . Let us start with the energy 
density. From the above expression, one deduces that 

Too = ^F 0i F Qj +-^5 kl F ik F jl +^m 2 Al+^m 2 5^A i A j . 

(379) 

Then, in order to evaluate the above equation, we have 
to separately calculate the four terms that appear in this 
expression. Let us first give the relevant results needed 
to calculate the first term. Upon using Eq. (376), it is 
easy to establish that 

1 f dk 
F> m A id0 A M = —J—^ 

3 

x £e,?e> 2 , ( 38 °) 



^(OldoA^AolO) = - 



dk 



(2tt) 3 J 2w{k) 



S' J 



(381) 



a=l 



and, 



^■(Ol^AoC^olO) 



dk 



(2tt) 3 J 2uj(k) 

3 

e o e o kikj. 



S' J 



(382) 



In the same manner, the second term, containing the 
space space component of the field strength, can be ex- 
pressed as 

9 ( dk 



(2tt) 3 J 2uj{k)' 

3 



X 



Y, - e^hh) . (383) 



a=l 

The third term can be written as 

and, finally, the fourth and last term takes the form 

3 



<^(0|^|0) 



1 



/dk . 



Having calculated the four terms appearing in the ex- 
pression of T o, one must now combine them. Using the 
relationships satisfied by the vectors basis 3 , one obtains 
the following expression for the vacuum energy density 



(P) 



1 3 



(2tt) 3 2 



dk tu(k). 



(389) 



For convenience, let us give the two equations needed to com- 
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Based on our experience with the calculation of the vac- 
uum energy of a spin 1/2 field, this expression is totally 
expected. As before, it is equal to the scalar field energy 
density times a multiplicative factor which takes into ac- 
count the spin of the particle. Here, for a massive vector 
field, there are three polarization states and, as a conse- 
quence, the multiplicative factor is three. 

In order to be exhaustive, let us now evaluate the pres- 
sure. The calculation proceeds exactly as before and it 
is easy to show that 



(390) 



from which, using the relations already established be- 
fore, we immediately deduce that 



(P) 



1 



(2tt) 3 6 



dfe 



u(k) 



(391) 



Of course, this is again the expected formula, namely the 
scalar field pressure times a multiplicative factor equals 
to three. Clearly, there is no need to repeat the discussion 
about the regularization of the above divergent integrals. 
It proceeds exactly as explained before. 



C. Mass-less Vector Field 

Finally, we consider the case where the vector field is 
mass-less, m = [16-22]. This is of course necessary if we 
want to realistically evaluate the vacuum energy density 
since we know that this type of field exist in Nature. If 
the mass vanishes, this means that the dispersion relation 
is now given by 



uj{k) = k, 



(392) 



as appropriate for photons. In that case, contrary to the 
massive case, we can no longer establish that d v A v = 0. 
Therefore, the equation of motion takes the form 



duff 1 A" 



d v (d^) = 0. 



(393) 



plete the calculation. It involves simple combination of the po- 
larization vectors. Using the explicit representation given by 
Eqs. (374) and (375), one arrives at 



3 



a = l 
3 



Ecu a 
«0 «0 



fc2 
rrt 2 



(386) 



(387) 



Then, these two expressions can be combined to lead to the fol- 
lowing relation which is what is needed in order to simplify the 
expression of the vacuum energy density 



k 2 u? 



(388) 



At this stage, we need to fix a gauge. As is well-known, 
it is convenient to work in the Coulomb gauge defined by 
the following conditions 



A = 0, d.A 1 = 0. 



(394) 



In particular, the last condition implies that Ai(t,x) is a 
transverse field. As a consequence, its Fourier expansion 
can be written as 



Ai{t,x) = 



1 



(2tt) 3/2 J VMk) a 

+ \ a k) e 



dfe 



a —iut+ik-x 

a k e 



(395) 



The two three- vectors e 1 and e 2 are perpendicular to the 
wave-vectors k. In the mass-less case, only two degrees 
of freedom remain. 

The calculation of the stress-energy tensor is now 
straightforward. Using the calculation of the previous 
section, in particular the expression of the time-time 
component of the stress-energy tensor (379), one finds 
that 



Too = ^doAdoA, + -^5 kl F. lk F oh 
and, since we now have, 



(396) 



» ,J Y, <■;<:; = 2, (397) 

a=l 

2 

E (^Thkj - etefhh) - 2k 2 , (398) 

a=l 

one finally arrives at the following expression for the vac- 
uum energy density 



(p) 



1 2 



(2tt) 3 2 



dfe oj(k), 



(399) 



that is to say the scalar field energy density times two, 
as expected since there are two polarization states. For 
the pressure, using the same approach, one immediately 
obtains 



(P) 



1 2 'dfe.* 2 



(2tt) 3 6 J u(k)' 



(400) 



where, as usual, the factor two in front of the whole ex- 
pression originates from the the two polarization states 
of a mass-less vector particle. 

We conclude by stating the main result established in 
this section. For any type of fields, the vacuum energy 
density can be written as 



(P) 



1 s 



(2tt) 3 2 



dk u(k), 



(401) 



where s represents the number of polarization states 
(s = 1 for a scalar field, ,5 = 4 for a spinor field, s = 3 



39 



for a massive vector field and s = 2 for a mass-less vec- 
tor field). The pressure is given by a similar expression. 
When these divergent integrals are regularized, one can 
show that, again for any type of field, this leads to the 
Lorentz invariant form of the vacuum equation of state, 
namely (p) = — (p). In principle, the results of this sec- 
tion makes the accurate evaluation of the vacuum energy 
density in the standard model of particle physics possi- 
ble. However, before discussing this point, we need to say 
something about the fact that /? vac is positive for bosons 
and negative for fermions. 



D. Application: Super-Symmetry 

1. Motivation 

The fact that the sign of the zero-point energy den- 
sity is different for fermions and bosons, see Eqs. (68) 
and (352), suggests a very simple way of solving the cos- 
mological constant problem. It is clear that one can de- 
sign a theory where the two contributions are equal in 
absolute value such that the final result is exactly zero. 
Let us see how it works in practice. For this purpose, let 
us consider again the Lagrangian of a massive real scalar 
field, see also Eq. (6) 



(402) 



where we have introduced a constant a in front of the 
kinetic term. Before we took a = 1 and it is interesting 
to discuss why we made this choice (after all, how do we 
know that a = 1?). As usual, see Eq. (61), one Fourier 
expands the scalar field as 



1 



dfc 



(2tt) 3/2 J v / 2u(k,a) 



c k e 



— iojt-\-ik-x 



„t -.icut — ik-x 



(403) 



with a new expression for the quantity u, compare to 
Eq. (62), namely 



u>(k, a) 



k +— > 
or 



(404) 



since this is the condition for Eq. (403) to be a solution 
of the "new" equation of motion [see also Eq. (60)] 



xij) + ,u)"i),oi),o - m 2 (j) = 0. 



(405) 



The new dispersion relation can easily be understood. 
We see that the quantity a simply renormalizes the value 
of the mass m. This makes sense since, if m = 0, we 
clearly see that the relativistic relation w = k is not af- 
fected by the presence of the constant a. 



The next step is to derive the commutation relations 
for the creation and annihilation operators. Let us as- 
sume for the moment that we have the standard rela- 
tion (63) 



[cfe.cj,] =S(k-p) 



(406) 



Then, from the expression of the Lagrangian, one can 
calculate the conjugate momentum. It reads II(t, x) = 
—a$>(t,x). This implies that the commutation relation 
between the field and its conjugate momentum now reads 



[<j>(t, x),w(t, y)} = ia6{x - y), 



(407) 



an expression which reduces to the standard commuta- 
tion relation only if a = 1. Therefore, if one wants 
to have [$, II] = iS(x — y) and, at the same time, 
[cfc,cj,] = 5 (fc — p) there is no other choice that taking 
a = 1. It is worth noticing that working with the two 
previous commutation relations (406), (407) and q ^ 1 
would lead to a different prediction for the energy density 
of the vacuum. Indeed, the stress-energy tensor is now 
given by 



T, 



which implies that 



1 -8 ii d i $di<5> + -m 2 § 2 . 
2 3 2 



From this expression, it is easy to show that 
a 1 



(P) 



(2ny 2 



dfc oj(k, a), 



(408) 



(409) 



(410) 



an expression that should be compared to Eq. (68). 

Of course, one can also demonstrate that the consis- 
tency can be re-establish if the commutation relations for 
the creation and annihilation operators are taken to be 
[cfe, ct] = S (k — p) /a. In this case the commutation re- 
lation of the field and its conjugate momentum have the 
correct form. Moreover, one has 



$ 2 |0) 



1 



dfc uj 2 (k,a) 



(2tt) 3 J 2u)(k,a) 

i r dfc 



x 1 (2tt) 3 J 2w(k,a) a 



(411) 



(412) 



|$ 2 |0) 



1 



dfc 



(2tt) 3 J 2u>(k,a)a' 



(413) 



[these formulas should be compared to Eqs. (64), (65) 
and (66), which implies that 

( P ) = -^lJdku J (k,a), (414) 
that is to say the "standard" result. 
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Therefore, in conclusion, the "1/2" coefficient in the 
expression of the vacuum energy density is only obtained 
if the commutation relations and the normalization of the 
field are properly and consistently chosen. In particular, 
if one chooses to work with a coefficient "1/2" in front 
of the kinetic term (canonical normalization) , then there 
is no other choice than working with [c^, c|J = <5 (fc — p). 
Or, in other words, the previous commutation relation 
implies that a "1/2" in front of the kinetic term must be 
present. 

We now apply the previous discussion to the cases of 
a complex scalar field and to a Majorana spinor. The 
Lagrangian for a complex scalar field should be written 

as 



(415) 



with, this time, no "1/2" in front of the whole expression. 
Indeed, if one writes that 



-t= {A l + iA 2 ) , 



(416) 



then the Lagrangian takes the form of the sum of two 
Lagrangian for the two real canonically normalized fields 
Ai and A 2 . 

The previous discussion has also important implica- 
tions for the Lagrangian of a Majorana spinor. A Majo- 
rana spinor is a spinor satisfying 

K=*u, (417) 
where the charge conjugate spinor ^ c is defined by 

^J C = CW. (418) 

In this expression the symbol "T" denotes the transpose 
matrix. The matrix C is the charge conjugation operator 
that can be represented by i-f 2 -f° and satisfies C-f^C" 1 = 
— 7 M . We have of course (vE ,c ) c = ^. From a Dirac spinor 
'J one can always write two Majorana spinors 



(419) 
(420) 



Of course, the above expressions can always be inverted 
and re- written as 



(421) 
(422) 



The two previous expression can be inserted in the Dirac 
Lagrangian (331) (taking m = 0). This leads to the 
following expression 



1 r^- .-^ 



1 r*r 



£ D irac = 7= (* M1 - M M2 ) *7"0M-7S ( f Ml + 



V2 



V2 



p Ml I ")l T Ml p M2 ' A 1 M2 



-m M2 rd^ ul +i^ Mi rd^ U2 ) .(423) 



But, using partial integrations, one easily shows that 

*m 2 7 m ^ m1 = - (S M W M2 ) 7 ^ mi = +* Ml7 <^* M2 , 

(424) 

so that the last term vanishes. Therefore, one conclude 
that the Lagrangian for a Majorana spinor is given by 

^Majorana = - -* M ~ m ) *M , ( 425 ) 

and that the expansion of a Majorana spinor operator in 
Fourier can be expressed as, see also Eq. (335) 



1 



dk 



ot — ik-x 



b r k u(k,r)e- iut+ik - x 
(426) 



where we notice that only the operators b r k are present 
which reflects the fact that, for a Majorana spinor, the 
particle and the anti-particle are the same. But the most 
important aspect is the coefficient "1/2" in Eq. (425). 
Clearly, it means that if we perform the calculation of 
(p) again, we are going to find 

{p) = -T^flj dku{k) > (427) 

that is to say an expression with an overall factor 2 in- 
stead of 4 for a Dirac fermion, see Eq. (352). 

We are now in a position to address the main point 
of this section. Let us consider the following model de- 
scribed by a complex scalar field and a Marojana spinor 

C = -g» v d^*d v § - mlM - ^ M (iVfy - m F ) * M> 

(428) 

where m B is the mass of the bosonic field and m F is 
the mass of the fermionic field. From the above con- 
siderations, it is obvious to calculate the vacuum energy 
density. It reads 



(p) 



1 2 
(2tt) 3 2 

1 2 
_ (2tt) 3 2 



dk, k 2 



dk 



^k 2 + m 2 F . (429) 



We have a factor two for the complex scalar field since 
we have seen that a complex scalar field is equivalent to 
two real scalar fields and a factor —2 for the Marojana 
spinor as shown before. In general, the quantity (p) is 
not zero because m B 7^ m F . But if we now assume that 



(430) 



then the zero-point energy density exactly cancels out 
and we have solved the cosmological constant prob- 
lem. The above model is in fact nothing but the sim- 
plest super-symmetric model, the so-called Wess-Zumino 
model [19, 70-74]. 
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2. The Vacuum Energy in Super-Symmetry 

Having motivated the idea of super-symmetry, we now 
study in more detail its implications. In particular, we 
would like to check that the type of cancellation studied 
before is generic in super-symmetry. In order to sim- 
plify the discussion, we start with considering the La- 
grangian (428) with m B = m F = 0. We will now use 
the standard super-symmetric notation with the dotted 
two-components Weyl spinors [19, 74]. In this formalism, 
the Majorana spinor can be written as 




(431) 



Then, using that the Dirac matrices can be expressed as 

n ( (^)A 

7 = ^ )d a ) , (432) 
we find that the Lagrangian (425) can be re-written as 

^Majorana = (^T" M *a - ^« (a^ a& fytf*, 

(433) 

where we used that (by definition) = ($„)* and 

if" = (<J/ a )*. The next step consists in using the relation 

^K)™^ = -0 /1 * d (s* i ) ia tf Q (the spinor com- 
ponents are Grassmann variables and, therefore, anti- 
commute, see below) and, then, to transform the result- 
ing expression by integration by part. In this way, one 
shows that the Lagrangian (433) is in fact twice the first 
term. As a consequence, the Wess-Zumino model is de- 
scribed by 

Av-z = -if"dp&d v * - (if ) iQ d^ a . (434) 

Let us now consider a super-symmetric transformation. 
In order to understand what it really means, let us use 
the analogy with a gauge transformation [19, 74]. For 
definitcness, wc consider a SU(2) transformation. Let v 
be a SU(2) doublet of spinor fields. After an infinitesimal 
transformation, this doublet v becomes a new doublet v' 
such that 

„' = C ^ J ^/ 2 V ~ (l 2 + i e i?p) v, (435) 

or 

Sv = ie j ^-v, (436) 

where Uj are the Pauli matrices, the generators of the 
group. The quantity e 3 is a vector (or a collection of three 
parameters) specifying the transformation. In super- 
symmetry, the transformation is characterized by a Weyl 
spinor £ Q which, therefore, plays the role of the quantity 
e> , A super-symmetric transformation transforms a bo- 
son into a fcrmion and vice- versa. As a consequence, one 
can write 

(5$ = V2C a ^ a , (437) 
S* a = -iV2(^) a gfd^, (438) 



where the factors i or y/2 arc introduced for future con- 
venience. Before showing that the Wess-Zumino La- 
grangian is indeed invariant under a super-symmetric 
transformation, one needs to introduce a last result. We 
would like the Lagrangian to be invariant without using 
the equation of motion (in other words, we would like the 
Lagrangian to be invariant "off-shell"). In this case, one 
needs to add a term such that the new Lagrangian reads 

C w _ z = -rTdp&duQ-We, {^ f a d u ^f a -F^F, (439) 

where F is a complex scalar field. It is easy to see that 
this term is harmless since the corresponding equation 
of motion is just F = 0. But this term is necessary if 
one wants to match the bosonic and fermionic degrees 
of freedom. A Weyl spinor has two-complex components 
and wc now have two complex scalar fields. Obviously, we 
also need to modify our super-symmetric transformation 
and they now take the form 

<5$ = V2C^ a , (440) 
6* a = V2^ a F-iV2(^) a$ fd^, (441) 
SF = iV2(d^ a )(^) a0 f . (442) 

Our goal is now to show that the Lagrangian (439) is 
indeed invariant under the three transformations (440), 
(441) and (442). 

In fact, the key to show that the Lagrangian is invari- 
ant is to be able to evaluate the above transformations 
for the complex conjugate quantities. We have 

= x/2(* a )*(rr=V2*df 

= V2l,X = -V2f**- (443) 

We cannot review here the complete dotted spinors for- 
malism in detail, see Ref. [74], but in order to make the 
calculations reasonably self-consistent, it is nevertheless 
interesting to recall a few things. The first point to dis- 
cuss is the definition of the dagger symbol for Grassmann 
variables. In fact, a very clear discussion can be found 
in p. 35 of Ref. [74] and there is no need here to en- 
tirely repeat it. In brief, one needs to remember that 
we deal with quantum fields. Usually, the dagger sym- 
bols means complex conjugation and transposition. For 
a scalar field, it simply amounts to take the complex con- 
jugate of the mode function and put (or remove) a dag- 
ger to the creation and annihilation operators. For the 
spinor component $ a , this is the same (here it is impor- 
tant to remember that one discusses one component only 
and not the spinor viewed as a column and, therefore, 
we "cannot" take the transpose of the matrix -there is 
no matrix- !) and one defines as the complex con- 
jugate of the corresponding spinor component times the 
creation and annihilation operators with or without the 
dagger symbol. However, the question arises as how to 
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evaluate the Hermitian conjugate of a product of Grass- 
mann variables. Here, Ref. [74] shows that the most con- 
venient way is to define (£ a ^ a ) = ^qC"^, that is to 
say something similar to what we are used to in the case 
of two matrices even if, again, we do not deal with ma- 
trices here but with components of Weyl spinors. This 
definition justifies the first step in the previous equation. 
Then, it is also interesting to discuss how the two last 

equations are obtained. Let us write ^ = a, "Jj = b 
_i 2 

and £ = c, £ = d, where a, 6, c and d are (complex) 
Grassman variables. As a consequence, one has ^d^™ = 
ac + bd. On the other hand, £ d = e^n^, where 



c a/3 



1 

-1 



from which we deduce that £j = d and £ 2 
same manner, one has '5 = e a ^^ ^ with 



(444) 



In the 



that £ d \I> = — £ *q which is the last result used to es- 
tablish the expression (443) of 8& . In a more compact 
notation, one has 



= -1%, 



(446) 



the key point being the fact that the matrix e^g is anti- 
symmetric. 

We now need the transformation of (S^ a )* = S^a- 
one has 

ffia = V2^+iV2[(^) a p]\ d^, (447) 

= ^e/ + i^KU^^ t . (448) 

Finally, following the same lines, the transformation of 
8F* can be expressed as 



-1 

1 



(445) 



i 2 

from which one obtains ip = — b and ip = a. Using these 

results, one has £d^" = — db — ca = ac + bd = "Jd^", 
which is the wanted result. But, using again the expres- 
sions of the components obtained above, one also notices 



si* = -iV2(ty (d^ a y [(an a$ 



(449) 
(450) 



We are now in a position where one can compute the 
variation of the super-symmetric Lagrangian. Collecting 
the above results, one obtains 



6/L 



-rfd„,6&d v & - rfdn&dvS® - iS^a (^T" d^ a - (o^ 1 )"" d^ a - 5F^F - F^SF 
-iV2F% (a»f a d^ a - iV2F^ (8^ a ) K) &p f - zV2^ & {a»f a ^F 



(451) 



+iV2f {a») B6i (d^ & ) F + V5ffa^ -V2(^f a {& v ) a f) V^dpdu* - T/2rrt a d v * a d li & 
+v / 2 (a u )aa i^f a tfd^dvtf. (452) 



Each term in this expression needs to be studied in de- 
tail. The first term is also equal to iy/2F^d^ a {cr^) adt £° 
and, therefore, exactly cancels the second one. The third 
one is, after integration by parts, iy/2 (d^a) (ct) qq £ a F, 
which is also equal to -i^/2£, a (<r^) a& [d^W^j F. This 

cancels exactly the fourth term. As a consequence, all 
the term involving the field F give a total derivative con- 
tribution only. It remains the four terms involving the 
scalar field $ and & . We first focus on the terms in- 
volving $ only. For this purpose, we use the following 
relation [19, 74] 



where 



<7 = — 1(7^(7 — a tr^j = —(T . 
4 



(454) 



Therefore, the second term involving $ [i.e. the sixth 
term in the above formula (452)] can be re-written as 



= -V2 



(455) 



the second term in this equation leading to a vanishing 
(453) contribution since ~5^ v is anti-symmetric. The next ma- 
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nipulation is to integrate by part. This gives 
-V2(a»f a (a^Vctfd^ 

(456) 

Therefore, this term exactly cancels the first term involv- 
ing $ [or the fifth term in the formula (452)]. Clearly, 
this also works for the two terms involving We have 
thus shown, without using the equation of motion, that 
the variation of the Wess-Zumino Lagrangian (439) is 
a total derivative. As a consequence, the Wess-Zumino 
model is indeed super-symmetric, that is to say invariant 
under the super-symmetric transformations (440), (441) 
and (442). 

We have seen before that the constant spinor £ a was 
the equivalent of the vector e J in the case of the SU(2) 
transformation considered at the beginning of this sec- 
tion. But what is the equivalent of the generator 0j/2 
for a super-symmetric transformation? Moreover, we 
know that the generators obey the formula [cr;/2, dj/2] = 
*eiifeCfc/2. What is the corresponding equation satisfied 
by the super-symmetric generators? We now turn to 
these questions. 

Using again the analogy with the SU(2) transforma- 
tion, we expect that an infinitesimal super-symmetric 
transformation to have the form 

S^d^Qa+I^)^- (457) 

Here Q a is the super-symmetric generator, the equivalent 
of (Jj/2. Clearly, since £, a Q a must be a scalar, Q a is a 
Weyl spinor as the notation indicates. 

In order to find the super-symmetric algebra, we first 
compute the quantity (S n S^ — S^S V ) ^? using Eqs. (440), 
(441) and (442). It is straightforward to show that 

(V? - « $ = ZCVaF - 2i£ a (a^^rjh^ 

-2 v a Z a F + 2i V a K) a/ j fd^ (458) 

= -2iC W) a pTpdn* + 2ir, a (a^fd^, (459) 

since t; a r] a = r] a ^ a - But, of course, we can evaluate the 
same commutator using the expression (457) for the in- 
finitesimal super-symmetric transformation. This leads 
to 

(Ve - <W * = - (v a Q a + VaT) 

x (Vq^+^Q' 3 ) (460) 
Then, expanding this expression (carefully taking into 



account the fact that the spinor components are Grass- 
mann variables), we arrive at 

(S^s - S e 5 v ) $ = v a {Q a Qp + QpQ a )Z p <P 

+ v a (QaQ 13 + Q P Qc^I $ <s> 
- C (QaQ^ + Q^Qcr) ^$.(461) 

One must now compare Eq. (459) with Eq. (461). First 
of all, we notice that, in Eq. (459), there is no term with 
two undotted or two dotted spinors £ and 77. This implies 
that 

{Qa,Qp}= {Q & ,Q $ } =0. (462) 
Then, a comparison of the first "cross-terms" leads to 

{Q a ,Qp}* = -aK)^8^ (463) 

and a similar relation for the second cross-term. Re- 
membering that Pfj, = —idfj. (the minus sign is in fact 
just a convention; the problem with super-symmetry is 
that the literature contains many different conventions 
(A nice treatment of this question can be found in the 
first problem of Ref. [75], Appendix C), one arrives at 

{Q a ,Q $ } = 2(<jn a$ P,. (464) 

This equation shows that two successive super-symmetric 
transformations are in fact equivalent to a space-time 
translation. As we will see, this has far-reaching con- 
sequences. However, for the subject discussed here, the 
above formula has also very important implications. In- 
deed, working it out explicitly, one has 

QiQi + QiQi = 2P 3 + 2P , (465) 
Q2Q2 + Q2Q2 = -2P 3 + 2Po, (466) 

from which we deduce 

4P = QiQi + QiQi + Q2Q2 + (467) 

Therefore, since the Hamiltonian is nothing but the time 
translation generator P , this implies that 

(0 \H\ 0) = - (0 \QxQi + QiQi + Q 2 Q 2 + Q2Q2I 0) > 0. 

(468) 

We have learned two things: firstly, if the vacuum 
state is not super-symmetric, i.e. if super-symmetry is 
spontaneously broken (the underlying theory is super- 
symmetric but the solution or the state in which the sys- 
tem is placed is not), then the vacuum energy is necessar- 
ily positive. Secondly, another remarkable consequence 
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of the above results is that, when the vacuum state is 
super-symmetric, that is to say when 



VIII. THE VACUUM ENERGY DENSITY IN 
CURVED SPACE-TIME 



?o|0)=0, 



(469) 



one has (0 \H\ 0) = 0, i.e. the vacuum energy is automati- 
cally zero. We have thus proven that the cancellation dis- 
cussed at the beginning of this section is in fact not a co- 
incidence. Moreover, we have identified under which con- 
ditions it is valid: it is always true for a super-symmetric 
system. This is the reason why super-symmetry has the 
potential to solve the cosmological constant problem. It 
gives a concrete example of a symmetry which forces the 
vacuum energy to vanish. 

Unfortunately, we know that super-symmetry must be 
broken in the real world. This comes from the fact that 
each boson (fermion) of the standard model is not ob- 
served to have a super-symmetric partner of the same 
mass. This means that vacuum energy is now given by 
Eq. (429). Let us evaluate it with our "incorrect" method 
which consists in regularizing the divergent integrals with 
a cut-off M, see Eq. (75). We find 



(P) £ 



8^ 
8^ ' 



1 + 



M 2 

2 

h M 2 
- m 2 ) 



(470) 
(471) 



We see that, even if super-symmetry is broken, the most 
divergent part oc M 4 of the energy density cancel out. If 
we now properly estimate the energy density by means 
of dimensional regularization (96), one obtains that 



(P)sus 



32tt 2 



In 



4 / 4 

m I m 



32tt 2 



.(472) 



Therefore, roughly speaking, one can estimate that vac- 
uum energy is now given by ~ -^Cf USY where M SUSY is the 
typical difference between the mass of the super-partners, 
which is also the super-symmetric breaking mass. 

A last remark is in order. In this section we have 
considered super-symmetry in flat space-time while the 
cosmological constant problem is formulated in curved 
space-time. In order to be consistent, we should there- 
fore treat super-symmetry in this last context. As is 
well-known, this leads to super-gravity [19, 75]. Regard- 
ing the vacuum energy question, super-gravity presents 
important differences compared to super-symmetry. In 
particular, one loses the result that vacuum energy is al- 
ways positive. But, this also shows that discussing the 
zero-point energy density issue in curved space-time can 
bring new aspects to the problem. This is the reason 
why, in the next section, we address this question. 



In the previous sections, we have computed the vacuum 
energy density in flat space-time. This may seem inap- 
propriate since we have argued that the problem occurs 
only in curved space-time. Here we show how to justify 
the previous approach. Intuitively, this is justified be- 
cause we deal with ultra-violet divergences for which the 
large scale structure of the curved manifold should not 
play a too important role. In other words, in the ultra- 
violet regime, one only probes the local properties of 
space-time and, locally, one cannot distinguish a curved 
manifold from the Minkowski space-time. However, we 
need to put this intuitive reasoning on solid grounds. For 
this purpose, let us consider our toy model again where 
matter is simply represented by a scalar field. The only 
difference with the previous sections is that, now, this 
quantum field lives in a curved space-time. Therefore, in 
order to have meaningful Einstein equations, one must 
assume that the quantum average of the stress energy 
tensor sources the Einstein equations. In other words, 
we start from [53] 

Rimv - ^Rg^ + A B flV = K ( T ^) , (473) 

Then, by analogy with Eq. (3), we define the effective 
action W to be 



sw 



(474) 



The effective action can be written in term of the source- 
less generating functional Z[0], namely 



W 



-tlnZ[0], 



(475) 



where we recall that the explicit expression of Z[0] is 
given by, see also Eq. (97), 



Z[0] 



V<t>e 



(476) 



the quantity ^matter being the action of the scalar field. 
The above result can be proven in the following way. 
Varying Z[0], one obtains 



6Z[Q] = 

As a consequence, 
2 5W 



V(j)iSS ma Uere 
(0 |i<5S mat ter| 0) 



6Z[0] 



1 



9 Z[0] 5g^ 

2 ^S^attcr 



Z[0] 

1 

Z[0] 



-g 6gW 
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(478) 
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(0|T^|0) = (T^), (481) 
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FIG. 4: Sketch of the Riemann coordinates used in order to derive an approximate expression for the Green function, see 
Eq. (503). The plot aims at illustrating the fact that the Riemann coordinates are a local concept. Since it is sufficient to 
establish the Green function locally (and not in the entire curved manifold), one can always endow the neighborhood of a point 
x' (denoted here by the region surrounded by the closed dotted line) with special coordinates x M such that the calculation is 
simplified. It is important to notice that this does not restrict the generality of the obtained results. 



since the mean value is denned by (• • • ) = (0| • ■ ■ |0)/(0|0) . 

The next step is to calculate W explicitly in our case, 
i.e. for a scalar field. We now show that the effective ac- 
tion can be expressed in terms of the Green function. The 
Green function is defined by the following expression [53] 



-3^V M V„ + m 2 )G F (a;,2/) 



1 



-.6(x ~ y). (482) 



It is easy to show that G F (x, y) = iD F (x — y), where the 
propagator has been defined in Eq. (101). Integrating by 
part the action of the scalar field $(x), the generating 
function can be written as 



Z[J] 



X><£ cxp 



% - / d d xy^^{x) (-^"V^V* + m 2 ) $(x) + i / d d xy^y~J(x)<S>(x) 



(483) 



where we have considered the d-dimensional case for fu- 
ture convenience. This path integral is a Gaussian inte- 
gral and, therefore, can be evaluated using the standard 
techniques. For J = one has 



Z[0] = (27r)°°/ 2 (detA')~ 1/2 
= (2 7 r)°°/ 2 e5 tr ( ln ^" 1 ), 



(484) 
(485) 



where K is the following operator, appearing in 
Eq. (483), 

K(x, y) = isf—g {-gTVpVv + m 2 ) 5 d (x - y). (486) 

But the inverse of the operator K is in fact the prop- 
agator D F (x,y). A pedagogical, but simplistic, demon- 
stration of this fact is the following one. The generating 
functional can be written as 



Z[0] 



(487) 



where, in order to reproduce Eq. (486) (or, rather, to 
mimic this formula), the matrix K xy is taken to be 



K xy — K5 xy and the "number K is defined" by K = 
i\J—g (— ff^V^Vy + m 2 ). Then, the inverse of the ma- 
trix K must obey 

Y^K^K^ =8 xy , (488) 

Z 

that is to say, using the definition of K xyi 

KK-y 1 = S xy . (489) 

But this equation is nothing but Eq. (482). This shows 
that, as announced, K X y = {D p ) . As a consequence, 
one can write 



Z[0] = (2TT)°°/ 2 e^< lnD r), 



W = — tr(lnZ) p ) -iln (2vr) 



oo/2 



(490) 



(491) 



We have reduced the evaluation of the effective action to 
the evaluation of the Green function (or the propagator) 
in curved space-time. We now turn to this question. 



4G 



For this purpose, we introduce the Riemann normal 
coordinates x^ with origin at the point x', see Fig. 4 and 
Refs. [76-78]. As it will become gradually clear, the mo- 
tivation for such a choice is as follows. Since we want 
to see how the Einstein equations are modified, one only 



needs to evaluate the Green function locally. In this case, 
one can introduce coordinates that drastically simplify 
the problem, i.e. the calculation. This is exactly the role 
played by the Riemann normal coordinates. In these co- 
ordinates, the metric can be written as [76, 77] 
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Then, we insert this metric into the Green function equation (482) and retain terms with coefficients involving four 
derivatives of the metric or fewer. A lengthy calculation leads to 



Tf v dpd v G - m 



G - -R a v d v G + -R^ a v x a x^d^d v G + -R, a x a G + ( ~Ra v .,p + jR afi > v ) x a x p d v G 



+ ^R" a v p;-,^x t d lt d v G + ^R- a0 x a x^G 



r> A r-> I r>K A r> 

~30 60 a 



60 



I * a Auk 



120 



-R, 



a/3 



1 

40 



V ^d^d u R aP I x a xf ) G 



3 1 

_ TjV j TD \U 

2Q a >Pi ' 20 7 



^R K a u pR K x + ^R K a xpR K % x ) x a x^x^d v G 
60 15 
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(493) 



where all the indices are raised or lowered by the Minkowski metric and where we have defined a rescaled Green 
function according to G F (x, y) = g^ 1 ^ 4 G(x, y). The next step is to Fourier transform the above equation. In this way 
the local expansion which, in real space, looks like an expansion in the coordinates x^ is transformed, in momentum 
space, as an expansion in terms of the inverse of fc M , as appropriate for a local expansion. More precisely, defining 

1 



G(x,x') 



d d ke"> a ' 3k « x ?G(k), 



(494) 



we can expand the Green function as G(k) = Go(k) + Gi(fc) + • • • where Gi{k) is the quantity appearing in front of 
a coefficient involving the i-ih derivatives of the metric tensor [76, 77]. Let us notice that the Green function is a 
function of x' in the sense that x^ represents the coordinates of a point x defined such x' is at the origin. Another 
remark is that we define the Fourier transform of the Green function with an overall factor (2-K)~ d such that the 
coefficient of the Dirac function in the right hand side of Eq. (493) is minus one and does not contain factors of 7T, 
thanks to the formula S d (x) = l/(27r) d J d d ke lkx . 

Introducing this expansion into Eq. (493) and consistently identifying the terms of a given order leads to the desired 
result. This iterative procedure gives 
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where the coefficient a a p can be expressed as 

3 1 1 

a a /3 = -—R:a/3 - —V^d^Rap + —R a R\f} 

Inserting this last expression into Eq. (494) leads to 
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with 



fi(x,x') = f + ^R- a x a ~^a aP x a ^ (499) 
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-a x 
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The final step of the calculation consists in expressing 
the Green function in the so-called proper time formal- 
ism [76, 77]. Firstly, we replace the quantity (k 2 + m 2 )~ 1 
with 



1 



k 2 



idsc- ls ( fe2+m2 ), 
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and, secondly, wc make use of the result 
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where a = rj a ^ x a xp /2. This leads to our final expression 
for the Green function, namely 



G(x,x') 
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1 + is ft + (is) 2 h + 



(503) 



Of course, one could always push the expansion to higher 
orders in inverse of the momenta in order to obtain higher 
orders in s. But, in the present context, this will not be 
necessary. 

We now have all the tools to evaluate the effective 
action (491). We now seek an explicit expression for 
the quantity ln(D F ). As already mentioned, D F can be 
viewed as the inverse of the "continuous matrix" K which 
means that, formally, one can write 



ii<> 



ids, 
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where we have used again Eq. (501). This expression 
allows us to express the matrix element of c~ lKs . In- 
deed, given that D F = —ig~ 1 / 4 G, a direct comparison of 
Eq. (503) with the previous relation leads to 
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Now, let us consider the following integral 
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In the limit A — s- 0, up to an infinite constant that can be 
ignored, only the logarithm survives. Therefore, we can 
define the logarithm of the operator K as 



kiK = ln£> 
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Then, using Eq. (505), one obtains 
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As a consequence, we have now an explicit expression for 
the effective action, namely 



or 
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In the coincident limit x' — > x, the determinant of the 
metric tensor has been taken to one since it is calculated 
in x' , that is to say for x — 0. For the same reason, we 
also have considered that a = 0. The above expression 
can be evaluated explicitly in terms of Eulcr function and 
one obtains 
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where we have introduced the scale fi in order to main- 
tain the correct dimension of the effective action. The 
above expression can be analyzed in a dimensional regu- 
larization scheme. In this case, one is led to the following 
expression 
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which is our final result. The i factor can be absorbed 
into a constant shift so we do not need to worry about 
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it. If we compare the effective action with the original 
Lagrangian (1), we see that the first divergent term in 
the above expression can be absorbed into a redefinition 
of the cosmological constant (the second term can be 
viewed as a redefinition of the Newton constant and the 
third term leads to the appearance of new terms in the 
action, see Ref. [53]). Explicitly, one has 
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that is to say exactly the same expression as before 
(except an unimportant factor 3/2), see for instance 
Eq. (95). As a consequence in a MS scheme [16], we 
can remove the divergence and obtain exactly the same 
expression derived before, see Eq. (96), namely 
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The calculations of this section justify our previous ap- 
proach in flat space-time. Indeed, we have just proven 
that all the results derived before in flat space-time are 
in fact valid in the more rigorous approach where the 
curvature of space-time is properly taken into account. 
This confirms that the vacuum energy does not scale 
as the cut-off to the power four, see Eqs. (75) and (78) 
but is in fact proportional to the mass of the particle to 
the power four times a logarithmic factor depending on 
the renormalization scale. In some sense, the calculation 
in fiat space-time can be considered as a computational 
trick: it is sufficient to calculate the vacuum energy in 
this simple framework since, from the previous consid- 
erations, we know that this is also the result that con- 
sistently emerges from an approach where the curvature 
of space-time is properly included. Of course, physically, 
the reason for this success is not so surprising and was 
already mentioned before. Since we deal with the ul- 
traviolet behavior of the theory, a local analysis should 
necessarily lead to a correct and consistent result. 

It is also worth mentioning that Eqs. (96) and (514) 
recently attracted lot of attention, sec Refs. [47, 51]. In 
particular, as we discussed at length in Sec. IV A and as 
was emphasized in these articles, these equations are the 
consequence of the fact that the regularization scheme 
used to tame the divergence of the vacuum energy must 
respect Lorentz invariance. In fact, these equations have 
been known for a long time, see in particular Eq. (6.50) 
of the standard textbook [53]. Therefore, the fact that 
the cosmological constant is proportional to the mass of 
the particle to the power four, and not to the cut-off to 
the power four as usually claimed, see again Eqs. (75) 
and (78), is in fact not a new result. Of course, it mod- 
ifies a lot our estimate of vacuum energy. For instance, 
the correct result gives zero for the photons while, for 
the same situation, the result based on the wrong regu- 
larization scheme gives infinity. Clearly, this means that 
the correct regularization scheme leads to a cosmological 
constant much smaller than the one obtained from the 



wrong approach. As a consequence, one could hope that 
the cosmological constant problem is in fact just an arti- 
fact due to the use of an incorrect regularization method. 
We turn to this question in the next section but, unfor- 
tunately, we will see that this is not the case. 



IX. THE VALUE OF THE COSMOLOGICAL 
CONSTANT 

We are now in a position to conclude the first part 
of this review. The considerations of the previous sec- 
tions lead to the following expression for the vacuum en- 
ergy [51] 
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where p B indicates the contribution coming from the 
bared cosmological constant and where the dots means 
that other phase transition could contribute. Let us now 
evaluate the first term more precisely for the standard 
model of particles physics [51]. In this case, one has one 
scalar field, the Higgs boson with n H = 1, m H ~ 125 
GeV, the six quarks (fermions) for which n uarks = —4 
and m t ~ 171.2GeV, m b ~ 4.2GeV, m c ~" 1.27GeV, 
m s ~ 0.104GeV, m u ~ 0.24GeV and m d ~ 0.48GeV, the 
leptons (fermions) with n lopton3 = 4 and m c ~ 0.51 lMcV, 
to m ~ 105MeV, m T ~ 1.77MeV, the neutrinos (fermions) 
the mass of which is so small that we can ignore them 
and, finally, the gauge bosons (massive vector fields), 



3, TO, 



91GeV and n 
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80GeV. 



Of course, we also have the photon but since to 7 = 0, 
as already mentioned before, it does not contribute to 
the vacuum energy. The only piece missing in order to 
calculate p vac in Eq. (515) is the renormalization scale 
fi. It was argued in Ref. [51] that, because we use the 
photons coming from the supernovae to determine the 
cosmological constant and because these photons couple 
to the metric the expansion rate of which is characterized 
by the Hubble constant, one should take fi ~ E^E^^ 
with E ~ H ~ 3.7 x 10~ 41 GeV and with the en- 
ergy of the photons corresponding to the wavelength 
A ~ 500 ran. This leads to \i ~ 3 x 10~ 25 GeV and implies 
that Eq. (515) reads 
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If one does not want to estimate fi, one can simply plot 
p vac as a function of fi as done in Fig. 5. Clearly, re- 
gardless of the precise value of fi, we are very far from 
p vac ~ 10 72 GeV 4 always mentioned in the literature. 

At this stage, Eq. (516) can be seen as the "predic- 
tion" of the standard model for the vacuum energy. At 
first sight, the calculation is straightforward and, in order 
to obtain the number in Eq. (516), we have used tech- 
niques of regularization that are known to work in similar 
contexts (i.e. when one computes a cross-section in par- 
ticle physics) and to lead to a very good agreement with 
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FIG. 5: Evolution of the vacuum energy density [more precisely the first term in Eq. (515) versus the renormalization scale fi. 
In the range considered here, the vacuum energy density is negative. The "divergence" observed around log 10 n ~ 5 does not 
correspond to a new physical effect but just signal that p vac becomes positive. 



experiments. Therefore, we now turn to the question of 
what is known experimentally about the vacuum energy. 



X. MEASURING THE COSMOLOGICAL 
CONSTANT IN COSMOLOGY 

Having discussed how the vacuum energy can 
be calculated, we now address the question of 
how it is determined and/or constrained cxperimen- 
tally/observationally. Here, we should distinguish the 
case of cosmology, treated in this section, and the other 
experiments, treated in the following sections. The dif- 
ference is that, in cosmology, one claims a detection of 
A while, as we will see, the other experiments can only 
put an upper limit on its value. Moreover, in cosmology, 
there are in fact many different and independent observ- 
ables that can be used to measure the cosmological con- 
stant [23, 24, 79-82]. All these techniques are reviewed in 
other articles [40, 41] and, therefore, it is not our purpose 
here to discuss this question in details. However, for com- 
pleteness and also because this has some impact on the 
theoretical status of the cosmological constant problem, 
we now discuss how A can be determined. We will do 
so from a theoretical perspective leaving, as mentioned 
before, the technical aspects to other articles [40, 41]. 



A. The Accelerating Universe 

As is well-known, the claim that vacuum energy has 
been measured (and not only constrained) was first es- 
tablished from the discovery that the expansion of the 
universe is accelerated (although one can find many 
claims about the value of the cosmological constant in 
the history of cosmology) [23, 24]. Since the above men- 
tioned claim has clearly very important consequences for 
physics, one needs to be very precise at this stage and this 
is why, in what follows, we will review the exact origin of 
this result. 

In order to construct a cosmological model we apply 
general relativity to the Universe as a whole [83]. The 
cosmological principle implies that the Universe is, on 
large scales, homogeneous and isotropic. As a conse- 
quence, the metric is given by the following expression 

ds 2 = -dt 2 + a 2 (t)^ ) dx i dx j , (517) 

(3) 

where 7^ is the metric of the three-dimensional space- 
like sections and where a(t) is the scale factor. The 
expansion is characterized by the Hubble parameter, 
H = a/a where a dot means a derivative with respect to 
cosmic time. In the standard model of cosmology, mat- 
ter is assumed to be a collection of TV perfect fluids (at 
least in the simplest version) and, as a consequence, its 
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stress-energy tensor is given by the following expression 



N 



Tftv = ^2 T $ = (Pt + Pt + PtSV > ( 518 ) 

i=l 

where p T is the (total) energy density and p T the (to- 
tal) pressure. These two quantities are related by the 
equation of state 



Pi 



w{p T ), 



(519) 



[in general, there is an equation of state per fluid consid- 
ered i.e. pi = Wi(pi)]. The vector is the four velocity 
and satisfies the relation u^u^ = — 1. In terms of cos- 
mic time this means that = (1, 0). The fact that the 
stress-energy tensor is conserved, S7 a T a)l = 0, amounts 
to 



Pt 



3- (p T 
a 



■Pi 
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This expression is obtained when one takes p — 0. The 
case n = i does not lead to an interesting equation for 
the background. Moreover, if we assume that the fluids 
are separately conserved, then the above relation is in 
fact true for each fluid, i.e. pi + 3H(pi + Pi) = 0. 

We are now in a position to write down the Einstein 
equations (10) for the Universe. We arrive at the follow- 
ing non-linear, second order, differential equations 



k 



N 



a a k 
2- + — + — 

a a z a z 



+ ( 521 ) 

i— 1 
N 

«5>i-A eff , (522) 



From the previous formulas, as already signaled in 
Eq. (11), we see that we can also define p A by p A = 
A b /k and p A ^ by p A ^ = — p A ^ = — A b /k. In this case, 
the cosmological constant is described as fluid with a con- 
stant energy density and pressure and with an equation 
of state p A + p A = 0, i.e. 



(523) 



This is of course consistent with Eq. (5). It is worth 
noticing that the vacuum energy is described by a fluid 
with a negative effective pressure. 

Combining the two Einstein equations, one gets an 
equation which permits to calculate the acceleration of 
the scale factor 
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Let us notice that, in order to obtain the above formula, 
we have not assumed anything about the curvature k. 



This equation is especially interesting because it gives the 
condition leading to an accelerated expansion, namely 
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Since the energy density of matter must be positive, we 
see that the above condition requires a negative pres- 
sure, i.e. some exotic form of matter. But, this is of 
course exactly what happens when the vacuum energy is 
the dominant fluid in the universe. Then, Eq. (524) can 
be written as a/ a ~ 2A ett /3 which is positive for a pos- 
itive effective cosmological constant. Therefore, in the 
framework of the standard model described before, the 
discovery that a > means that a fluid with a negative 
pressure is driving the expansion. 



B. The Hubble Diagram 

In this section, we explain how the measurement of 
the Hubble diagram can be used to infer the value of the 
cosmological constant. 

Let us first consider a source in flat space-time whose 
absolute luminosity is given by L (energy emitted per 
unit of time). This quantity can be written as 



L = 



At 



(526) 



where n E is the number of photons emitted during the 
time interval At. Suppose that we observe the source 
from a distance <i L using a telescope, the surface of the 
mirror being AS*. The number of photons received is 
given by 



(527) 



On the other hand, the flux of photons (or apparent lu- 
minosity i.e. energy per unit of time and surface) through 
the telescope is, by definition, 



AS At 



from which we deduce that 
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We have just recovered the well-known fact that the ap- 
parent luminosity of a source decreases as the inverse 
squared of the distance. On the other hand, we can use 
this relation as an operational definition of the distance, 
namely d L = v /L/(4tt$ r ). 

The previous considerations can be used as a definition 
of distance in cosmology. For this purpose, let us now 
consider the same situation but in a curved space-time 
with the following metric [a form more explicit than the 
one given in Eq. (517)] 



ds 2 



-dt 2 + a 2 (t)(dr 2 + r 2 dn 2 2 ) 



(530) 
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We assume that the source is located at the co-moving 
coordinate T\ while the observer is located at the origin 
of the coordinates. In order to reproduce the above fiat 
space-time calculation in this new situation, we first need 
to calculate the surface of the sphere t = t R =const., 
r = n =const. surrounding the source. The surface is 
just given by the integral of the covariant volume element, 
namely 

S=fy/gd 2 x=[ J \Jalr\ sin(9d6>d(p = \ita\r\ . 

(531) 

This means that the number of photons received in the 
telescope can now be expressed as 



AS 



(532) 



But we have also to take into account the two fol- 
lowing phenomena. Firstly, the energy of the photons 
has changed during the propagation from the source 
to the observer. If A is the wavelength, we have A = 
(2n/k)a(t) oc v~ x , from which we deduce that E It /E E = 
v R Jv E = a E /a R , that is to say we have the redshift due 
to the expansion. Secondly, if the photons are emitted 
as a burst of duration At E , they will be received as burst 
of duration At R = At E a R /a E . We are now in a posi- 
tion where we can calculate the apparent luminosity. We 
obtain 
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where z is the redshift of the source. Now if we use the 
operational definition of the luminosity distance (529), 
we find 



d-Az) = a R ri(l +z) = a R (l + z) 
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We see from the last expression that the distance of a 
source depends on the behavior of the scale factor be- 



tween the time of emission and the time of reception. On 
the contrary, if we measure the luminosity distance ver- 
sus the redshift, we can learn about this behavior. Mea- 
suring the apparent luminosity is quite easy and, clearly, 
the difficulty lies in estimating the absolute luminosity of 
the source. It is also convenient to work in terms of the 
distance modulus, m — M, since the data are often pre- 
sented in this way. The quantity m (or sometimes m\, \) 
is the apparent bolomctric magnitude. It is related to 
the apparent bolometric luminosity I by 



I = l(T 2m/5 x 2.52 x 10" 5 erg/cm 2 
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For instance the sun has an apparent bolometic magni- 
tude of —26.85. On the other hand, the absolute bolo- 
metric magnitude M is the apparent bolomctric magni- 
tude that the source would have at a distance of 10 pc 
and is given by 

L = 1(T 2M/5 x 3.02 x 10 35 erg/sec . (536) 
For instance, the absolute bolometric magnitude of the 
sun is 4.72. As a consequence, using these two definitions, 
the distance modulus is defined by 

d L (z) = 4r< ( z ) = 10 1+(m ~ M)/5 x l(T 6 Mpc , (537) 

such that d L is dimensionless and c/Hq = 3000 /i _1 Mpc 
is of course nothing but the Hubble length (we have re- 
established the speed of light for convenience) . From the 
above expression, one immediately deduces that 



H(z) = m - M = 5 log 10 d L (z) + 25 + 5 log 10 I — J , 

(538) 

where c is expressed in kmx s _1 , Hq in kmxs -1 x Mpc -1 
and c? L in Mpc. 

In practice we only have access to d L over a limited 
range of rcdshifts. Therefore, one can derive a general ex- 
pression by Taylor expanding the scale factor. For small 
redshifts, we have 



i(t) = a(t R ) + a(t R )(t -t R ) + -a(t R )(t - t R f + ■■■ = a(t R ) 
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where q = — ad/(a) 2 is the acceleration parameter. A positive acceleration parameter corresponds to a decelerating 
universe whether a negative one corresponds to an accelerating universe. Using the previous expression, the luminosity 
distance becomes 
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As expected the result only depends on the time of flight t E — t n . In order to obtain our final result, we must 
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FIG. 6: Luminosity distance d L (z) (left panel) and distance modulus p(z) (right panel) in a matter dominated Universe. The 
solid line corresponds to the exact result. The dashed line is the linear law d L ~ cz/Ho and the dotted dashed line represents 
the approximated result at second order, see Eq. (544). We have used h — 0.72. 



express this time of flight in terms of the redshift. The 
redshift of the source is denned by z = a(t R )/a(t E ) — 1. 
Using the Taylor expansion of the scale factor, we end up 
with 

z = -H R (t E -t R ) + H 2 R (l + ±q u yt E -t R ) 2 + - ■ ■ , (542) 

which can be easily inverted to give the expression of the 
time of flight, namely 



z 1 



t E -t R =-- + -{l + - qR )z< 
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The last step is to replace the above expression into the 
Eq. (541). This gives 



(544) 



where we have taken into account that the time of re- 
ception is the present time (denoted with the subscript 
"0"). We see that the measurement of the luminosity 
distance versus the redshift allows us to determine the 
Hubble parameter today. If one goes sufficiently far in 
rcdshifts, one can determine whether the universe is ac- 
celerating or not via the accelerating parameter go- It is 
worth noticing that this last conclusion is independent of 
the underlying gravity theory (provided this is a metric 
theory). The typical behavior of the luminosity distance 
and of the modulus is illustrated in Fig. 6 for a matter 
dominated universe. Let us also notice that the method 
that we have just presented for pedagogical reasons is of 
course not used in practice because we have now mea- 
sured up to redshifts of 2 ~ 1.4 for which the Taylor 
expansion is not valid (and, of course, in the real world 
the analysis is always more complicated that the simpli- 
fied version discussed above). In this case, one must use 



an exact expression, that is to say one must compute the 
integral in Eq. (534) exactly. 

If we now use general relativity (i.e. a specific theory 
for gravity), one can go further and relates the value of 
the acceleration parameter to the matter content of the 
Universe. The critical energy density p cri , defined by 
Peri = 3H 2 1 k, is used to evaluate the contribution of 
each species through the parameter Qi = pi/p cl -i where 
Pi is the energy density of the species i at present time. 
Then, the Friedman equation can be re-written as 
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and the acceleration parameter can be re-expressed ac- 
cording to 



1 N 

= - &i (1 + 3u>,-) 



(546) 
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The universe contains pressure-less gas (cold dark matter 
and ordinary baryonic matter) for which w = but also 
radiation (w = 1/3) and, finally, vacuum energy (w = 
— 1). In this case, one obtains 



q = -fi m + £! 7 - fi v 



(547) 



Today the contribution of radiation is negligible since we 
know from the measurement of the CMB that il^h 2 ~ 
2.48 x 10 -5 [81]. Therefore, the acceleration parameter is 
in fact given by qo ~ fi m /2 — ^ vac - On the other, we also 
know that the universe is spatially flat which means that 
O m + f2 vac ~ 1 [81]. As a consequence, from the mea- 
surement of qo — —0.67 ± 0.25 [23, 24], one can deduce 
the two quantities f2 m and r2 vac . On finds that f2 m ~ 0.3 
and Q ~ 0.7. This is how one reaches the conclusion 
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that our universe is accelerating and that vacuum energy 
dominates the present energy content. Remarkably, as 
announced at the beginning of this section, it is a detec- 
tion and a measurement of p vac . One finds 
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10 
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GeV 4 , 



(548) 



a number that should be compared with Eq. (516). 

At this stage, several remarks are in order. Firstly, we 
notice the large mismatch between the theoretical expec- 
tation (516) and the above, observationally determined, 
number (548), something like 54 orders of magnitude 
(but much less than the 122 orders of magnitude often 
quoted in the literature). This is of course nothing but 
the cosmological constant problem although we will be 
more accurate with regards to its definition in the sub- 
sequent sections. Something is clearly wrong but, as will 
be discussed, it is difficult to identify where the mistake 
is. Secondly, it is clear that the above conclusion rests 
on many assumptions. For instance, we have assumed 
that the universe is homogeneous and isotropic [42] and 
that gravity is well described by general relativity. This 
means that, a priori, the cosmological constant problem 
is present in this context only. Thirdly, we have also im- 
plicitly assumed that the reason for the accelerated ex- 
pansion is vacuum energy. This is clearly a natural idea 
since there is a term in our theory that must be present 
since it satisfies all the required properties (covariancc, 
energy conservation) and that precisely leads to the ob- 
served phenomenon. A priori, nothing more can be asked 
to a new theoretical framework. But, given the mismatch 
discussed before, one can also doubt the identification of 
the source of the accelerated expansion with vacuum en- 
ergy. As a matter of fact, it is possible to construct mod- 
els where the universe accelerates because of some new 
source of matter (often named "dark energy" ) . This is 
for instance the case of quintessence and/or galileons [43] 
models where a scalar field is responsible for the acceler- 
ation. In this case, one can no longer claim that we have 
detected the vacuum energy and that its value is given 
by Eq. (548). However, this does not solve the cosmo- 
logical constant problem because, even if the source of 
the acceleration is another fluid, we still have the con- 
straint p w < p CTl , Clearly, the number (516) does not 
satisfy this inequality which means that the problem is 
still present. In this sense, the constraint deduced from 
cosmology is very important: even if dark energy is not 
the cosmological constant, the observation that the en- 
ergy density today is the critical energy density severely 
limits the value of A off . Fourthly, in principle, it is pos- 
sible to check experimentally whether the reason for the 
accelerated expansion is vacuum energy or something else 
like quintessence [84-99] (more recently, in Ref. [100], 
the current acceleration was also explained by a scalar 
field experiencing quantum fluctuations during inflation) . 
Measuring the equation of state parameter w is the clue 
since, if this is a constant equals to —1, then one can 
be sure that the source of the acceleration is the cos- 
mological constant. Unfortunately, this measurement is 



difficult [101]. 

One more general grounds, it is clear that it would 
be better to have constraints, or even a measurement, of 
vacuum energy in an experimental context which is is not 
cosmology where, as reminded before, measurements are 
difficult and always subject to many biases and system- 
atic errors. In the next sections, we explore whether this 
is possible. 



XI. MEASURING THE COSMOLOGICAL 
CONSTANT ELSEWHERE THAN IN 
COSMOLOGY? 

In the previous section, we have written the Ein- 
stein equation in an homogeneous and isotropic universe. 
Here, we investigate the influence of the cosmological 
term in a spherically symmetric situation. This will allow 
us to study how the presence of A eff affects the motion of 
the planets in our solar system and how it modifies the 
energy levels of the Hydrogen atom. 



A. The Static and Spherically Symmetric 
Gravitational Field in Presence of a Cosmological 
Constant 



Let us analyze the influence of the cosmological 
constant in a static and spherically symmetric situa- 
tion [102]. This means that, in spherical coordinates, 
one can write the metric tensor as [102] 

ds 2 = -B(r)dt 2 +A(r)dr 2 +r 2 (ck9 2 + sin 2 6>d^ 2 ) , (549) 

where A(r) and B{r) are two free functions. We seek 
solutions of the Einstein equations (10) with a vanish- 
ing stress energy tensor. It is easy to show that these 
equations can be expressed as 
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The next step is to calculate the components of the Ricci 
tensor for the metric tensor given by Eq. (549). Straight- 
forward calculations lead to the following expressions 
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the other components being zero. In the above equations, 
a prime denotes a derivative with respect to the radial 
coordinate r. 

We are now in a position where one can determine the 
solution of the field equations. If we take the combination 
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Rtt/B + Rrr/A = —A eS B/B + A e? A/A 
Einstein equations (550)], one obtains 
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[from the mass m given by 
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which implies that B(r) = C/A(r), where C is a con- 
stant. Using this expression, one can re-express Rgg as 
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and, as a consequence, the corresponding Einstein equa- 
tion now reads 



rB'(r) +B(r)=C- CA a!! r 2 



(557) 



This equation can easily be integrated and one obtains 



B(r) = D + j-DA^j, 



(558) 



where D is another integration constant. The constants 
C and D are fixed in such a way that, in absence 
of a cosmological constant, one recovers the standard 
Schwarschild solution. This amounts to take D = 1 and 
C = 2GM, where M is the mass of the central body and 
G the Newton constant. This completes our calculation 
since the two functions A(r) and B(r) are now completely 
specified, namely 



A(r) 



1 



B(r) 



B{r) = 1 - 



2GM A 



(559) 



We can now compare this solution to the standard 
Schwarschild solution. The fact that the time-time com- 
ponent of the metric tensor is inversely proportional to 
the r-r component is still true but the radial dependence 
of the function B(r) is modified by the presence of a 
cosmological constant. In the Newtonian limit, we have 
gtt — 1 + 2V(r), where V(r) is the Newtonian potential. 
In the present context, this implies that 



V(r) 



GM A 



6 



(560) 



Therefore, we obtain the usual gravitational potential 
corrected by a term proportional to the cosmological con- 
stant. This term looks like an inverted harmonic oscil- 
lator. This corresponds to a force felt by a test body of 



GMt 



-u 



mA „r 

cif 



u. 



(561) 



where it is a vector directed towards the central body 
of mass M. The cosmological constant force is directed 
in the other direction and can therefore be viewed as a 
kind of "anti-gravity force" . We also notice that it is pro- 
portional to r which means that its effect will be more 
important on large scales. This last property is the clue 
to understand why measuring the cosmological constant 
is more efficient on large scales (the best example being 
of course cosmology). This means that the modification 
of the planet orbits and/or the energy levels of the atoms 
will certainly be a small effect. On the other hand, we 
know that measurements in, say, atomic physics are ex- 
tremely accurate and, a priori, there is the hope that 
this accuracy could compensate the smallness of the ef- 
fect. Unfortunately, as we are now going to study, this 
will not be the case. 



B. Planets Orbits 



Our goal in this section is to study how the planet or- 
bits are modified by the presence of vacuum energy [103- 
105]. We will follow the treatment of Ref. [103]. We as- 
sume that the gravitational field created by the sun and 
the cosmological constant is given by Eqs. (549) and (559) 
with M = Mq. The motion of a planet is represented 
by the motion of a test body in this background gravita- 
tional field. The corresponding geodesic equation can be 
written as 



dr 2 



,, dx u dx x 



dr dr 



(562) 



where r is an affinc parameter (not necessarily the proper 
time). Then, after having evaluated the Christoffel sym- 
bols for the metric (549), one can express these equations 
explicitly. One obtains 



■55 




dr 2 



2A 



dr 



2A 



r fd6 
A \dr~ 

2 dr d(9 



dr 2 4 

• 2 

r sm 



.4 



S' dt dr 
~B drdr 

dy> 
d7 



d^? 

dr 2 r dr dr 

d 2 u? 2 dr dtp 

— - H + 2 

dr 2 r dr dr 



— cos a sm 



dip"' 
dr, 
cos d# d<£ 
sin dr dr 



(563) 
(564) 

(565) 
(566) 



Let us now solve these equations. Firstly, we can assume 
that the orbit stays in the 9 = tt/2 plane. This auto- 
matically solves Eq. (565). Then, if we divide Eq. (563) 
by dt/dr and Eq. (566) by dip/dr, we obtain two new 
equations, namely 
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= 0. 



(567) 
(568) 



The first equation can be used to define the affine pa- 
rameter. We choose it such that it satisfies dt/dr = 1/B. 
Then, the second equation tells us that 



2 dl P T 

r d7 = J ' 



(569) 



where J is a constant. Finally, Eq. (564) can be re- 
written as 

2 

= 0. (570) 



dr 2 



2A 



J 2 

Ar 3 



2AB 2 



If we multiply this last equation by 2^4dr/dr, one arrives 
at 



_d_ 



.4 



dry 
dW 



1 

B 



0. 



(571) 



This means that the quantity inside the square brackets 
is a constant. Let us call this constant — E. Then, this 
gives an equation expressing how the radial coordinate r 
varies with the angle ip, 



dr 



A(r) 



1 



1 



J 2 B{r) 



E 
J 2 



(572) 



It is easy to express the constants E and J in terms of 
the aphelion and perihelion r + and r_ , which are the 
two points on the orbit such that dr/df = 0, see Fig. 7 
(where the elements of the orbits are represented) . This 
gives 



1 

J 2 

E = 



B- 



B 



-l ' 



JB+ - r 2 _/B. 



(573) 
(574) 



with B± defined by B± = B(r±). 

So far, we have been very general and we have never 
specified the function B(r). From now on, we are going to 
use the explicit form of this function, see Eq. (559). Let 
us notice in particular that, in absence of a cosmological 
constant and in the limit where GM/r <C 1, the two 
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above expressions reduce to 
E - 1 GM 



, J 2 ~ 2GMD, (575) 



where D = r+r-/(r++r-) (of course, not to be confused 
with the integration constant introduced in the previous 
sub-section). As is shown in the following footnote, these 
are exactly the relations obtained in the Newtonian case. 
More generally, Eq. (572) takes the form 



dr^ _ 1-E 4 2GME 3 2 



dip) 



+ - 



J 2 J 2 
A „r 4 EA „.r 6 



3J 2 



2GMr 

(576) 



Performing the standard change of variable, u = 1/r and 
differentiating with respect to ip, we obtain 4 



This is the differential equation that needs to be solved 
in order to find the trajectory of the planet. It should be 
compared with its Newtonian counterpart, see Eq. (583). 
One can check that, in this regime (without a cosmologi- 
cal constant), the two equations are indeed identical since 
according to Eq. (575), E — s- 1, A cff = and the term 
3GMu 2 can be neglected. The rclativistic term SGMu 2 
is responsible for the precession of the perihelion and the 
last term represents the correction due to the cosmolog- 
ical constant. 

We now turn to the solution of Eq. (585). Let us first 
review how one calculates the precession of the perihelion 
without a cosmological constant [102]. The idea is to 
restart from Eq. (576). For A off = 0, the polynomial on 
the right hand side can be written as 



E 



J 2 



r(r — r + )(r — r_)(r — e) 



(586) 



since we know it has to vanish at r = r±. The only 
quantity that remains to be determined in this expression 
is e. Comparing the above relation with Eq. (576), we 
deduce that 



4 It is interesting to compare this result with the Newtonian anal- 
ysis. In this case, the Lagrangian is given by 



m , ., .,. GMm 
L=-(r 2 +r 2 V 2 ) + , 



(577) 



where a dot denotes a derivative with respect to time. Again we 
have chosen the orbit to be in the plane 9 = tt/2. This leads to 
the following equations of motion 



mr — mr<p + 



GMm 



0. 



mr (p + 2mrrip = 0. 



(578) 
(579) 



The second equation gives ip = J/r 2 while the first one can be 
re-written as 



d /r 2 J 2 GM\ 
di V ~2 + 2T 2 ~r / 



(580) 



The quantity between the parenthesis is a constant and we call 
it —E. Then, this equation can be re-expressed as 



/dr 



2E , 2GM , , 

r H r s — r . 

J 2 J 2 



(581) 



From this equation one can deduce an expression for E and J in 
terms of the aphelion and the perihelion. One gets 



GM 
r+ + T- 



J 1 = 2GMD. 



(582) 



Introducing the variable u = 1/r and differentiating once more 
with respect to ip, one arrives at 

d 2 u GM , , 

^ +U =1P- (583) 

This equation can easily be integrated. The solution reads 

, GM GM , , . 

u(ip) = Fcos [ip - ipQ) + -j^- = -j^- (1 - e cos (/3 ) , (584) 

where F and ipo are tw° arbitrary constants. Clearly, one can 
always take ipo = 0. The above trajectory represents, as is well- 
known, an ellipse with eccentricity e. 



1 - E 

- 1 



-e = 2GM. 



(587) 



Using the expressions derived before, it is easy to show 
that e = 2GM/ (1 — 2GM/D). As a consequence, one 
can write 



J 



<P+-<P- = 



1 



VT=EJr_ y/(r-r+)(r-r-) 
1 dr 



^Y^Ifr r 



(588) 



There are several methods to evaluate this integral. As a 
matter of fact, it can be integrated exactly in terms of El- 
liptic functions [48, 49]. But, being given that e is a small 
quantity, it is better to expand the result in this parame- 
ter. This leads to easier calculations that are, moreover, 
more explicit. In addition, the method can be general- 
ized to the case of a non- vanishing cosmological constant, 
which is of course our main goal. In fact, if we restart 

from Eq. (572) and use that B^ 1 ~ 1 + 2GAf/r H , 

we see that the term in the squared bracket in Eq. (572) 
is a quadratic function of 1/r which vanishes at r = r±. 
This means that one necessarily has 

1 *-*«c(i-l-)tt-±).<m 



J 2 B(r) J 2 r 2 



r r_ / \ r r. 



Then, writing Eq. (572) in terms of u rather than r, one 
arrives at 



U \ ~^2^(u-u_)(u-u+). (590) 



d<^ / A(r) 



Now, when A otf ^ 0, as is clear from Eq. (576), the 
term in the square bracket of Eq. (572) is no longer a 
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quadratic function of 1/r even if the expansion of the 
function i? _1 (r) is used. However, this term still vanishes 
at r = r± and the integral will be dominated by the con- 
tributions coming from the vicinity of these points, see 
also Eq. (588) . The idea is then to continue to work with 
Eq. (589) and to take into account the cosmological con- 
stant in the definition of C [103]. This should provide a 
reasonable approximation of ip + — </?_. Using Eqs. (573) 
and (574), one can write 
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J 2 B(r) 
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7 



(BZ 1 - BZ 1 ] 



.2 



c 



1 



1 

r+ 



(591) 



Then, if we differentiate both sides twice with respect to 
u, this gives 



C 



(u+ — u~) (u+ — u-) A" (it) 



-1 



2(u + -U-)A'(u) 
uA" 



1. 



.4' 



(592) 
(593) 



u—L 



where L = 2D is the semilatus rectum. In the above ex- 
pression, we have used A + — A_ ~ (u + — U-) A' (L^ 1 ). 
So far, the method is completely general provided that 
A -1 = B. We now use the explicit form of A^ 1 . Expand- 
ing everything in GMu and in A off , one obtains [103] 



C~ -1 



4GM A„L 3 4A„L 3 



GM 



(594) 



The final step consists in using the above expression in 
Eq. (590). One finds [103] 



A l / 2 {u)du 



u- [C (u — u+) (u — u_)] 1//2 
A 2G.M . A .. //' 
2~~L 



7T 1 



2GM 



,(595) 



where we have expanded A 1 / 2 C x l 2 and where we have 
used the fact that J™ + du/[(u— u + )(u— it-)] 1 / 2 = ln(u + — 
u_) — ln(u + — u_) + 2i7r/2, the complex i being exactly 
cancel by the i coming from the fact that C, which ap- 
pears in a square root, is negative [its expansion starts 
with —1, sec Eq. (594)]. As a consequence, we obtain 



Aip = 2 (ip + - ipJ) - 2tt 



67r — j- + Ai Pa 



(596) 



where the first term is the standard rclativistic perihelion 
precession while [103] 



GM : 



(597) 



is the contribution due to the cosmological constant. The 
above equation is the main result of this sub-section. It 
gives the perihelion precession due to the presence of vac- 
uum energy. Let us notice that the scaling Atp A oc L 3 in- 
dicates that the precession is larger for planets that are 
far from the Sun. It is of course totally consistent with 
the fact that the effect of the cosmological constant is 
more important on large scales. 

We now use the above result to constrain the amount of 
vacuum energy present in our universe. For this purpose, 
let us now calculate this perihelion shift for, say, Mercury 
(as mentioned before, the effect is larger for other planets 
but this does not modify the conclusions obtained below; 
estimates for Venus, Earth and Mars can be found in 
Ref.[105]). In this case, we have r+ ~ 68.8 x 10 6 km and 
r_ ~ 46 x 10 6 km. This implies that D ~ 27.7 x 10 6 km. 
We also have GM ~ 2.95km. Finally, it is interest- 
ing to calculate in terms of the vacuum energy den- 
sity observed in cosmology, see Eq. (548). For the cos- 
mological constant, this leads to the following number 
A obs ~ 8.81 h 2 x 10~ 48 knr 2 . This leads to the following 
expression 



Aip A ~ 1.6h 2 x 10 



-24 



A 

aZ 



rad/revolution, (598) 



or, since Mercury completes 415 revolutions each century 
and there are 360 x 60 x 60/(27r) arc-seconds per radian, 



2.7h 2 x 10 



-16 



A 

aI 



arc-second/century, 



(599) 

which is, of course, completely unobscrvable. The reason 
for such a result has already been mentioned. The solar 
system has a characteristic scale (say a few astronomical 
units) which is too small for the cosmological constant to 
leave a sizable imprint. However, even if one cannot mea- 
sure the vacuum energy density with this method, it can 
nevertheless be used to constrain its value. Indeed, the 
uncertainty in the precession of the perihelion of Mercury 
is about 0.1" par century, which means that 



Aip A ~ 2.7h 2 x 10~ 16 [j L -)< 0.1, 



implying that 



Pvac £ 10 15 p crl - 3 x 10 _aa GeV 



32 , 



(600) 



(601) 



Of course, this limit is not competitive with what one 
has obtained in cosmology. However, when compared to 
Eq. (516), it leads to an interesting piece of information 
which allows us to define what the cosmological constant 
problem is more accurately. Indeed, we see that there is a 
contradiction between the theoretical expectation (516) 
and the observations of the planet trajectories in our so- 
lar system. In other words, the predictions of Eq. (516), if 
true, would lead to orbits drastically different from what 
we observe. Therefore, even without the cosmological 
observations and only from the motion of the planets in 



■58 



the solar system, one sees that our calculation of vacuum 
energy must be incorrect. Put it differently, it would 
be incorrect to say that the cosmological constant prob- 
lem originates only from the recent astrophysical obser- 
vations. In some sense, the problem is much worse. Only 
by observing the planets in the solar system, a class of 
observations which seems to be much more straightfor- 
ward than the cosmological observations, we know that 
our regularization method of the vacuum energy must be 
flawed. 

C. Modifications of the Atomic Levels 

We have just seen that observing the motion of plan- 
ets in our solar system can tell us something about the 
value of the vacuum energy. A priori, another possibil- 
ity to observe the influence of the cosmological constant 
is through the modification of the atomic energy levels 
that its presence would cause. In this section, we inves- 
tigate this question in the case of the Hydrogen atom 
since this is technically simple [106— 1 1 1] . We will follow 
the treatment of Ref. [106]. In particular, since the prob- 
lem is still spherically symmetric, one can use the metric 
tensor (549) obtained before. 

The motion of an electron around the nucleus, living in 
a static and spherically symmetric manifold, is described 
by the following Dirac equation 

(T^ + to) * = 0, (602) 

where VE" is a Dirac spinor, see Sec. VII A and to the mass 
of the electron. In the above equation, are the Dirac 
matrices satisfying the relation 

{T",T u } = 2g' iV . (603) 

Let us notice that we have already presented this formula, 
see Eq. (332). Here, we slightly change the notation and 
now denote the Dirac matrices in curved space-time by 
(notice that we have also included a factor i in the def- 
inition of r M ). The matrices r M are related to the Dirac 
matrices in flat space-time, 7 M , by = e a ^ a , where 
e a ^ are the vierbein fields such that g^ v = e a ^ v T) a p. 
The symbol Z? M denotes the covariant derivative acting 
on the spinor, namely 

= (0„-av)tt, (604) 
where the spinorial connection can be expressed as 

k-V = -\l a lpe a y g vX (V M e" A ) + iqA„, (605) 

with V M e^A = d^e 13 \ — T^^er a an d A^ is the vector field 
describing the electromagnetic field acting on the electron 
(q is the charge). The Dirac equation can also be re- 
written as 

ido^ = HV, (606) 
where the effective Hamiltonian is given by 

H = -i (.g 00 ) _1 r°r (di* - uji$) + iT y 

-i ( 5 00 )" 1 r°TO*. (607) 



In this case, the Dirac equation looks like an effective 
Schrodinger equation. 

It is clear that the typical dimension of an atom is very 
small in comparison to the typical curvature scale of the 
manifold. Therefore, as was done for the calculation of 
the Green function in Sec. VIII, one can perform a local 
analysis. For this purpose, we now use the Fermi normal 
coordinates [78]. In these coordinates, the metric tensor 
takes the form 

.goo = -f - R omm x e x m + • • • , (608) 
So* = —Rotimx'x™ + ••• , (609) 

g%j = kj - ^Riej m x £ x m H , (610) 

where the dots indicate higher order terms. This implies 
that the veirbein fields can be expressed as 

e o = o o - zomx x , (611) 
e i = d i — -R HraX x , (612) 
and the Christoffel symbols can be written as 

^ij = g (Roijm + Rojim) X™ , (613) 

rg< = RoiomX™, (614) 

Fjfc = — (Rjikm + Rkijm ) X™ , Tqq = 0, (615) 
r 0j - = Ro m jiX m , r 00 = Ro i0m x m . (616) 

Using the above results, this finally leads to the following 
expression for the spinorial connection [106] 

= \poijR J oo m x m + -■ykljR kj omX m 

+iqA , (617) 

w» = \loljR° 3 im x m + ^k^R k3 im x m 

+iqAi. (618) 

The above formulas will allow us to calculate the effec- 
tive Hamiltonian (607) explicitly. But before reaching 
this stage, we must also determine vector field A^. The 
Lagrangian of the corresponding gauge field is given by, 
see also Eq. (369) 

S = - J d 4 x^ {^g af) g^F^F u[3 - fA^ , (619) 

from which we deduce that the equations of motion read 

g»- v V ll V v A 1 -R%A (T = -j 1 . (620) 

At first order in the curvature, these equations can be 
expressed as [106] 
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S^didjAo + ^R ejm x e x m d l d 3 A + ^R mm x m d l A Q + 2R k Mm x m d l A k - ^R im x m d l A Q = -j 
^ 9i9j A k + + \^ + ^A a - 2 -R* k A„ - lR» 00k A p 



RioomX m d l Ak — —R 



2 

— llprri >* 



l d p A k 



3k- 



(621) 



(622) 



For a one electron atom with a nucleus of charge Ze, one 
can take 



first order perturbation SA k obeys the equation [106] 



j° = -ZeS(r), j k = 0. 



(623) FfkdiSM 



In what follows, we determine the vector perturbatively 
in the curvature. At zeroth order, in flat space-time, the 
equation for the time component reduces to 



S^didjAo = ZeS(r), 
whose solution is of course 
,(o) 



Ze 1 
47r r 



(624) 



(625) 



Then, we write Aq = A^ + <L4o, where 5Aq is first order 
in the curvature. By inserting this expression into the full 
equation (621) and using the zeroth order solution (625), 
one arrives at [106] 

PdidiSAo + (3R°mm ~ Rim) x e x m = 0. (626) 

The solution of the above equation reads [106] 

5A = ±-^(R + 4R 00 ) r 



1 Ze , o 
~12 4^ 3i? mm ~ Re 



(627) 



In the same manner, one can determine the space compo- 
nent. Obviously, at zeroth order, one has A[ 0) =0. The 



2Ze 1 2Ze 



whose solution can be written as [106] 



1 Ze 1 Ze n x e x m 
5A k = o~rRokr + - — R\ km . 

2 An 6 47r r 



0, (628) 



(629) 



This completes our calculation of the electromagnetic 
field. 

Having determined the vicrbcins fields (i.e. the quanti- 
ties describing the gravitational force) and the gauge field 
(i.e. the quantity describing the electromagnetic force), 
we are now in a position where we can calculate the 
Hamiltonian (607). In order to make contact with the 
usual treatment of the Hydrogen atom, it is more conve- 
nient to introduce the standard Dirac matrices a, and (3 
related to the flat space-time Dirac matrices by 



7o = if3, 7, = -if3cti. 
As is well-known, they satisfy 

OLiOLj + OLjOii = 2Sij, 

otifi + /3a t = 0, P 2 = 1. 



(630) 



(631) 
(632) 



Then, we use the expressions of the spinorial connec- 
tions (617), of the gauge field (625), (629) and of the 
Dirac matrices in the equation (607) of the Hamiltonian. 
This leads to the following expression 



H = -ia l di + mp - - - ^R omm x i x m a i d i - % -R ie , jm x l x m a> d l - l -R mm x l x m a> a 1 d t - ^-R oeim x e x m d l 
r 2 6 6 2 

+^a i oPa k R jkim x m + l -a l a 1 R 0]ml x m + a 3 Ro mlJ x rn - l -a 3 R ojQm x m + ^(R oeml x e x m y + ^Ro^r 

+ 777C (R + 4i?oo) r - — C (Ri m + ZRoiOm) -R mm x e x m pa l + -Ro i0m x e x m p. (633) 

12 12 r b I 

I 



In this expression we have defined ( = Ze 2 /(47r) and correction due to space-time curvature. Since the above 
we have taken q = — e. The three first terms are the expression is quite complicated, it is interesting to es- 
usual Dirac Hamiltonian. The other terms represent the timate the order of magnitude of each term. The cur- 
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vature is given by D~ 2 where D is a typical distance. 
Each x 1 is typically of the order of the Bohr radius, 
i.e. x l ~ ao ~ m~ 1 e~ 2 ~ m -1 £ -1 - As a consequence, 
<9i is of order to£. This also means that the typical veloc- 
ity is given by mQ/m ~ £. Of course the matrices a 1 and 
/? are dimensionless but, in the non relativistic limit a l di 
must reduce to ~ mv 2 which implies that a 1 ~ v ~ £. In 
the following wc consider that £ = Za where a is the fine 
structure constant is a small number. Then, the domi- 
nant term in the above expression is the last one, which 
means that the Hamiltonian can be approximated as 



H 



Hq + H. nt ~ a ■ p + m/3 
1 



-mR eom^x m /3, 



(634) 



where we have used pi = —idi. If we use the following 
representation for the Dirac matrices 



er' 
a 1 



h 

o -h 



(635) 



and write the Dirac spinor as 



* = e 



iEx° ( 4> 



(636) 



then the Dirac equation reduces to two equations that 
can be expressed as 



a px+ [ m h \mR omm x l x m 

r 2 



71 + r + \ mR oeomX e x m ) 



From the second relation we can express x m terms of 
a ■ p<p and then we can insert this expression into the 
first equation. This leads to 



C 1 

E + m-\ h -mRoiomX x r> 

r 2 



71 h ^mR omm x e x m 



(639) 



or, using a ■ p 4> = p 2 <j> and expanding the denominator of the first term in the above equation 



1 - 



i c 



E + mr 2(E + m) 



tt-ommX X 



p 



E + m 



C 1 

- - -mR MQm x l x m U = (E-m) 
r 2 



(640) 



But, if we now perform an additional expansion assum- a cosmological constant. This leads to (the arguments of 
ing non-rclativistic velocities for the electron (which, in E indicates the value of the quantum numbers n and £) 
practice, is a good approximation), one arrives at the 

A£(1,0) = ^mR oi0m (is \x £ x m \ Is) (644) 



following expression for the term p /(E + m) 



E + m m + m + p 2 /{2m) + 



+ (641) 



As a consequence, we conclude that the non relativistic 
limit is given by 



mRo Mm -6 im (ls\r 2 \ls) (645) 



ii?oo(ls|r 2 |ls) 



(646) 



The ground state wave function is given by (at this order 
(642) of perturbations theory, there is no need to take into 
account the modification on the wave-function; one can 
where the non-relativistic Hamiltonian can be written use the standard one) 



H m cf>=(E-m) 



[106] 



p^_C 

2m r 



-mR omm x £ x m /3. 



\ls) 



1 



(647) 



Ho + H, 
1 

+ 2 



(643) 



This is our final result. Endowed with this equation, 
one can calculate, by means of perturbations theory, the 
modification of the ground state due to the presence of 



where we recall that ao is the Bohr radius. This implies 
that 

(ls pi ls \ 4f r V<V->dr = = ia l 

X 1 1 1 ™§ 7 Trag 4 

(648) 
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The vacuum Einstein equations are, see Eq. (550), = 
A off g^u and, therefore, at first order in the cosmological 
constant one Rqq = — A cff (here, of course, we only com- 
pute the corrections due to the cosmological constant and 
not those originating from the Earth gravitational field) . 
This means that 



A£(1,0) 



1 

-— mc 
2 



2 A 



ff "0 ' 



(649) 



where we have re-established the speed of light for con- 
venience. In fact, this result could have been anticipated 
from the Newtonian approach. Indeed, we have seen 
that the Newtonian potential, in presence of a cosmolog- 
ical constant, is given by — A eff mc 2 r 2 /6, see Eqs. (560) 
and (561). Using the perturbation theory with this po- 
tential exactly leads Eq. (649). Let us now evaluate the 
correction. Using that the reduced mass of the atom is 
essentially the electron mass, m ~ 9.1 x 10 _31 kg and 
that the Bohr radius is ap = 0.52 x 10~ 10 m, one obtains 



A£(1,0) 



-6 x 10 



-69 l2 



A 
Al 



eV, 



(650) 



a quantity completely unobservable, probably for ever. 
For comparison the Lamb shift is of the order 10 -6 eV. At 
this stage this does not come as a surprise. As discussed 
in the case of planet orbits, the size of the system is 
so small that the influence of the cosmological constant 
is negligible. However, requiring that the shift due to 
the cosmological constant be smaller than the Lamb shift 
leads to 



p mc < 1.6 x 10 62 /9 cri 



1.3 x 10 16 GcV 4 



(651) 



This result is less good than the one obtained from planet 
motion by about 50 orders of magnitude. Moreover, the 
theoretical expectation (516) satisfies this inequality. Let 
us notice in passing that the result based on the wrong 
regularization scheme, p vac ~ does not satisfy the 
previous bound. 

We now have a clearer view of what the cosmological 
constant problem is. On one hand, one can calculate 
the vacuum energy density by means of, apparently, well 
known techniques of regularization. These techniques 
usually lead very robust results, fully consistent with the 
observations. For instance, the radiative corrections to 
the calculation of cross-sections have proven in very good 
agreement with various high energy physics experiments. 
One the other hand, various observational results (and 
not only in cosmology) indicate that the previous calcula- 
tion is not correct. The theoretical expectation value for 
the vacuum energy density is so large that it should have 
already be seen by many experiences. Moreover, in cos- 
mology, under a certain number of hypothesis, one even 
measures p vac which is found to be many orders of mag- 
nitude smaller than the number derived from the theory. 
Since, as mentioned before, the theoretical framework 
seems to be robust (let us recall again that it involves 
calculations that can be found in the first pages of any 



textbook on quantum field theory) and since the obser- 
vations seem to be convincing, we are facing a genuine 
mystery. However, before accepting this conclusion, it is 
worth checking that no loophole is present in the previ- 
ous reasoning. This will be the goal of the next sections. 
A first question that can be asked is whether the zero- 
point fluctuations that are at the origin of the problem 
are a real physical phenomenon or just an artifact of the 
formalism of quantum field theory. We turn to this issue 
in the next section. 



XII. DO THE VACUUM FLUCTUATIONS 
REALLY EXIST? 

We have seen in the previous sections that the zero- 
point fluctuations manifest themselves as the vacuum en- 
ergy. Since this leads to conclusions that seem to be in 
contradiction with observations, it is legitimate to ask 
whether these fluctuations really exist in Nature. In par- 
ticular, we want to investigate whether the zero-point 
fluctuations could manifest themselves in other physical 
phenomena. The answer to this question is usually pos- 
itive and one experiment which is considered as a proof 
that the vacuum fluctuations are real is the measurement 
of the Lamb shift [18, 112]. In the next sub-section, we 
discuss this phenomenon. 



A. The Lamb Shift 

Let us consider the "motion" of an electron in an atom. 
Its "location" is described by the vector r. The vector r 
is going to fluctuate because of the interaction between 
the electron and the zero-point fluctuations of the electro- 
magnetic field. This interaction will slightly modify the 
position of the atomic levels leading to the Lamb shift. 
Let us study in details how this effect can be derived. 
Let V(r) be the Coulomb potential which determines the 
properties of the atom. At the point r + 5r, the potential 
can be Taylor expanded according to 

dV 1 d 2 V 
V(r + Sr) = V(r) + -^Sxi + - ^ ^ SxiSxj H . 



dxi 



2 dxidxj 



(652) 

Then, if we time average the previous expression, one 
obtains 



(V(r + Sr)) 



dV 

V(r) + ^-(Sx t 
oxi 



1 d 2 V 

2 dxidx 



(SxiSxj) 



(653) 

and, because the problem is spherically symmetric, we 
have (Sxi) = and (SxiSxj) = (8r 2 )5ij/3. The first term 
in Eq. (653) is the unperturbed Coulomb potential and 
the last one (since the second one vanishes) can be viewed 
as a small perturbation of the Hamiltonian, namely 



AH 



1 d 2 V 

6 dxj dx, 



(Sr 2 ) 5 K 



(654) 
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In an atom with atomic number Z, the Coulomb poten- 
tial is given by 



V 



Za 



(655) 



where a ~ 1/137 is the fine structure constant. This 
means that 5ijd 2 V/(dxidxj) = ATrZa5(r). 

Let now us assume that the atom is placed in the 
quantum state \£,m). As we have already discussed in 
Sec. XI C, at first order, the time-independent perturba- 
tions theory tells us that the correction to the energy of 
the level {£, m) is given by [57] 



AE(n,£) = (£,m\AH\£,', 



(656) 



Using Eq. (654), one obtains the following expression for 
the energy displacement 

AE(n,£) = ^-Za\y ne (0)\ 2 (Sr 2 ) , (657) 

where ty n £ is the wave-function in the corresponding 
quantum state. Since the wave-function at the origin 
can be expressed as [57] 



l*n*(0)| 



1 



.1/2 



mZa 



3/2 



(658) 



where m is the electron mass (m ~ 0.511 MeV) one finally 
arrives at [18] 

AE(n,£) = ^^^(Sr^6 i0 . (659) 
12 n a 

Therefore, we have reduced the problem to the calcu- 
lation of the quantity (<5r 2 ). In order to evaluate this 
quantity, we must now evaluate the displacement of the 
electron. The electron "moves" under the influence of 
the electric field within the atom and, as a consequence, 
Sri obeys the following equation 



(660) 



where E^ = u v = —Ai is the electric field (or, rather, 
its vacuum fluctuating component). Using the usual ex- 
pansion of the vector potential in terms of creation and 
annihilation operators, see Eq. (395), one can write 



Ei(t,x) 



1 



(27T) 



3/2 
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In the same manner, one has to Fourier expand Sri which 
is considered as an operator in this context. One obtains 



Sr t (t,x) 



1 



dfc]Te?(fc) 



\3/2 

+Sr*(k,t) (a%Ye 



5r{k,t)a%e ikx 

(662) 



The Fourier amplitude of the displacement operator can 
be calculated from the equation of motion (660). This 
leads to the following expression 



Sr{k,t) 



1 



(663) 



From this expression, one can now evaluate the quantity 
(<5r 2 ). One finds 



(Sr 2 ) 



(0 \S ij SnSrj\0) 

i r 2 

33 y dfe^^e?e^r(fe,t)<5r*(fe,i) 



(27T) 

2a 



(664) 
(665) 



where we have used a = e 2 / (47r). Therefore, inserting the 
above formula into Eq. (659), one obtains the following 
equation for the energy shift [18] 



AE(n,£) 



4mZ 4 a 5 
37m 3 



— 0£ 

U! 



(666) 



This equation is the main result of this section. It gives 
the energy levels displacement due to the presence of the 
electromagnetic zero-point fluctuations. 

As usual this expression is divergent and needs to be 
regularized. In a simplified treatment, one simply as- 
sumes, see Ref. [18], that the wavelengths in the above 
sum must be larger than the Compton wavelength of 
the electron, which implies lo < w max ~ m and that 
u > w m in ~ 1/ao — l/[l/( aTO )] — ma i where ao is the 
Bohr radius. In this case, the previous results reduces to 



AE(n,£) 



4mZ 4 a 5 1 / 1 
37m 3 V a 



(667) 



Let us now evaluate this quantity for the Hydrogen atom. 
Clearly, there is no effect for n — 1. But for n = 2 and 
£ = 0, one has 



. „. . 4ma 5 



668 MHz, 



(668) 



where we have used m ~ 0.511MeV = hv with h = 
4.135 x 10 _15 eV x s. 

This shift has been experimentally observed and is usu- 
ally taken as a proof that the zero point fluctuations are 
real since they lead to an observed physical phenomenon. 
Based on this result, it now seems difficult to argue that 
the zero-point fluctuations are just an artifact of the for- 
malism of quantum field theory. Moreover, there exists 
another experiment, the Casimir effect, where the effects 
of the vacuum fluctuations can be observed. In the next 
section, we investigate this case. 
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FIG. 8: Sketch of a Casimir cavity of width L. 



B. The Casimir Effect 



tions, namely 



In this section, we study the Casimir effect [113-119]. 
Since its discovery, this effect has been studied in great 
details and it is clear that, here, one cannot review all the 
literature. We will give a description of this phenomenon 
in a simplified context only and will refer to Refs. [113— 
115] for readers interested in learning more about this 
interesting subject. 



1. The Casimir Force 

Let us consider an experimental situation where we 
have two conducting plates in the (x, y) plane separated 
by a distance L long the z direction, see Fig. 8. For sim- 
plicity, instead of considering the electromagnetic field in 
the cavity (which is the case in the real world, especially 
when measurements are performed), we will just treat the 
case of a real massive scalar field $(i, x). This field obeys 
the Klein-Gordon equation already studied and solved 
before, see Eq. (60), except that now the boundary con- 
ditions are modified by the presence of the two plates. 
Let us write $(t,x) = X (x)Y (y) Z (z)T \t) . Then, in- 
serting this anzatz into the Klein-Gordon equation, one 
obtains 



f 1 d 2 X 1 d 2 Y 1 d 2 Z 
1 1 1 

T X dx 2 



Y dy 2 



Z dz 2 



m 2 = 0, (669) 



where m is the mass of the scalar particle. Clearly, one 
can separate the variables and, as a consequence, the 
above equation can be split into four differential equa- 



T 2 
— — m 
T 

1 d 2 X 

X~dx 2 ' 

1 d 2 Y 

Y dy 2 

1 d 2 Z 



k' 



(670) 
(671) 
(672) 
(673) 



Z dz 2 

where the constants C, k x , k v and k z are related by 

C - k 2 x - k 2 - k 2 z = 0. (674) 



The solution in the x and y direction can be expressed 
in terms of standard plane waves, X(x) = A x e lk:cX + 
B x e~ lk * x and Y(y) = A y e lk y y + //,,(■ The quan- 

tities A x ^ y and B x y are just integration constants. Of 
particular interest is of course the solution along the z- 
axis. It can also be written as 



Z{z) = A 2 



+ B z e 



-ik z z 



(675) 



where A z and B z are two integration constants. But 
the difference with the case of a free field is that the 
boundary conditions "feel" the presence of the two plates. 
These boundary conditions arc given by Z(0) = and 
Z(L) = 0. The first one implies that A z = —B z while 
the second one gives 



or, equivalently, 



2ism{k z L) = 0, 



k z = n— , 



(676) 
(677) 



G4 



where n is an integer. Finally, the equation for the func- In this expression, k x and k y take continuous values, 
tion T(t) can also be easily solved and the solution reads 



T(t) = A t e iu;t + B t e~ lult , 



(678) 



where lu 2 — C 2 + m 2 . Therefore, upon using Eq. (674), 
the frequency u takes the form 



k 2 + k 2 + n2 -j^2 + m2 - 



Based on the previous considerations, one can now 
write the field operator as an expansion in terms of cre- 
ation and annihilation operators. But this expansion en- 
codes the fact that, along the z-axis, the wave-number is 
(679) discrete. Concretely, we have 

I 



dk, 



(2ny/ 2 J (2ny/ 2 ^\ 



dk v / 2 1 

7172 2^/-^ 



in (njz) (c k ^ kv , n c-^ t+lk * x+lk yy + cl x!ky J ut - ik * x - ik ^ . (680) 



The operators Ck x ,k y .n and k ^ n obey the usual commutation relations for bosons. It is interesting to check that the 
above expansion is indeed consistent with the commutation relations between the field and its conjugate momentum. 
For this purpose, let us compute the following quantity 



/ Ax j d„jf d, (mf.) =/^/A|/d./d» 
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Ck x ,k y 



-iuit+ik^x+ikyy 



+ C fe k n e *" 



t-ik*x-ik yV \ e -rp*x- PyV J dz gin (m,l Z ^ sin (thj-z) (681) 



Using the fact that J Q L dzsm(nirz/L) sia(rmrz/L) = (L/2)5 mn (as can be checked directly) and J dxe^ fex p ^ 1 
(2n)5(k x — p x ) one arrives at 



J dx J dy dz$(t, x) e-v^-vyv sin (m^z) = 



-iuit , t „iurf 
Px,— p y ,m 



(682) 



In the same manner, one can also evaluate the following integral involving the conjugate momentum Tl(t, x) = <f>(t, x) 
(a dot means a derivative with respect to time) 



J dx J dy J dzll(i, cc) e 



ip x x-ip y y ■ 



7T 

sm | m—z 



2tt. L 
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From the two above expression, one can deduce a formula expressing the creation and the annihilation operators in 
terms of the field operator and its conjugate momentum. One obtains 



■T. ,P V ,m 
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we" m 
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- j dxdy j dzH^ p " x+ip y y sin (m^z 



We are now in a position where we can compute the commutator of a creation and an annihilation operators. It can 
be written as 



3 Px,P H ,™: c L,«? B ,m = (2-k) 2 L — J y 7 J V J 

+—n(t,x)e~ 



$(t,x)e- ip * x - ip y y sin ( n-z 

L 



— II(i, x) e - ip * x - ip yy sin ( n^z) , $(t, ») e^*-*** sin f m ~z 

\ L J \ L 



■H(t, x) e^+^f sin (m^z) 



(684) 
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Then, we make use of the canonical relation [$(t, x), H(t, y)} = iS 3 (x — y) in the previous formula. This leads to the 
following complicated expression 



(2tt) 2 ~ 



dsdy / dz / dxdy / dz< c 



-ip x x-ip y y+i 



v^+iivv sin 



sin ( m-z ) [$(t, a:), n(i, sb)] + _ e - < ** x -*»»+ < *» ie + < 9*» sin ( n-z ) sin ( m- 



[n(t,a:),$(i,a)] 
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L/ to z 



(686) 



We have thus reached our goal, namely we have shown 
that the canonical commutation relation between the 
field and its conjugate momentum implies the usual com- 
mutation relation between the creation and the annihila- 
tion operator. 

We are now in a position where we can evaluate the 
energy density of the field in the Casimir cavity. It is of 
course given by the same expression as in the free case, 
see Eq. (67). This means that one should first calculate 
the three terms that participate in the expression of Too, 
see Eqs. (64), (65) and (66). Let us start with (0|$ 2 |0). 
Using the Fourier expansion of the field and the commu- 
tation relations established above, one obtains 

(oi$ 2 io> = i^i^f-ii- 

x 1 17 / 2tt / 2tt ^ L 2uj 

J J n=l 

x sin 2 (n^-z) . ( 687 ) 
and, if we integrate over the cavity volume 

where D 2 is the area of the plates (of course, the two par- 
allel plates have the same surface). In the same manner, 
one finds 



(688) 
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Finally, the last relevant term reads 

Lj*^<!t! A -t±b <«»> 

Putting everything together and using the definition (67), 
the expression of the energy of the field within the empty 



cavity takes the form 



da; (0|T 00 |0) 
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Of course, this expression has exactly the expected form. 
Roughly speaking it can be written as iw" , see 
Eq. (679), the only difference being that, along the z-axis, 
the sum takes into account the fact that the momentum 
is discrete. Another common point with our previous 
considerations is that the above number is actually di- 
vergent and must be regularized. In the following, we 
assume m = for simplicity and use dimensional regu- 
larization. Therefore, the expression of the energy can 
now be written as 



E = 



y n—l 

y n—l 

d cy d/2 P °° 1 
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2tt J r(d/2) 



V- 
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where d is the "dimension" of the plates (i.e. d = 2 in 
the real world) and k± denotes the wave-vector living 
in the sub-manifold defined by the plates. Then using 
the change of variables y = Lk± / (wr) and the following 
definition [48, 49] of the Euler's integral of first kind 



B(x,y)=2 
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the energy can be re-expressed as 

(D\ d *W r s* +1 r(-d/2-l/2) d+1 

W 2 UJ r(-i/2) ^ n • 

(697) 

The above formula can be expressed in term of the Rie- 
mann zeta function defined by the following expression 



(698) 



n=l 



Moreover, upon using the two following equations 
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one can re- write the energy as 



E 
D 3 



The Casimir force arises because of a shift in the vacuum 
energy when the plates are present compared to the sit- 
uation where they are absent. So the relevant quantity 
is in fact AE/D d = (E — E Q )/D d where E /D d is the 
above quantity in the limit L — > +oo. As a consequence, 
for d = 2, one arrives at 



AE _ 7T 2 

L) 2 " _ ~ 1440L 3 ' 



(702) 



where we have used £(4) = 7r 4 /90. Therefore, the corre- 
sponding force is given by [113-116] 
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dL \D 2 J 480L 4 : 



(703) 



which is one half the result obtained in the case of the 
electromagnetic field since we have only one state of po- 
larization instead of two. We notice that the force is 
attractive (this could change in the case of more compli- 
cated geometries). 

It is also interesting to check that the final result is 
independent from the regularization scheme used to ob- 
tain it. In order to test this idea in the simplest situation 
(i.e. where the calculations are easy), let us consider again 
the energy of a Casimir cavity but, this time, with d = 0. 
In this case, our zeta-function regularized result reads 



E 



i r(i)C(2) 

L 4tt 



= — — , (704) 
24L 



since T(l) = 1 and £(2) = 7r 2 /6. One the other hand, we 
can return to the expression of E in Eq. (697). In the 
case d = 0, the energy density p{L) = E/L reads 



P(L) 



7T 

2L2 



(705) 



In order to regularize this divergent quantity, we now 
proceed differently and introduce an exponential cut-off 
such that (previously, we argued that introducing a cut- 
off was not a good method because this breaks Lorentz 
invariance. Here, of course, this objection is no longer 
valid as the presence of the two plates along the z-axis 
obviously breaks translation invariance; as a consequence 
working with a cut-off seems legitimate in the present 
context) 



/L 



The sum is easy to perform and one obtains 



p(L,a) = ^sinh- 2 (|L 

7T 7T 
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Now, the regularized Casimir energy is the shift energy 
in the vacuum after having sent the cut-off to zero; in 
other words 



E = L lim p(L,a) — lim p(L,a) 

a— >-0 L L— >oo 

~24L' 
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that is to say exactly the same expression obtained with 
another regularization scheme. It is reassuring to notice 
that the final result is independent of the regularization 
method used. 



2. The Casimir Stress-Energy Tensor 

In this section, we aim at calculating the stress energy 
tensor of a scalar field in a Casimir cavity [113-116, 119]. 
Our goal is of course to compare this tensor with the vac- 
uum stress-energy tensor = —p mc g^ iu . Here, we use a 
method based on the calculations of the Green function. 
This is yet another method of calculation of the Casimir 
effect, different from the two ones exposed before and, 
therefore, it will be interesting to compare this approach 
with the treatments used previously. 

Let us start with the calculation of the Green function. 
It obeys the following equation 



fjT^8 v G (x,x) = S i4) (^ - a") . 



(711) 



It is interesting to compare the above formula with 
Eq. (482). It is clear that this is the same equation but 
with g^y = rj^n (flat space-time) and m = (mass-less 
field). In order to solve this equation, it is convenient 
to Fourier transform G{x 1 x) in frequency and transverse 
momentum. Explicitly, one writes 



G(x, x) 



duj dkj__ _j (1) (t_t)+i* a! (a!-x)+ifc v (i/-p) 
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Inserting this expression into the Green equation (711), 
one obtains 



(713) 



dz 2 

where we have defined 
A 2 = to 2 - 



Ti+A g(z,z) = S(z- z) 



(714) 



Our goal is now to determine explicitly the Green func- 
tion. Let us start with the solution within the plates, 
i.e. < z,Z < L. When z ^ ~z, the solution can be 
written 



g(z, z) = A(z) cos (Xz) + B(z) sin (\z) , 



(715) 



where A(z) and B(z) are two unknown function of z. We 
must now take into account the boundary conditions. If 
z < ~z, the condition g(z = 0,z) = implies that A = 
and therefore g = Bs'm(Xz) (for z < z). On the other 
hand, if z > ~z, the condition g(z = L,~z) = means that 



g = A cos(Az) 



. cos(AL) . . 
A . ; : sin (Xz , 
sm(AL) 



z > z. 



(716) 



We have two unknown quantities and, therefore, we need 
two equations. One is provided by writing the continuity 
of the Green function when z = ~z. The other is obtained 
by integrating Eq. (713) from ~z — e to ~z + e. This gives 



_ f dz 2 
which leads to 
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Therefore, taking the limit e — > 0, the solutions of these 
two equations (i.e. the one expressing the continuity of 
the Green function and the one we have just derived) 
read 

A{z) = ~sin(Alz), (719) 
A 

B(z) = Icos(A^)-^M^). (720) 
A sm(A-L) 

As a consequence, the Green function between the plates 
can be expressed as [116] 



g{z,z) 
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z < z, 



z > z. 



(721) 



The same method can be used to determine the Green 
function outside the cavity. Let us apply it for z and z 



larger than L, i.e. on the right hand side of the cavity. 
In this case, it is more convenient to re-write the general 
solution (715) as 



g{z,z) = C(z)e^ z + D(z)c 



-iXz 



(722) 



When z > ~z, we require only one branch to be present 
which amounts to take D = 0. For z < z, the bound- 
ary condition is still g(L,~z) = 0. Then, straightforward 
manipulations lead to [116] 
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! iA( >- L )sin(Az-AL), z < z, 
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which completes our determination of the Green func- 
tion. The calculation of the Green function outside the 
cavity, but on the left hand side, proceeds exactly in the 
same way. The corresponding expression is obtained from 
the one above with L — 0, reflecting the fact that the 
boundary condition is now g(0,z) = 0. 

Endowed with the Green function of the problem, one 
can now turn to our main goal, namely the calculation of 
the stress-energy tensor. We have already seen that, see 
also Eq. (101), 



G(r,r) = i (0 |T [$(r)$(r)]| 0} 



(724) 



But we know that the Too component of the stress energy 
tensor can be expressed as, see Eq. (7), 



Too = \ [(5 $) 2 + (&$) 2 + (d y $>) 2 + (c\$) 2 ] . (725) 



Therefore, if we apply the following operator to the Green 
function in Eq. (724) 



[<9 (% + d x dx + dydy + d z dz] , 



(726) 



where c^j means a derivative with respect to x^ 1 , to 
G(r,r) and, then, take the space-time points to be the 
same r = r, one should obtain an expression for (Too). 
Using Eqs. (712) and (714), the result reads 



(T„„) - i / 



du dk x dk. 



2ir (2tt) 2 
xg(z,z) 



k 2 x + k 2 y + d z dz) 



(727) 



Then, we use the expression of the Green function inside 
the plates (721) and one arrives at 
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(r 00 ) = ttt 



1 f duj dk x dkv 



1 



2i J 2tt (2tt) 2 A sin (AX) 
1 f duj dk x dk y 1 



(lo 2 + k 2 + fej,) sin (Xz) sin (AX — AX) — A 2 cos (Xz) cos (AX — Xz) 
\uj 2 cos (AX) - k 2 cos (2Xz - AX)] . 



2i J 2tt (2tt) 2 A sin (AX) 

In order to evaluate this integral, we perform the two Wick rotations u — > iQ and A — > in. This leads to 

1 



Irp \ I TdC, f dk x dk y 



[C 2 cosh (kX) + /c 2 cosh (2kz - «X)] . 



2J 2tt y (2?r) 2 Ksinh(KX) 
Then, we introduce polar coordinates in the plane (£, A;), i.e. £ — kcos8 and k = ns'md. This gives 



(T ° 0> - 

1 



tt/2 



d<2^ 



sm( 



sinh (kL) 



[cos 2 cosh (kL) + sin 2 cosh (2nz — kL) 



12tt 2 



dKK 3 



cosh(ftX) i cosh(2«;z — kL) 



sinh(ftX) 



sinh(«;X) 



r 



(728) 
(729) 

(730) 

(731) 
(732) 



Finally, writing that 

cosh(ftX) 
sinh(KX) 



-r 



(733) 



the expression of (Too) can be re-expressed as 



(Too) = - — 



IHH 



1 

6tt 







1 



(734) 



The last step of the calculation consists in evaluating 
the last term [we call it g(z) in what follows; not to be 
confused with the Green function] of the above equation. 
One has 



i i 

"oTT 2 16X 4 

1 1 
"OTT 2 16X 4 

+c 



dyy 3 J : (735) 



n = 

a y(l—n—z/L) 



°° poo 

£ / d ^ 



C-f - 1 



j){z/L-n) 



(736) 



where we have expanded the denominator in a geometric 
series. The two above integrals can be expressed in terms 
of the Hurwitz zeta function defined by [48, 49] 



1 



m=Q 



(m + a) 



and this leads to 
1 



g(z) = 



16tt 2 X 4 



C 4,- +C 4,1-- 



(737) 



(738) 



time-time component of the stress energy tensor, namely 

7T 2 1 



(To, 



1440X 4 6tt 2 



o 



dn—+g(z). (739) 



Of course we notice that the first term in the above equa- 
tion is similar to Eq. (702). 

Let us now evaluate (Too) outside the cavity. The cal- 
culation proceeds along the same lines but we now need 
to use the Green function given by Eq. (723). Straight- 
forward manipulations lead to 



(Too) 



1 

e^ 2 



oo ^3 

dn- 



2 16^ 2 ( z - L) 



(740) 



We obtain the same structure as inside the cavity, namely 
a divergent term and a surface divergence term, i.e. a 
term which is divergent only on the plate. As explained in 
Ref. [116], this type of terms can be removed by restoring 
conformal invariance. 

Let us now calculate the quantity (T zz ) inside the cav- 
ity. Clearly, following the previous considerations, it can 
be written as 



^ = h I 



1 f dco dk x dk y 



2ir (2tt) 2 

xg(z,z) 



(dt&t - d x dx - dydy + d z &z) 



(741) 



Following exactly the same procedure as before, we arrive 
at 



(T zz ) 



1440X 4 6tt 2 



dn- 



(742) 



If we use the fact that J Q °° k 3 / (e 2Ki - l)dn = Let us now calculate (T zz ) outside the plates. As we did 
7r 4 /(240X 4 ), one arrives at our final expression for the for the time time component, we must now use the Green 



G9 



function (723). This gives 

/T zz ) = - [ du kxdky ° iKz ~ L) y 



2i J 2tt (2tt) 2 A 
+i cos (Xz — XL) 



3 



-cIk. 



sin (Xz — XL) 

(743) 
(744) 



Of course, this quantity is divergent. 

Finally, in order to complete the determination of the 
stress-energy tensor, one must calculate (T xx ) = (T yy ). 
It can be expressed as 

( T xx) = j. J ^^0fr ( d A + d x dx- dydy- d z &z) 



xg(z,z) 



(745) 



dw dk x dk y 



1 



2tt (2tt) 2 Asin(AT) 
x sin(Az) sin(AL — Xz) 

+X cos(Az) cos(AT — Xz) 



(746) 



A first step consists in dealing with the terms propor- 
tional to k 2 — k 2 . In fact these terms are of the form 

y 

J dk±(k 2 — k 2 )J r (k±), where T(k±) represents the other 
terms in the integral which are function of fc 2 + k 2 only. 
Therefore, one can evaluate them by going to polar co- 
ordinates in the plane (k x ,k y ). This leads to an expres- 
sion of the form J dk^dtpk^cos 2 ip — sin 2 tp)T(k±) = 0, 
which shows that the corresponding contribution van- 
ishes. Then, following the same steps as before, one re- 
duces the above expression to 



(T xx ) — 



1 



12tt 2 



sinh(/-tT) 



cosh(2«;z — kL) 



-2 cosh(re,2) cosh(«T — kz) 



= -(T 00 ), (747) 



where we have used that 2 cosh(Kz) cosh(fiT — kz) = 
cosh(/tT) + cos1i(2kz — kL). Of course, as already men- 
tioned, this also shows that (T yy ) = — (Too). All the other 
components being zero, this completes our calculation of 
the stress energy in the cavity. Only the calculations 
of (T xx ) = (T yy ) outside the plates remains to be done 
and we now turn to this question. Considering again 
Eq. (745), we have 



it \ 1 f duj dkxdk y e lA(z ~ L) 

{Ixx) ~ 2l J 2n (27r) 2 A 



uj 2 sin(Az — AT) 



-iX 2 cos(Az — AT) 



= -<T 



oil/ 



(748) 



We are now in a position where we can write the full 



stress energy tensor. Inside the cavity, it takes the form 



(Pvac +Pc™) 



/l ON 
0-100 
0-10 

\0 3y 



-g(z) 



a o o o\ 

o-i o o 

0-10 

^oo oo/ 



(749) 



where we have defined the two quantities p and p c 
by 



1 



Pvac 



dn — 

o 2 



1440T 4 



(750) 



Outside the cavity (more precisely on the right hand 
side), the results derived above imply that 



( T A^)out S1 do = Pv„, 



a o o on 
o-ioo 

0-10 
i0 3; 



16tt 2 (z - L) 4 



(\ ON 
0-100 
0-10 

Vo o 0/ 



(751) 



As discussed in Rcf. [116], the terms proportional to g(z) 
and (z — T)~ 4 are surface divergent terms. They can 
be removed in the case where conformal invariance is re- 
established and, for that reason, are not present in elec- 
tromagnetism, i.e. when one considers an electric field 
inside the cavity rather than a scalar field. In this case, 
we see that the difference between the stress-energy ten- 
sors inside and outside the cavity reads 



A(T^) 



/l ON. 
0-100 
0-10 

V0 3y 



(752) 



which is a finite result. The infinite part has been re- 
moved by subtraction. Finally, the stress energy tensor 
can be written in a covariant form, namely [119] 



air (43 



■<T), 



(753) 



where z^ in the unit vector in the z direction. 

At this point, several comments are in order. Firstly, 
we notice in Eqs. (752) and (753) that the energy den- 
sity (i.e. the time-time component) is similar to the ex- 
pressions already derived in Sec. XII B 1 by other meth- 
ods. This confirms that this calculation is independent 
of the regularization scheme used to obtain the result. 
Secondly, it is clear that the stress-energy tensor (753) 
is not similar to the cosmological constant stress-energy 
tensor T^ = —p mc g^ u . This illustrates the difference 
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between the Casimir case and the cosmological constant 
case. In the latter case, one measures the absolute value 
of the vacuum energy while, in the former case, one is 
only sensitive to the difference of vacuum energy from 
one side of the plates to the other. 



Thirdly, the Casimir Force has been observed in the 
laboratory and this is usually taken as another evidence 
in favor of the existence of the zero-point fluctuations. 
After all, since we observe the Casimir force and since 
this force is due to a change in the structure of the vac- 
uum that is perturbed by the presence of the plates, 
the fact is that the vacuum must exist as a real phys- 
ical phenomenon and is not an artifact of the quantum 
field theory formalism. In addition, this also validates 
the methods of regularization used to calculate p vac . In 
the context of the Casimir effect, these methods seem 
to work very well and to lead to predictions that are in 
agreement with the experiments. Therefore, why would 
the very same techniques, used in a similar context, fail 
to regularize the cosmological constant? 



Fourthly, the claim that the calculation and the ob- 
servation of the Casimir force are evidences in favor of 
the reality of vacuum fluctuations has been challenged 
in Refs. [120, 121]. It is indeed possible to derive the 
expression of the Casimir force without referring at all 
to zero point fluctuations, by means of the "source the- 
ory" approach to quantum field theory. However, to our 
knowledge, there is no attempt to do the same for the 
Lamb shift (but this would certainly be an interesting 
exercise). Therefore, it seems that zero-points fluctua- 
tions are nevertheless necessary to quantum field theory 
even if the above remark should lead us to tone down the 
claim that they have been seen in the laboratory. Let us 
also notice that, very recently, another situation where 
the vacuum fluctuations can lead to an observable effect 
has been studied in Ref. [122]. This article investigates 
whether the quantum zero point fluctuations could cause 
the rotation of a small nano device. This "rotational vac- 
uum effect" has in fact a similar origin than the Casimir 
effect. 



Finally, it should be clear that, if the Lamb shift and 
the Casimir effect seem to indicate that the vacuum fluc- 
tuations are real, they do not say anything with regards 
to their gravitational properties. In other words, they 
do not say how the vacuum fluctuations weigh (it is 
worth noticing here that the gravitational properties of a 
Casimir cavity have been studied in Refs. [123-129]). It is 
clear that there exists the possibility that the zero-point 
fluctuations do exist but have non-standard gravitational 
properties, i.e. that they abnormally weigh. It is there- 
fore interesting to investigate if there is an experiment 
which measures whether this is true or not. In the next 
section, we turn to this question. 



XIII. DO THE QUANTUM FLUCTUATIONS 
GRAVITATE? 



In this section, we study the gravitational property of 
the Lamb shift [12, 130]. We will follow the treatment of 
Ref. [130]. In Sec. XII A, we showed that the interaction 
of the electron with the "vacuum electric field" inside 
the atom leads to a shift of the atomic levels. In some 
sense, this means that the zero-point fluctuations affect 
the mass (energy) of the atom. Therefore, by studying 
how the atom falls down, one can put constraints on the 
gravitational properties of the vacuum fluctuations. This 
is the idea that we pursue in this section. However, since 
the corresponding effect is small for an atom and is larger 
for a nucleus (of course, the Lamb shift also exists in this 
case), it is better to study the latter case. For this reason, 
we start with recalling some basics facts about nuclear 
models. 

The main idea is that the interaction between a given 
nucleon and the other A — 1 nucleons can be mimicked 
by a Woods-Saxon potential given by 



^w B (r)=-Vb 



1 — exp 



r- R 



(754) 



where a describes how abrupt the potential is while R 
is the radius of the nucleus. As is well-known one can 
take R = r^A 1 / 3 where r ~ 1.3 x 10~ 15 m = 1.3 fm. 
The quantity Vq represents the depth of the potential 
and can be taken to be Vo — 45MeV. In fact, in order 
to have a reasonable model, one must add to the Woods- 
Saxon potential a spin orbit interaction term which can 
be expressed as 



-f(r)L 



(755) 



where L is the angular momentum of the nucleon and 
s its spin. With a good approximation, one can write 
f(r) ~ 24A~ 2/3 MeV. Then, the Hamiltonian of the nu- 
cleus can be written as 



H 



A 

E 

3=1 



^ r +V ws (r j )-f(r j )L j 



(756) 



The total wave-function can be taken as the product of 
the individual wave functions (in fact, more precisely, one 
should consider a Slater determinant in order to take into 
account the fact that the total wave function must be 
anti-symmetrized) , that is to say 



*(*•!,••• ,r A ) \[ \' :>,.:• 



(757) 



Each individual wave-function can be obtained from a 
Schrodingcr equation given by (for the moment, one ig- 
nores the spin orbit term since it can be treated as a 
perturbation) 



-^-v 2 + y ws (r,) 



* j {r j )=E j * j {r j ), (758) 
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FIG. 9: Energy levels of the nucleus. On the left-hand side are represented the levels in the approximation where the Woods- 
Saxon potential is approximated by a parabola. In this case, each level (which includes a series of sub-levels) is separated by 
the energy Tilo. On the right hand side are represented the levels where the spin orbit interaction (755) is taken into account. 
This breaks the degeneracy of the sub-levels and reproduces the magic numbers. 



(759) 



where (i is the reduced mass. 

11 1 
fjj m p (A — l)m p ' 

or /i ~ m p (A — I)/ A ~ m v ~ 1.6 x 10 _27 kg, m p being 
the proton mass. Unfortunately, the Schrodinger equa- 
tion cannot be solved analytically for the Woods-Saxon 
potential and, therefore, we are forced to make approxi- 
mations. The simplest one is just to describe the Woods- 
Saxon potential with a parabola, i.e. 



(760) 



with ujq = ^/2Vb/(/ii? 2 ). In this case, the question be- 
comes solvable. We can use the spherical symmetry of 
the problem and write the wave function as 



A r j) 



Ye™ (8j,(Pj) , 



(761) 



where Yg m is a spherical harmonics [48, 49]. The radial 
function u n { is controlled by the equation 



2n dr 2 



1 



1 



2fir 



1 1 

fiuj Q r 



u n ,e 



(Ej + V )u n .i, 



(762) 



which can be solved in terms of Laguerre polynomials 
L" [48, 49]. Explicitly, one has 

*(r,) = Cry^V^L^T (^V) n™ {9 ilVi ) , 

(763) 

where C is a normalization constant to be determined. 
In the above expression, one has n > 1 and I > 0. The 
energy levels are given by the expression 



E(n,£) = -Vo + huo 



2(n-l) + £+- 



(764) 
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The factor 3/2 = 3 x 1/2 originates from the fact that 
we have a three-dimensional harmonic oscillator. The 
quantum numbers have been chosen so that the ground 
state is n = 1, I = which allows the spectroscopic 
notation "Is", see Fig. 9. 

The only thing which remains to be done is to deter- 
mine the constant C . Obviously, this is done by normal- 
izing the wave-function, i.e. 



J = 1, (765) 



Using the solution obtained above and performing the 
change of variable p — (fiuio/h)r 2 , the previous condition 
reduces to 

JCf fiiuoyt-a/i 
2 V h J 

POO 

x / ^V^)^) = 1. (766) 



Then, using Eq. (7.414.3) of Ref. [49], one obtains 

|C|^2(n-l)!(^) f+3/2 r-(n + £ + i). (767) 

We are now in a position where we can compute the value 
of the wave- function at rj = for I = 0. This gives 



2 1 

47r' 



(768) 



Using the fact that L^O) = (\ J^) 
l/2)/[V7r(n- l)!/2], one obtains 



I* 



n,£=0 



(0)| : 



2 r(n + l/2) /^o\3/2 
7T 2 (n-1)! V ft J 



(769) 



This completes the first part of the calculation. 

We have established on very general grounds that the 
shift in energy due to the interaction of a charged particle 
(in Sec. XII A, the electron in the atom, here the proton 
in the nucleus) is given by Eq. (657). The means value 
of the square of the displacement has been calculated in 
Eq. (665). Applying these formulas to the case of the 
protons in the nucleus, one obtains [130] 

aE( „, 0) ^ ?z <.) |.JUp)fjT£. <™> 

As usual, this infinite integral must be regularized and 
usually leads to logarithmic corrections, see for instance 
Eq. (667). We will ignore this factor in what follows 
since this does not affect too much the final result. As 
explained in Ref. [130], the sum over "(«)" runs only over 
s-wave since this is the only way to get \& n ^(0) 7^ 0. 
Moreover, still following Ref. [130], it is reasonable to 
approximate the field experienced by a proton by a field 
produced by the other protons in inner shells. This is 



why Z^> represents the total number of protons present 
in shells inner than the proton "(i)" under consideration. 

Our goal is now to test how the vacuum fluctuations 
weigh. For this purpose, we write the gravitational mass 
of a given nucleus as [130] 



rjAE, 



(771) 



where m; is the inertial mass. The parameter 77 is a priori 
a free parameter which controls to which extent the weak 
equivalence principle would be violated. If rj = 0, there is 
no violation of the universality of free fall. A convenient 
measure of potential violations of the weak equivalence 
principle is given by the Eotvos ratio 



ra g (l) _ m g (2) 
mi(l) m(2) ' 



(772) 



for two bodies 1 and 2. Using the previous calculations, 
this ratio can be expressed as 



4a 2 
'3to^ 



« |viv(2,i) 



*£fio(0) 



o(0) 



(773) 



where we have used m\ ~ Am p , a relation valid since we 
work at first order in the parameter 77. Then, we use the 
fact that 



(/.two 



.3/2 



(2m v V \ 3/4 



1 



A 1 / 2 



(774) 



to express the Eotvos ratio as [130] 



4a 2 

7/(1,2) = 7/ — 



2U 



3/4 



.4 



1 y y(i) r(n + l/2) 
A\ /2Z ? 1 (n-1)! 
1 {i) T(n + 1/2) 



(n-1)! 



(775) 



Let us evaluate this number for aluminum (Z = 13, 
A = 27) and Platinum (Z = 78, A = 195). In the case 
of aluminum, since we only have 13 protons, the only s- 
shell protons are on the 1 Si/ 2 level (in fact, we have two 
protons in lsi/2, four protons in I753/2, two protons in 
lpi/2 and five protons in 1 d 5 / 2 > scc Fig. 9). As a conse- 
quence, there is no inner protons for those protons (since, 
obviously, there is no inner level than the ground state 
1 S1/2) an d the corresponding contribution in the Eotvos 
ratio vanishes. In the case of Platinum, all the levels 
arc filled up until the level 1 hn/2 which contains height 
protons, see Fig. 9. Therefore, we have two protons on 
1 s 1 / 2 ; two protons on 2 S1/2 and two on 3 S1/2, see again 
Fig. 9. This means that the 3 S1/2 protons see four inner 
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protons and the 2 s\/2 protons see two inner protons. As 
a consequence, one can write 

(776) 

Using the values of Vq and ro chosen previously this leads 
to 

J7(A1, Pt) ~ -5.6 x lO" 10 r) (777) 

Given the fact that |j?(Al,Pt)| < 10~ 12 [131, 132], we 
reach the conclusion that 

77<1.7xl0 _3 . (778) 

Therefore, the vacuum fluctuations seem to have stan- 
dard gravitational properties to a very good accuracy. 

The conclusion of this section is that it seems difficult 
to argue that the zero-point fluctuations do not gravitate. 
Clearly, this makes the cosmological constant problem 
more acute. However, the results obtained here are by no 
means a final proof. For instance, it seems reasonable to 
assume that the stress energy tensor of the Lamb shift is 
not of the cosmological constant type. Therefore, the fact 
that a nucleus falls down normally does not necessarily 
imply that the cosmological constant weighs according to 
the Einstein equations. Nevertheless, no sign of inconsis- 
tency appears in the above situation and one is obviously 
tempted to assume that zero-point fluctuations couple to 
gravity as any other type of matter. Again, this makes 
the cosmological constant problem more mysterious. 

XIV. THE WEAK EQUIVALENCE PRINCIPLE 
IN QUANTUM MECHANICS 

In the previous sections, we have argued that the zero- 
point fluctuations seem to be real and seem to gravitate 
normally. Therefore, these two reasons cannot be in- 
voked to avoid the cosmological constant problem. How- 
ever, there is another (but related) issue that it is in- 
teresting to investigate. When we consider the vacuum 
stress energy tensor in the Einstein equations, we im- 
plicitly assume that the zero-point fluctuations obey the 
weak equivalence principle, i.e. that they "fall" in a grav- 
itational field as any other type of matter. Clearly the 
zero point fluctuations are of quantum-mechanical origin. 
The weak equivalence principle is well established (and 
tested) in classical physics and general relativity is based 
on this principle. But what is the status of this princi- 
ple in quantum mechanics? If it is not valid in quantum 
mechanics, maybe we do not have the right to couple the 
quantum zero point fluctuations in the standard way in 
the Einstein equations? Could it be a way to avoid the 



problem? In order to discuss this issue, we will consider 
two situations, the Collela, Overhausser and Werner ex- 
periment (COW) and the quantum Galileo experiment. 

A. The Collela, Overhausser and Werner (COW) 
Experiment 

The COW experiment [133-135] consists in sending 
neutrons in an interferometer placed in a weak gravi- 
tational field, see Fig. 10. The presence of the gravi- 
tational field affects the wave-function of the neutrons 
and changes the relative phases between the two beams. 
When the interferometer is rotated, the gravitational 
phase shift is modified because the height difference be- 
tween the two arms is changed. This effect can then be 
detected (and was detected) in the interference pattern. 
The COW experiment is conceptually very important. 
This was the first experiment to measure the effect of 
the gravitational field on the wave-function. In fact, it 
shows that the gravitation is just an ordinary force as far 
as the Schrddingcr equation is concerned. 

Let us now describe the experiment in more detail. The 
neutron interferometer is represented in Fig. 10. Let us 
assume that, along one of the two paths, a phase shift Aip 
is introduced. Our first goal is to evaluate the intensity 
seen by the detectors O2 and O3. After diffraction, the 
wave function of the neutron takes the form 

(r) - <PT{9)c lk ° r + <S>D{d) C - ikG r , (779) 

where fco is the incident wave-number and fee the wave 
number after Bragg diffraction (which, therefore, must 
satisfy the Bragg condition). The quantity <E> is the am- 
plitude of the incident wave and T = T{9) (for "trans- 
mitted") and D = D{9) (for "diffracted") are two co- 
efficients describing the amplitude of the two branches 
of the wave- function after the neutron has emerged from 
the crystal. Let us now consider the the path ABDO3. 
At point -B, the wave function is given by (the origin of 
the coordinates is chosen to be at point A) 

ty(AB) = 3>T{9)e ik ° rB , (780) 

since the wave was just transmitted. Then, at point B, 
we have a diffraction and, therefore, the wave-function at 
point D reads 

V(ABD) = <PT(9)c lkorB D{8)c ikG{r °- rB \ (781) 

Finally, at point D, we have another diffraction but this 
time with opposite angle. The leads to 

V(ABD0 3 ) = <S>T(e)e iko ' rB D(6)e iko < r °- rB) 

xD{-6)c ikG( - r - rv) . (782) 

We also assume that, along that path, there is a phase 
shift Atp that we do not specify for the moment. As a 
consequence, the total wave function can be written as 

V(ABD0 3 ) = <S>T(6)c ik °- rB D{9)c ik °- {rD - rB) 

xD(-9)e iA<p e iko < r - rD \ (783) 
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FIG. 10: Representation of the Collela, Overhausser and Werner (COW) interferometer. The interferometer is placed in a 
vertical gravitational field (represented by the green arrow) and is rotated around the axis AB. 



where we have introduced the factor e lAip . gives the expression of the wave function at detector O3, 

It is straightforward to follow the same procedure namely ^(Os) = ^(ABDOz) + ^(ACDOj,), 
and to establish the expression for ^(ACDOs). This 



*(0 3 ) = $ \T(6)e lko rB D(e)c lko irD - rB) D{-6)c lAlp + D(9)c lk ° rc D{-9)e lka(r °- rc) T{9) 



J,k -{i — t-d) 



(784) 



Using this expression, one can calculate the expected intensity. Straightforward manipulations leads to the following 
expression 



/(O3) = |*(0 3 )| 2 =a(l + cosA^) 



(785) 



where a is a constant that can be determined from the above expression and which depends on the incident flux and 
the crystal structure. Then, one can repeat the same analysis and determine the wave function at detector 2 - It 
reads 



*(0 2 ) = $ T(9)e lkorB D< y 9)e lkG ^ D - rB '>T(-e)c lA ' f ' + D(9)e lkcrc D(-9)e iko < rB - rc) D{9) 

I 



a ifca-(i — r D ) 



(786) 



The intensity at detector O2 can be computed along the 
same lines and one is led to 

I(0 2 ) = /3 -a cos Aip, (787) 

where a and /3 are just constants that do not play an 
important role in what follows [of course, the constant a 
is the same as in Eq. (785)]. 

Having estimated the form of the signal expected for 
any phase shift Aip, we now turn to its calculation in 
the case of the COW experiment. For this purpose, we 



return to the Schrodingcr equation which can be written 
as 

V 2 * + w 2 1' = 0, (788) 

with 

- 2 = I? ( E V ) ■ ( 789 ) 
Now, let us assume that the potential is given by 

V = V + AV(r), (790) 
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where Vo is a constant and AV a small space-dependent 
small perturbation, i.e. AV <C Vq. If the perturbation 
does not change abruptly, the Wentzel-Kramers-Brillouin 
(WKB) wave- function 



^ ex exp 



i / w(r)dr 



(791) 



is a very good approximation to the actual solution. In- 
troducing the notation fc = 2m(E — Vo)/h 2 and using 
that p = mv = hko, an expansion of the phase in AV/V) 
leads to the following expression 

* oc exp (ik -f-^J AVdtj , (792) 

where the subscript "P" indicates that one integrates 
along the non perturbed path. From the above equation, 
we conclude that 



1 

Atp = -- 



AVdt. 



(793) 



The same result can also be expressed differently. Using 
energy conservation, on can write 



P 

2m 



Vo + AV = E, 



and, since p = po + Sp, one obtains 
Po 



2m m 



5p + V + AV = E, 



(794) 



(795) 



or. 



v ■ Sp = -AV. 



(796) 



As a consequence, the phase shift can be re-written as 
the following expression 



A(p 



1 



5p ■ dr, 



(797) 



where the integral should again be calculated along the 
path of the neutron. 

Let us now determine the phase shift for the COW ex- 
periment. Let us consider two point Pi and Pi that corre- 
spond to one beam transmitted and one beam diffracted 
in the interferometer, see Fig. 11. The potential differ- 
ence between those two points is given by 



V(P 2 ) - V(Pi) = mgdsmijj, 



(798) 



where d can be expressed as d = Apcos(7r/2 — 26). 
Given that the quantity x in Fig. 11 can be written as 
x = AP\ cos 9, one deduces that d = 2xsin#. As a con- 
sequence, 



V(P 2 ) — V(Pi) = 2m g rav<7a;sin0sin^>, 



(799) 



where we have now carefully written that the mass which 
appeared in the above expression is the gravitational 
mass (since this is the "coupling constant" to the gravi- 
tational field). This implies that the phase shift can be 
written as 

2 f L dx 
Aip = — — I 2m g gox sin simp -, (800) 



ri 
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where we have taken into account the fact that the line 
P\Pi is "moving" with a speed vcosB. The factor two 
in front of the whole expression comes from the fact that 
the phase shift accumulates and is, along the total path, 
twice its value in the fist part of the interferometer. The 
integration over the variable x is easily done and we ob- 
tain 



hv 



L 2 tan 9 sin tjj. 



(801) 



Finally, noticing that the total area of the interferometer 
is nothing but A = 2L 2 tan 8 and that the de Broglie 
relation is v = h/(rrii n X), where m-m is the inertial mass 
and A the wavelength of the incident beam, we arrive 
at [133-135] 

Aip = — ■j^m sia , v m in goAXsm if). (802) 

This expression is quite remarkable as it involves the 
gravitational and the inertial mass together with the 
Planck constant. Moreover, we deal with an observable 
quantity which docs not depend on the ratio of the grav- 
itational mass to the inertial mass, as is usually the case, 
but on the product of these quantities. We also see that 
the intensity of the signal at detectors O2 and/or O3 de- 
pend on the orientation of the interferometer if). 

However, in practice, one must deal with two other 
effects [133-135]. The first one is the Sagnac effect due 
to the rotation of Earth. Classically, the Hamiltonian of 
a neutron is given by 



H = 



2m; 



Wgravfl 1 ' r - OJ • L, 



(803) 



where L = r x p is the angular momentum of the neutron 
and u> the angular velocity of Earth. Upon using the 
Hamilton equations, this implies that 



p = m;n r + m m oj x r. 



(804) 



Therefore, one obtains a Sagnac phase shift which can be 
expressed as 



— p <p (w x r) ■ dr, (805) 



where we have used Eq. (797). Using Stokes theorem, 
this can also be written as 



47rm in 
Aip„ = — - — u) ■ A., 

T Sagnac 1 ' 



(806) 



where A is the vector associated to the interferometer 
area. If T is the incident neutron beam west of due south 
and S the co-latitude, the Sagnac phase shift takes the 
form 



A^ S a gna 



■w A (cos ip cos S + sin ip sin T sin S) 



(807) 



The COW experiment is such that the incident beam is 
directed due south which implies that T = 0. In the case, 
one obtains 



Ap 



uj A cos 5 cos 



(808) 



a dependence in ip which is different from the gravita- 
tional phase shift, see Eq. (802). 

Finally, there is a shift due to the fact that the inter- 
ferometer bends and/or warps under its own weight. The 
COW experiment claim that this effect can be described 
by the following equation 



A( ^bo„d = -'7bc„d sin </'- 



(809) 



This expression can be justified by noticing that the 
bending effect depends on the rotation angle. Clearly, 
this effect is difficult to estimate from first principles. 

Therefore, the total phase shift expected in the COW 
experiment is the sum of the three contributions dis- 
cussed before, namely the gravity, Sagnac and bending 
shifts. This leads to the following expression 

A( P = -?grav Sil1 V> + <7 S ag„ac C0S ^ ~ 9bo„d Sm (810) 

where, in order to take into account Eqs. (802) and (808), 

^goAX/h 2 and 



-2i:m 



g rav 77ii n 



we have defined q grav e 

cognac = Airm^ojA cos S/H. The above expression can 
be re-written as 



where q 2 

<?Sag„ac/ (?g r 



Aip = qsm(ip - ip ) , 

(<W +9bc„d) 2 + ffLnac 



(811) 



and tan ipo 



<7 bond ). In practice, the Sagnac effect is 
only 2.5% of the gravitational effect but we see it leads 
to a global shift of the oscillatory pattern. 

Let us now describe the result of the experiment. Ac- 
cording to the previous considerations, the intensity ob- 
served at detector O3 is 



1(0 3) =a{l + cos [gsin(^ - Vo)]} ■ 



(812) 



By Fourier transforming this signal, one can extract the 
frequency q rav of this oscillations. If one repeat this pro- 
cedure for different wavelengths, one can fit the depen- 
dence of q rav with A which allows us to determine the 



product m g 



1. Let us recall that the neutron mass is 



usually obtained by mass spectroscopy on the deuteron 
which leads to the following value 

TP 

m n = run - m v + -1 = 1.6747 x 10~ 24 g, (813) 

where E 1 is deuteron binding energy. On the other hand, 
the COW experiment has found [133] 



(m grav m in ) 1/2 = (1.675 ± 0.003) x lO"^ 4 g. (814) 

This shows that the weak equivalence principle seems to 
be satisfied even in a purely quantum mechanical situ- 
ation. This also shows that gravity can be "coupled" 
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to the Schrodinger equation in a standard way. The 
COW experiment is therefore very important since this 
was the first time that this was demonstrated. More- 
over, if we consider that m n = m m then the two previous 
equations shows that the COW experiment proves that 
TOi n = m grav even in the quantum domain. 



B. The Quantum Galileo Experiment 

It is well-known that, classically, the motion of a body 
in a constant gravitational field only depends on the ratio 
TOgrav/win (which is measured to be unity). However, if 
one writes the Schrodinger equation for the same situa- 
tion, namely 



ih 



& &(t, z) 



H 2 d 2 V{t,; 

dz 2 



2m; 



m grav gz*(t, z), (815) 



it is apparent that m- ln and m grav no longer cancel out. 
Therefore, it is interesting to study how the universality 
of the free fall is recovered in quantum mechanics. We 
investigate this point in the following sections. 



1. The Salecker-Wigner- Peres Clock 

In this section, we present a simple model of a quantum 
clock, the so-called "Salecker-Wigner-Peres clock" [136- 
140] . We introduce this system because we want to study 
how a quantum particle falls down in a constant gravita- 
tional field. One way to to do that is to measure the time 
of flight of the particle and to compare it with its clas- 
sical analogue. Hence the question of how to define the 
time of flight becomes important. As is well-known, the 
measurement of time (or, to be more accurate, of time of 
arrival) is problematic in quantum mechanics [141-145]. 
In particular, there is no time operator because time is 
not a dynamical variable. One way out is to introduce a 
quantum system, the "clock" , that is coupled to the par- 
ticle and changes its state when the particle crosses the 
starting and the ending points along its "trajectory" , see 
Fig. 12. Reading these two states then gives information 
about the time of flight. We now study how this general 
idea works in practice. 

We assume that the clock has an odd number N = 
2j + 1 of states. The Hamiltonian of the clock is given 
by 



d 

#clock =UJJ = -ihw — , 

o9 



(816) 



where the quantity 9 can be viewed as describing the 
position of the clock's hand. The corresponding states of 
the clock are represented by the following wave-function 



U n (6) 



1 



J 1)1) 



'2tt 



(817) 



that are eigenfunction of -ff c iock since H c \ oc \ l u n (9) 
nhu>u n (6). It is also easy to show that 



(818) 



which gives the time evolution of the system. 

It is also convenient to introduce another basis v s (9) 
defined by (s = 0, ■ • • , N - 1) 



v s {6) 



— y 



c — —u n (6) (819) 

n=-j 

1 sm.[N(6-2-Ks/N)/2] 



sin [(6 - 2irs/N) /2] 



(820) 



If N is large enough, these functions peak at 9 = 2ns /N 
and this corresponds to the hand of the clock pointing to 
the s-th hour with an angle uncertainty of ±tt/N (hence, 
as expected, if N is large, the accuracy is good). More- 
over, using Eq. (818), we have 



- lH ^ t/h v s {9)=v s (9-ujt). 



(821) 



This means that if we evaluate the above equation for a 
time r such that ojt = 2n/N, then one obtains 

e-ur— r /% s (6)=v s+HmodN) (6), (822) 

which confirms the interpretation given before. 

Let us now consider a free particle traveling along the 
x-axis. We want to measure the time of flight between 
x = and x = L, see Fig. 12. In order for the clock to 
record this time, we must of course couple it to the parti- 
cle. Following Rcf. [137], we assume that the Hamiltonian 
of the total system "particle+clock" can be written as 



H 



P 

2m 



V(x)H clock , 



(823) 



where V(x) is the projector operator which is one if 
< x < L and zero otherwise and m is the mass of 
the particle. We write the eigen wave-functions of the 
system as 

*(t,x,0) = iP(x,9)c- lEt/h 



1 



r Mx>n(0)c- lEt/ \ (824) 

n=-j 

and inserting this form into the Schrodingcr equation, 
one obtains 



E 

n=~j 



h 2 dV» 

2m dx 2 



nhujV(x)ip n (x) 



u„(e) = o. 



(825) 



Therefore, we just have to solve the time-independent 
Schrodinger equation in presence of a rectangular bar- 
rier. Of course, this exercise is solved in any textbooks 
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x=0 x=L 



FIG. 12: Sketch of a quantum "Peres" clock. The clock is coupled to the particle traveling from the point x — to the point 
x — L. When the particle crosses x = 0, the clock changes its state and starts running. When the particle reaches x — L, the 
clock stops and a comparison of the two clock states permits a measurement of the time of flight of the particle. 



on quantum mechanics, see for instance Rcf. [57]. It is 
straightforward to show that 



Ae ikx ■ Be "" . a;<0, 
Ce lpx + De~ tpx , 0<x<L, 
Ee lkx , x > L, 



-ikx 



where we have defined 



s/2mE 



\/2m{E — nftju) 



= kVT 



(826) 



(827) 



(828) 



with e = nhw/E, a quantity which measures how the 
clock perturbs the particle. The quantities A, B, C, D 
and E are integration constants that are fixed by requir- 
ing the wave-function and its derivative to be continuous 
at the barrier. Notice that, after the barrier there is only 
one branch; this is of course because we assume that the 
particle is coming from the left, see Fig. 12. Matching 
ip n and dipn/dx at x = and x = L, one obtains 



.4 
B 



1 

4p 
1 

4p 



(k + p)S k - p 1 L - (k - p)S k+ ^ L 
(k + P )S k -^ L — {k — p)S k+ ^ L 



1 

4k 
1 

4^ 



(k + p)c l( - k - p)L +(k - p)S k+ ^ L 
(k + p)S k ~^ L + (k - p)S k+ ^ L 



(829) 
(830) 



Then, we assume that the clock does not perturb the 
particle too much and, as a consequence, that e«l. In 
this case, B ~ and A ~ c l ( k ~p) L E. Indeed, we have 
p = ky/1 — e which, at leading order, simply gives k ~ p 
(but, of course, one should use the next to leading order 
when one estimates the argument of the exponentials, 
see below). Then, initially, this means that the wave- 
function can be written as, see Eq. (824) 

-, n=j 

tp(x,6) = -y= Ac lkx u n {9) = Ac lkx v {9), x < 0, 



n=-j 



(831) 

because k docs not depend on n, see Eq. (827). We see 



that, initially, the clock points to zeroth hour, as ex- 
pected. After the barrier, using the previous considera- 
tions, the expression of the wave-function can be written 
as, see Eq. (824) 



ip(x,e) ~ —e lkx e- z{k - p)L u n (9), x » L, 

n=-j 



where we have used Eq. (819). But we have 



(832) 



(833) 
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which implies that (for x 3> L) 



A 



n=j 



^ikx ^ ^ ^inO — inuL/ (2E f m) 1 ^g^^J 



Ae ikx v Q 



ujL 



y/2E/r 



(835) 



This means, see Eq. (821), that the clock now indicates 
the time 



t 



L 



L 



{2E/m) 1 / 2 v ( 



las 



(836) 



where i> c i as = {2E/m) 1 / 2 is the classical velocity. There- 
fore, the clock has measured a time of flight which is 
nothing but the classical time of flight L/v c \ as [137]. 

We now consider the situation represented in Fig. 13, 
see Ref. [140]. An observer, located at x = —L, sends a 
particle towards a reflecting barrier of height Vq located 
at x = 0. When the particle leaves the Peres clock is 
turned on. Then the particle bounces back and when it 
comes back at the observer position, the Peres clock is 
turned off. With this experimental set up, one can mea- 
sure the out-and-back time of flight. Following the same 
approach as before, the wave-function can be written as 



where we now have 
k 




Be~ lkx , x < 0, 
x > 0, 



V 



v /2m(E - nhuj) _ \/2mE 

h ~~ h 

y/2m(Vo + nhoj - E) 



E 



y/2m(V - E) f 

n V ' v -E 



1/2 



(837) 



(838) 



(839) 



Then, we have to match the wave-function and its deriva- 
tive at x = 0. This gives 



D 

G 



k — ip 
k + ip 

2k 
k + ip 



A = RA, 
A. 



(840) 
(841) 



As before, we assume that, initially, the clock was is a 
given state, here vq(9). Therefore, the wave- function in 
the region x < can be expressed as, see Eq. (824) 



A — J 

ip(x, 9 ) = -j= ( e * kX + Re~ %kx ) u n {0)- (842) 



But = 1 and its phase is given by —2 arctan(p/fc). As 
a consequence, one has 



i/j(x=-L,e) = 



n—j i- 

— y 

n=-j 



—ikL 



| g— 1% arctan(p/fc)gifcL 



u n (9). (843) 



The next step is similar to the case of a free particle, 
namely we expand the phases of the two wave-function 
branches in terms of the parameter e. This gives 



4>{x = -L,0) 



A 

+ in | 6 - 



]T cxp 

n=-j 

ujL 



iL 



V2jnE + in [0 + 

n V v clas 



exp 



Vq-E\ iL r—— 

— e — / ^ ~K 



fclas y/E(V - E) 



(844) 



A iL , \ I ujL \ A 

—= cxp — — v2mE ) v \ -\ H — = exp 



-2i arctan 



Vo - E \ iL 



— V2mE 
n 



luL 



XVq 



fclas y/E(V - E) 



(845) 



and this allows us to directly read the time indicated by the clock. For the first branch, it is —L/v c \ as while for 
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x=-L 



x=0 




FIG. 13: Measurement of an out-and-back time of flight. An observer, located at x = — L sends a particle towards a square 
potential step located at x = 0. The physicist turns on the clock when the particle leaves and turn it off when the reflected 
wave (i.e. the reflected particle) comes back and reaches x = — L. 



the second one, it is given by L/v c \ &s + hj yj E(Vq — E). 
As a consequence, the time of flight of the particle can 
be expressed as 

At = — + - H (846) 

Vclas y/E(V ~ E) 

Therefore, the time of flight is the classical one, 2L/v c \ as , 
but, this time, there is an additional contribution. This 
one can be understood as follows. The quantity At can 
be re-written as [140] 

2{L + d) 



At 



(847) 



where d — 1 /p, see Eq. (839) (evaluated at e = 0) is the 
penetration depth into the potential step. This means 
that there is an additional delay due to the tunnel ef- 
fect and the fact that the particle has a non vanishing 
probability to be below the barrier [140]. In this regime, 
the "velocity" is nevertheless given by its classical value 
v clas [140]. 

This concludes this section on the Peres clock. We 
have seen that this is a useful device to measure times of 
flight. We now use it to study the motion of a quantum 
particle in an uniform gravitational field. 



2. The Gravitational Case 

In this section, we study the status of the weak equiv- 
alence principle and of the universality of the free fall in 
quantum mechanics [139, 140, 146-150]. 

Following Ref. [140], let us consider the situation where 
a particle is sent upwards in a uniform gravitational field, 
see Fig. 14. If z denotes the vertical coordinate, the corre- 
sponding classical situation is described by the following 
Lagrangian 

^in . 9 
Z — 1 

2 



L(z,z) 



where g is the gravitational field and where we have made 
the difference between the inertial mass mi n and the grav- 
itational mass (the "gravitational charge") m g . The so- 
lution to the classical equation of motion reads 



z(t) 



1 TO grav ^ 

2 m- m 



gt + v ini t + L, 



(849) 



where z — L is the initial altitude (at t = 0) and v- nli is 
the initial velocity, related to the energy E of the particle 
by the following expression 



2ff- 



E 



%rav<? 



L 



1/2 



(850) 



Of course, the motion of the particle only depends on the 
ratio TOgrav/wm. The maximum altitude, z max , reached 
by the particle, see Fig. 14, is given by 

(851) 

,9 



E 



which means that, from its initial position to z max , the 
particle can rise a distance of E/(m glav g) — L. After 
having reached z max , the particle falls back and the up- 
to-down classical time of flight reads 



At 



TTlty 



1/2 



E 



- L 



''''gravy 



1/2 



(852) 



If TOg ra v = ?7i; n , this equation can also be written as At = 
2v ini /g. 

Then, we would like to consider the same situation 
but from a quantum-mechanical point of view [140]. In 
particular, we would like to measure the quantum time 
of flight. Given the discussion of the previous section, 
it is clear that we must couple the quantum particle to 
a Peres clock located at z = L. Then, the next step is, 
of course, to solve the Schrodingcr equation. Following 
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FIG. 14: Sketch of a quantum particle in an uniform gravitational field. The particle is sent upward, the initial time being 
recorded by a quantum Peres clock. The particle reaches a maximum altitude and then falls back to the detector at z = L 
where the Peres clock measures this return time. 



the previous considerations, we just have to consider the 
equivalent of Eq. (825), but for the case of an uniform 
gravitational field. It reads 



dz 2 



2mj n (E — nhw) 2m- u 



■f9 



r, 2 



r, 2 



i> n = o. 



(853) 

It is straightforward to obtain that the solution can be 
expressed in terms of Airy functions [48, 49] Ai(z) and 
Bi(z). If we retain only the branch that leads to a 
bounded wave-function, the result is given by 



71= J 



<P(z,9) = ^= V Ai 

n=-j 



-b 



u n {9), 



(854) 



where A is constant fixed by the normalization of the 
wave-function and a and b are defined by 



2m in m giav g 
E — nfko 



1/3 



(855) 
(856) 



The next step consists in evaluating the wave-function 
at x = L, with the assumption that the return point is 
located far from the maximum altitude. This means that 
L — b is a negative quantity and that \L — b\ is large [more 
precisely, one has L — 6 = L — E / (m glav g) + Ee/(m glav g) 
and the quantity L—E/ (m grav g) is minus the distance the 
particle can rise, see the remark after Eq. (851). There- 
fore, L — b < and, if this distance is large, \L — b\ is 



indeed a large quantity]. In this limit, one has 




u n {6), (857) 



where £ = (2/3) [(6 - i)/a] 3 / 2 . Explicitly, this quantity 
can be expressed as 



3a 3 / 2 
2 



E 



E 



(l-c)-L 

3/2 



3/2 



(858) 



3a 3 / 2 \ 



E 



a 3 / 2 \m glav g 



- L 



1/2 



(859) 



In order to calculate the difference between the two times 
indicated by the clock and corresponding to the two 
branches of the wave-function, it is in fact sufficient to 
estimate the part of the phase shift that depends on 
n. Therefore, it is clear that the phase shift ±7r/4 in 
the expression (857) will not contribute because it is n- 
independent. Moreover, the term [(z — 6)/a] -1 / 4 is n- 
dependent but is common to the two branches of the 
wave-functions. As a consequence, it will cancel out when 
the phase shift is determined (but would participate to 
the "absolute" time indicated by the clock when the par- 
ticle leaves z = L and returns to z — L). Summarizing, 
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the up-and down time of flight is given by 




that is to say exactly Eq. (852), see Ref. [140]. Sev- 
eral remarks are in order here. Firstly, although it is of 
quantum-mechanical origin, the above equation does not 
contain h that cancels out. Secondly, this result seems 
to indicate that the weak equivalence principle holds at 
the quantum level. Thirdly, it may appear surprising 
that the quantum At is not corrected by a term corre- 



sponding to the penetration depth, see Sec. XIV Bl, in 
particular Eq. (847), as it is the case for a square poten- 
tial. The interpretation given in Ref. [140] is that there 
is a finite probability that the particle tunnels below the 
barrier but there is also a finite probability that the parti- 
cle is scattered back before the classical turning point. If 
these two probabilities cancel then one obtains the clas- 
sical prediction. Fourthly, the above result is valid only 
if the return time is measured far from the maximum 
altitude. What happens if, on the contrary, the particle 
starts from a position which is close to the turning point? 
To investigate this case, we re-start from Eq. (854) and 
write the solution as 
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N ^ 3 I 

n=-j 



^ /3 Jl/3 (0 + C-^ /3 J-l/3 (C)] + [(l - C W3 ) Jl/3 (0 + (l - e- M/3 ) J-l/3 (C)] } MO), 



(862) 

where y = (b — z)/a > and C = 2y 3 / 2 /3. In this way, we have identified the incident (first term in the squared 
bracket) and the reflected wave (second term in the square bracket). Indeed, using the asymptotic behavior of the 
Bessel functions at infinity, one can check that each branches give precisely e ± ^ +i7r / 4 - ) if z = L, see Eq. (857). If we 
now consider the limit y <C 1, then one arrives at 
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the two first terms corresponding to the incident wave 
and the two last ones to the reflected wave. The two 
branches are of course complex conjugate to each other. 
Therefore, if we denotes by H the complex number cor- 
responding to the two first terms in the above equation, 
then the wave-function takes the form 
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argH = arctan 
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The rest of the calculation proceeds as before. We first 
expand the phase in e and, then, take the limit E/(m g g)— 
L — >• 0. This leads to 
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As a consequence, one finds that the time of flight is given 



by 



At 
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(867) 
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in rough agreement with Ref. [140]. This time, we no- 
tice that h has not canceled in the final expression and 
that the time of flight is different from its classical coun- 
terpart (852). This also means, as already mentioned, 
that the cancellation which leads to a quantum time of 
flight equals to the classical one is in fact valid only if the 
particle starts far from the turning point. 

We conclude that the weak equivalence principle (or 
the universality of the free fall) can be extended to the 
quantum regime. The linear potential seems to possess 
precisely the shape which leads to identical times of flight. 
It is important to realize that this property depends on 
the shape of the potential and would not be obtained with 
another shape, see Ref. [140]. Moreover, this seems to be 
true only far from the turning point. If one identifies the 



83 



wave-vector to k — {2m- va m S[ajV g / S 2 ) 1 / 3 , see Ref. [151], 
the de Broglic relation v ~ hk/m allows us to estimate 
a velocity. Then, from the quantum time of flight, see 
Eq. (868), we can construct a length ~ [fi 2 /(m 2 g)] 1 / 3 
which typically controls the validity of the approximation 
just mentioned. 



3. Atom Trampoline 

In this section, we want to briefly mention the case of 
the "atom trampoline" [150, 151]. Let us consider again 
the experimental set up of the previous section but let 
us now assume that a reflecting wall has been installed 
on the ground, at z = 0. Then, the boundary condition 
ip n (0) = leads to Ai(— b/a) = 0. If z n denotes the 
zero of the Airy function, this means that we now have 
a discrete spectrum of energy levels given by 
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defined according to 

*(t,z) = <p[t'(t),z'(t,z)] 
x exp 
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satisfies the equation [146] 



_d<p{t',z') h 2 d 2 <p{t',z') 
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2m dz 1 ' 



+ mgz'ip(t',z'). (874) 



The potential in the above equation is precisely of the 
form describing a constant gravitational field. In this 
sense, we find that the equivalence between a constant 
gravitational field and an accelerated frame propagates 
to quantum mechanics. 



5. Falling Composite Quantum Objects 



It is interesting to notice that this result depends on 



the ratio 



,/mi n and not on 



' gray 



Moreover, 



this spectrum has been observed (with ultra-cold neu- 
trons which leads to E n of the order of ~ fO -12 cV) in 
the gravitational field of the Earth [152, 153]. There- 
fore, this result confirms the discussion about the COW 
experiment described in Sec. XIV A. It experimentally 
establishes that the gravitational force "couples" to the 
Schrodinger equation in a standard way. 



4- Schrodinger Equation in an Accelerated Frame 

It is also frequent to refer to the weak equivalence prin- 
ciple as the property stating that, locally, the effect of a 
constant gravitational field can be mimicked by an ac- 
celerating frame. Therefore, it is interesting to study 
whether this claim holds in quantum mechanics. This 
question has been studied in Refs. [146, 150, 151]. Here, 
we follow the treatment of Ref. [146]. 

Let us consider the free Schrdndingcr equation (in one 
dimension to simplify the problem). It reads 
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Then, let us consider the following coordinate transfor- 
mation 



z = z 
t' = t 
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It is straightforward to show that the function (p(t',z') 



In the previous sections, we have studied how a quan- 
tum particle falls in a uniform gravitational field. For in- 
stance, we have investigated under which circumstances 
one can say that the weak equivalence principle holds in 
quantum mechanics. Here, we would like to study the 
same question but in the case of a composite quantum 
object. In particular, it is interesting to ask whether the 
motion of the composite system (the atom) can depend 
on its internal state. 

Let us start by recalling the Ehrenfest theorem [57] for 
a quantum particle. It states that the time derivatives of 
the mean position and momentum are given by (for sim- 
plicity, we consider the case of a one-dimensional system 
since we have in mind a particle falling down along the 
z-axis) 
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In the case of a linear potential, V(z) = m gTa , v gz, the two 
previous equations expressing the Ehrenfest theorem can 
be combined in such a way that we arrive at 



d 2 (z) 
dt 2 
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(877) 



This means that the average position (z)(t) evolves ex- 
actly as the classical z given in Eq. (849). However, let 
us also notice that the variance of the position is not a 
function of the ratio m grav /mi n only, see in particular 
Refs. [146-148]. 

After this short reminder, let us consider a bound sys- 
tem of two particles interacting through a potential V 
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which only depends on the distance between them. We 
also assume that this system is embedded in a region 
where there is an exterior field described by the poten- 
tial U(r). The corresponding Lagrangian reads 



t l ■ ■ \ m in.l -2 , m in,2 . 2 

L{ri,r 2 ,ri,r 2 ) = -^^r 



2 Jl ~ r ~ 
-U( n )-U(r 2 ) 



r 2 - V (n - r 2 ) 
(878) 



where r\ and r 2 are the positions of the two particles and 
77Mn,ij "^in,2 their mass. If we now introduce the center 
of mass position vector, 



m-in.l 



-n 



Win, 2 



min.i + m ini2 m in4 + m illi2 
and the relative position vector, 

r = n - r 2 , 
then Eq. (878) takes the form 
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£(r OJ r,r Q ,r) - yfj + ~ r 2 - V(r) 

-U( ri )-U(r 2 ), (881) 

where M = m m j + mj n 2 is the total mass and /i = 
Win,iTOi ni 2/(T7ii ni i + mj nj2 ) is the reduced one. Of course, 
in the above expression r± and r 2 must be viewed as 
function of r G and r 5 . From the above considerations, 
it is straightforward to obtain the Hamiltonian of the 
system. We arrive at the following expression 



H = ^ + ^ + V(r) + U(r l) + U(r 2 ), 



where p G = P\+p 2 and p = (/i/m ina )pi - {fi/m in>2 )p 2 . 
If we now assume that the external field is a uniform 
gravitational field (directed along the z-axis), then it is 
easy to show that 

U(r\) + U(r 2 ) = (TOgrav.l + "lgrav,2) 5 Z G 

1 / 
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U(z G ) + U(z), 



where, according to Eq. (880), z = z\ — z 2 . As a conse- 
quence, the Hamiltonian (884) takes the form 



H = H„ + H 



(887) 



5 The relation between these quantities is given by 
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where the center of mass and relative Hamiltonians can 
be expressed as 



P — 

2AI V gJ 

2// 



+ V(r) + U(z). 



(889) 



The physical interpretation of these equations is clear. 
The center of mass behaves as a free particle embed- 
ded in a uniform gravitational field. On the other hand, 
the internal states of the composite system are deter- 
mined by the potential V(r) but are also affected by the 
gravitational field through the term U(z). In practice, 
one expects that U(z) <C V(r) and, therefore, the per- 
turbations of the energy levels will be either very small 
or even totally negligible. Moreover, because [z,p Gi2 ] = 
[z G ,p 2 ] = 0, the two above Hamiltonians commute 



[H a ,H ti 



0, 
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which implies that the total state space is in fact the 
tensorial product of the center of mass state space and 
of the relative state space. 

Let us now study the behavior of (z a ). One has 
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In the same manner, one can now calculate the evolution 
of (p a ,z)- One obtains 
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Therefore, combining Eqs. (891) and (892), we arrive at 
the following expression 



d 2 (z G ) m 

grav.l ^grav,2 



At 2 
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Of course, this equation is similar to Eq. (877) and this 
implies that the center of mass mean value of the quan- 
tum composite system will follow the classical trajectory. 
In particular, this means that (z G ) does not depend in 
which internal energy state the system is placed. There- 
fore, if, for instance, we consider a falling Hydrogen atom, 
the behavior of the position mean value of the atom will 
not depend of whether the atom is in the state, say, Is 
or 2p. We also notice that the coefficient in front of g in 
Eq. (893) is the ratio of the total gravitational mass to 
the total inertial mass. In this sense, the weak equiva- 
lence principle is also satisfied by the composite quantum 
systems. 
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Finally, it is also interesting to repeat the gcdankcn 
experiment in Fig. 14 but with a composite system. It 
seems reasonable to couple the Peres clock to the center 
of mass of the system which means that the Hamiltonian 
of the system is still given by 



H = H r 



H 



(894) 



where, now, the center of mass and relative Hamiltonians 
can be expressed as 
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Then, following Eq. (824), we write the wave- function of 
the system as 



n=j 



-iEt/h 
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(897) 

The next step consists in inserting this wave-function in 
the Schrodinger equation in order to obtain the equiva- 
lent of Eq. (825). The result can be written as 



-—A G ip n (r G ) + U(z G )^ n (r G ) + nhwV(z G )t})„{r G ) - E G ^ n (r G ) 



+ ^n(r G ) 
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(898) 



where E = E G + E tc1 . This implies that only the center 
of mass part, ip n (r G ), of the total wave- function will de- 
pend on n. As a consequence, the internal part should in 
fact be written as x( r ) instead of Xn(f)- This means that 
x{ r ) cannot cause a rotation of the Peres clock. There- 
fore, the time indicated by the clock will not depend on 
the internal state of the composite system. We conclude 
that the time of flight of an atom, defined and calculated 
in terms of the Peres clock, does not depend in which 
quantum state this atom is placed. 

This concludes the section on the weak equivalence 
principle in quantum mechanics. The result is mixed. We 
have shown that, very often, it is satisfied in the frame- 
work of quantum mechanics. But, on the other hand, 
it seems sometimes to be modified by quantum effects, 
see in particular Eq. (868). Maybe this indicates that 
its application to the more complicated case of quantum 
vacuum fluctuations of a field is rather suspicious? 



XV. CONCLUSIONS 

We are now in a position where we can really formulate 
the cosmological constant problem. Vacuum fluctuations 
seem to exist in Nature and to have normal gravitational 
properties. Therefore, it seems natural to postulate that 
they participate in the value of the cosmological constant. 
However, when one tries to use well controlled techniques 
of quantum field theory to compute their energy density, 
one obtains a number which seems to be in contradiction 
with the measurement of A that one obtains in cosmol- 



ogy and with the constraints that one deduces from (for 
instance) the motion of the planets in our solar system. 
The difficulty of the problem stems from the fact that, 
in order to remove the above contradiction, one neces- 
sarily has to abandon something that is considered as ro- 
bust, i.e. renormalization in quantum field theory, weak 
equivalence principle, high accuracy measurements of the 
expansion of the universe, etc .... 

Even if this is not the purpose of this article, let us 
conclude this review with a few words about the solu- 
tions that have been proposed to solve the cosmological 
constant problem, for a complete overview of the subject, 
see for instance Ref. [154]. Of course, the most obvious 
solution is super-symmetry [19, 74, 75]. However, as we 
have discussed in Sec. VII D, super-symmetry has to be 
broken and this destroys the "miraculous" cancellation 
of the various terms participating the vacuum energy. In 
fact, the only common point between these various con- 
tributions is gravity. Therefore, this suggests that the 
mechanism that cancels the cosmological constant makes 
use of gravity. This idea has been explored in the so- 
called adjustment mechanisms [155-157]. These models 
work reasonably well but, unfortunately, they suffer from 
a severe disease. Indeed, in these scenarios, the New- 
ton constant is time-dependent and it turns out that the 
variation of the gravitational coupling is too strong to be 
compatible with the known constraints on \G/G\. There- 
fore, although the idea is very attractive at first sight, it 
seems that it does not yet exist a realistic realization of 
it. It is also worth signaling that back-reaction mecha- 
nisms have also been studied in Refs. [158, 159]. Here, 
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the idea is that the long wavelength (supcr-Hubble) cos- 
mological perturbations are described by a stress-energy 
tensor that exactly cancels the cosmological constant. 

Another class of solution is based on quantum grav- 
ity and quantum cosmology [160-162]. It is argued that 
the no-boundary Hartlc-Hawking wave-function is such 
that it peaks at A = 0. However, this approach suffers 
from various limitations among which is the fact that the 
path integral is not properly defined and that probabil- 
ities are not positive definite (this comes from the fact 
that the Wheeler-De Witt equation is in fact similar to a 
Klein-Gordon equation. The positivity of probabilities in 
quantum cosmology can be only be restored in the WKB 
approximation) . 

Recently, following the pioneered work of Ref. [163], 
there have been many attempts to solve the cosmological 
constant problem in the framework of extra dimensions, 
see chapter seven of Ref. [154] or Ref. [164]. Another 
popular solution, also based on high-energy physics con- 
siderations, is the landscape approach [12] in string the- 
ory [165]. This approach assumes that all the pocket 
universes of the landscape are populated during eternal 
inflation. Combining this fact with the anthropic princi- 
ple [166], it is argued that the most probable value of A 
is approximately the value observed today. This solution 
suffers from the difficulty of defining a measure on the 
landscape [167] and the use of the anthropic principle in 
this context has been criticized in Refs. [168, 169]. 

Finally, another possibility that has recently been in- 
vestigated is to modify the Einstein equations such that 
they become blind to a stress energy tensor of the form 
T^v oc g^u, see Refs. [170-172]. 

In order to be complete, let us mention that models 
of analogue gravity have also been developed to study 
the vacuum energy problem [173-175]. In these models, 
one can calculate the equivalent of the cosmological con- 
stant and study precisely its origin since the physics con- 
trolling the microscopic constituents is explicitly known. 
These models are based on many different types of ana- 
logue physical systems such as condensed matter sys- 
tems [176], super-fluid Helium [177], Fermi liquid [178], 



Bose-Einstein condensates [179, 180] etc ... Then, one 
can show that a simple calculation of the ground state 
may lead to an incorrect result, the only possibility in 
order to reach the correct answer being to fully take into 
account the microscopic theory. This would indicate that 
a correct calculation of the cosmological constant must 
necessarily be based on the fine structure of space-time, 
i.e. probably on a consistent theory of quantum gravity. 

As the variety of the subjects studied is this review 
shows, the cosmological constant problem is a very rich 
one. Given its difficulty, it is probable that it will have 
to wait for a new theory beyond the current standard 
model to find its resolution. This theory will have to 
describe the gravitational properties of the vacuum fluc- 
tuations and to regulate the ultra-violet infinities that 
appear in the calculation of vacuum energy density. This 
list of requirements almost constitutes an identity card 
for a theory of quantum gravity. But the most interest- 
ing aspect maybe goes the other way round. Given the 
fact that cosmology has enabled us to grasp observational 
signatures related to the vacuum energy, maybe these ex- 
perimental results will help us to deduce and to establish 
a convincing theory of quantum gravity? This is the rea- 
son why the cosmological constant problem appears to 
be so important and so interesting. 
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